Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

The Atacama Cosmology Telescope: A CMB lensing mass map over 2100 square degrees of sky and its cross-correlation with BOSS-CMASS galaxies

We construct cosmic microwave background lensing mass maps using data from the 2014 and 2015 seasons of observations with the Atacama Cosmology Telescope (ACT). These maps cover 2100 square degrees of sky and overlap with a wide variety of optical surveys. The maps are signal dominated on large scales and have fidelity such that their correlation with the cosmic infrared background is clearly visible by eye. We also create lensing maps with thermal Sunyaev-Zel’dovich contamination removed using a novel cleaning procedure that only slightly degrades the lensing signal-to-noise ratio. The cross-spectrum between the cleaned lensing map and the BOSS CMASS galaxy sample is detected at 10-σ significance, with an amplitude of A = 1.02±0.10 relative to the Planck best-fit LCDM cosmological model with fiducial linear galaxy bias. Our measurement lays the foundation for lensing cross-correlation science with current ACT data and beyond.

gravitational lensing↗
NASA NTRS2021

Positional Dependence of SNPP VIIRS Solar Diffuser BRDF Change Factor: An Empirical Approach

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership satellite regularly performs on-orbit radiometric calibration of its reflective solar bands (RSBs), primarily through observations of an onboard sunlit solar diffuser (SD). The on-orbit change of the SD bidirectional reflectance distribution function (BRDF) value, quantified by a numerical factor called the H-factor, is determined by the onboard SD stability monitor. Our previous study showed that the H-factor is solar angle and view direction dependent. In this study, we determine the dependence of the H-factor on the detector SD view footprint location. We fit an empirical model to the NASA Collection 1 SNPP VIIRS Level 1B (L1B) spectral reflectance difference across the detectors in an RSB over uniform Earth scenes of the Libya 4 desert and deep convective clouds. We apply the model predicted SD positional dependent H-factor to calibrate the RSBs. Under this new calibration scheme, the original unreal striping is removed from the homogeneous Libya 4 desert and the deep convective cloud images, as well as the original unreal striping from the Dunhuang desert image. The SD positional dependent H-factor has been used to calculate the SNPP VIIRS RSB radiometric correction factor for the NASA Collection 2.0 SNPP VIIRS L1B products.

SNPP VIIRS↗
NASA NTRSDate not supplied

The Application of Remote Sensing and Machine Learning to Improve Early Warning Systems for Harmful Algal Events in the Highland Lake Chain, TX

Beginning in 2019, harmful algal events have caused canine deaths in both Lady Bird Lake and Lake Travis located near Austin, Texas. These two reservoirs are part of the larger Highland Lakes chain, managed by the City of Austin Department of Watershed Protection (COA DWP) and the Lower Colorado River Authority (LCRA), which fulfill municipal, commercial, and agricultural water demands. Given the recent increase in favorable environmental conditions for algal events in central Texas, NASA DEVELOP partnered with LCRA and COA DWP to improve monitoring and early detection of algal events, utilizing satellite remote sensing and machine learning. Spatially and temporally varied chlorophyll a concentrations, cyanobacteria detections, turbidity, and water surface temperature products are used as environmental proxies. Landsat 8 Operational Land Imager (OLI) and the Sentinel-2 MultiSpectral Instrument (MSI) data were used to provide a combined revisit time of up to ~3 days and < 30 m per pixel products. Chlorophyll a concentrations were estimated using a pre-trained Mixture Density Network, and cyanobacteria detection was accomplished using the Broad Wavelength Algae Index, which can differentiate algal blooms from algal proliferations (mats). In situ data were used to validate remotely sensed measurements and quantify uncertainties. Preliminary results show a good fit between the modeled output and in situ observations, suggesting that remote sensing data can be used to retrieve biogeochemical properties and/or inherent optical properties (IOPs) of water columns in these inland human-made lakes. Uncertainties were introduced from the sensitivity to atmospheric correction, inherent mismatch between satellite and sampling data, and a relatively lower signal-to-noise ratio over water. The resulting products enable near real-time monitoring of environmental proxies relevant to algal event presence in the Highland Lakes chain, and will ultimately support water management, decision making, and risk communication.

Shuyu Chang↗
NASA NTRSDate not supplied

Dust Composition of Comet 81P/Wild 2 From JWST Spectroscopy Compared to Stardust’s Fine-Grained Materials and Gems-Rich IDPS

The Stardust Mission returned samples from the coma of Jupiter Family comet 81P/Wild 2 for detailed laboratory analyses. Here we present and discuss the best-fit thermal dust model for the coma dust of comet 81P/Wild 2 as observed by JWST. Comet 81P was observed through JWST GO 018xx, using NIRSpec (2.9–5.3µm, λ/∆λ≈1000) and MRS IFU (4.9–28.1µm, λ/∆λ≈3000) on 2023-03-20 UT and 2023-03-24 UT, respectively, at a heliocentric distance of 1.85 au and JWST distance of 1.43 au (phase angle of 32 degrees). The dust coma of 81P as revealed by JWST offers a salient compliment to laboratory studies of Stardust samples that typically are bigger than 2 µm and up to 60 µm in size.

D H Wooden↗
NASA NTRS2010

Kriging as a Means of Improving WAAS Availability

The Wide Area Augmentation System (WAAS), an augmentation of the Global Positioning System (GPS), provides safe and reliable use of GPS signals for airline navigation over much of North America. Currently the largest source of positioning error in the system is signal delay caused by the ionosphere. To allow the user to take account of such error, WAAS computes and broadcasts ionospheric vertical delays at a set of regularly-spaced grid points. In addition, WAAS computes and broadcasts a safety-critical integrity bound at each ionospheric grid point (IGP) called the Grid Ionospheric Vertical Error (GIVE). GIVEs are constructed to be sufficiently large to protect the user against positioning error due to the presence of ionospheric irregularity. In the initial operating capability (IOC) of WAAS, the vertical delay estimate at each IGP is determined from a planar fit of neighboring slant delay measurements, projected to vertical using an obliquity factor specified by the standard thin-shell model of the ionosphere. In WAAS Follow-On (WFO) Release 3, however, the vertical delay will be estimated by an established, geo-statistical technique known as kriging. Compared to the planar fit model, the kriging model is found, in general, to match better the observed random structure of the vertical delay. This paper presents the kriging methodology that will be used to estimate the vertical delay and its uncertainty at each IGP, and it assesses the subsequent improvement in WAAS availability enabled by kriging.

GPS↗
NASA NTRS1976

Pressure-induced absorption by H2 in the atmospheres of Jupiter and Saturn

Observations of the S(1) line of the pressure-induced fundamental band of H2 in the spectra of Saturn and Jupiter are analyzed by comparing the observed line shape with predictions of both a reflecting-layer model and a homogeneous-scattering model of the atmospheres. Upper and lower limits are derived for the values of relative intensity that occur between 5000 and 5100 kaysers, the percent absorption due to H2 at 5100 kaysers is estimated, and it is established that methane and ammonia cannot account for the broad absorption attributed to hydrogen. The comparison with the model atmospheres shows that the reflecting-layer model appears to give the best fit to the Saturn line profile for temperatures near 150 K, while both models provide good fits to the Jupiter data, but for widely differing temperatures. Difficulties encountered in determining the true continuum level, especially for Jupiter, are discussed.

Martin, T. Z.↗
NASA NTRS1986

Statistical analysis of failure data on controllers and SSME turbine blade failures

The expressions for the maximum likelihood functions are given when the failure data are censored at a given point or at multiple points, or when the data come in groups. Different models applicable to failure data are presented with their characteristics. A graphical method of distinguishing different models by using cumulative hazard fucnction is discussed. For the failure data on controllers the model is determined by cumulative hazard function and chi-square goodness of fit. Using the Weibull Model the maximum likelihood estimators of the shape parameter and the failure rate parameter are obtained. The confidence intervals, meantime between failures, and B1 are determined. Similarly, for the data on Space Shuttle Main Engine (SSME) blade failures the maximum likelihood estimators are obtained for the Weibull parameters. The variances, confidence intervals, meantime between failures, and reliability are determined. The analysis is performed under assumption of grouped data as well as randomly placed data.

Patil, S. A.↗
NASA NTRS2004

Invariant Theory for Dispersed Transverse Isotropy: An Efficient Means for Modeling Fiber Splay

Most soft tissues possess an oriented architecture of collagen fiber bundles, conferring both anisotropy and nonlinearity to their elastic behavior. Transverse isotropy has often been assumed for a subset of these tissues that have a single macroscopically-identifiable preferred fiber direction. Micro-structural studies, however, suggest that, in some tissues, collagen fibers are approximately normally distributed about a mean preferred fiber direction. Structural constitutive equations that account for this dispersion of fibers have been shown to capture the mechanical complexity of these tissues quite well. Such descriptions, however, are computationally cumbersome for two-dimensional (2D) fiber distributions, let alone for fully three-dimensional (3D) fiber populations. In this paper, we develop a new constitutive law for such tissues, based on a novel invariant theory for dispersed transverse isotropy. The invariant theory is based on a novel closed-form splay invariant that can easily handle 3D fiber populations, and that only requires a single parameter in the 2D case. The model is polyconvex and fits biaxial data for aortic valve tissue as accurately as the standard structural model. Modification of the fiber stress-strain law requires no re-formulation of the constitutive tangent matrix, making the model flexible for different types of soft tissues. Most importantly, the model is computationally expedient in a finite-element analysis.

Freed, alan D.↗
NASA NTRS2020

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗
NASA NTRS1994

The X-ray spectra of blazars: Analysis of the complete EXOSAT archive

We analyzed the 0.1-10 keV spectra of 26 blazars (21 BL Lac objects and 5 highly polarized quasars), on the basis of 93 exposures taken from the EXOSAT archives. Fits were performed first with a single power-law model. Indications are found that better fits can be obtained with models where the slope steepens at higher energies. We therefore considered a broken power law and found that in a large fraction of the spectra the fit is significantly improved. Fits with a power law + oxygen edge at 0.6 keV are also explored.

Sambruna, Rita M.↗
NASA NTRSDate not supplied

Electric Aircraft Propulsion Power System Impedance Modelling Methodology

This paper describes techniques to be used to model megawatt Electrified Aircraft Propulsion (MW EAP) components, specifically motor/generators and inverters, DC-DC converters, and long power leads; resulting in impedance-based models which require no information from the component manufacturers. A custom curve fitting method is described, and an example power supply model is fit using measured impedance data. A method called subscale measurement, employing measurement equipment which has much lower power capability than the MW EAP system under test, is described, simulated, and used to generate an impedance transfer function model of a load. Finally, stability analysis using a Nyquist based approach is performed on the combined supply and load system, using the generated models.

Megawatt electrified aircraft propulsion↗
NASA NTRS1984

Structure of the Uranian rings. I - Square-well model and particle-size constraints

A least squares fitting is added to the diffraction model for the Uranian rings and the model is then used to determine the relative optical depths of the rings at visible and near-IR wavelengths. Light transmission is assumed constant on a scale of a few radial km in the rings. The amount of starlight passing through the rings is calculated from occultations with account taken of a monolayer ring, diffraction effects, and multiple ring layers. A diffracted occultation profile is generated with consideration given to the angular diameter of the star, the impulse response of the detector, and atmospheric aberrations. Calculation procedures which remove noise are outlined and results are compared with the observational data base. The wavelength dependence of optical depth is formulated. It is shown that occultation data do not support the presence of a large fraction of submicron particles. No optical depth variations are projected for 0.88-2.2 microns wavelengths. Unchopped data are needed to verify the model predictions.

Elliot, J. L.↗
NASA NTRS2001

Time-dependent oral absorption models

The plasma concentration-time profiles following oral administration of drugs are often irregular and cannot be interpreted easily with conventional models based on first- or zero-order absorption kinetics and lag time. Six new models were developed using a time-dependent absorption rate coefficient, ka(t), wherein the time dependency was varied to account for the dynamic processes such as changes in fluid absorption or secretion, in absorption surface area, and in motility with time, in the gastrointestinal tract. In the present study, the plasma concentration profiles of propranolol obtained in human subjects following oral dosing were analyzed using the newly derived models based on mass balance and compared with the conventional models. Nonlinear regression analysis indicated that the conventional compartment model including lag time (CLAG model) could not predict the rapid initial increase in plasma concentration after dosing and the predicted Cmax values were much lower than that observed. On the other hand, all models with the time-dependent absorption rate coefficient, ka(t), were superior to the CLAG model in predicting plasma concentration profiles. Based on Akaike's Information Criterion (AIC), the fluid absorption model without lag time (FA model) exhibited the best overall fit to the data. The two-phase model including lag time, TPLAG model was also found to be a good model judging from the values of sum of squares. This model also described the irregular profiles of plasma concentration with time and frequently predicted Cmax values satisfactorily. A comparison of the absorption rate profiles also suggested that the TPLAG model is better at prediction of irregular absorption kinetics than the FA model. In conclusion, the incorporation of a time-dependent absorption rate coefficient ka(t) allows the prediction of nonlinear absorption characteristics in a more reliable manner.

Non-NASA Center↗
NASA NTRS2019

Retrieval of Snow Freeboard of Antarctic Sea Ice Using Waveform Fitting of CryoSat-2 Returns

In this paper we develop a CryoSat-2 algorithm to retrieve the surface elevation of the air–snow interface over Antarctic sea ice. This algorithm utilizes a two-layer physical model that accounts for scattering from a snow layer atop sea ice as well as scattering from below the snow surface. The model produces waveforms that are fit to CryoSat-2 level 1B data through a bounded trust region least-squares fitting process. These fit waveforms are then used to track the air–snow interface and retrieve the surface elevation at each point along the CryoSat-2 ground track, from which the snow freeboard is computed. To validate this algorithm, we compare retrieved surface elevation measurements and snow surface radar return power levels with those from Operation IceBridge, which flew along a contemporaneous CryoSat-2 orbit in October 2011 and November 2012. Average elevation differences (standard deviations) along the flight lines (IceBridge Airborne Topographic Mapper, ATM – CryoSat-2) are found to be 0.016 cm (29.24 cm) in 2011 and 2.58 cm (26.65 cm) in 2012. The spatial distribution of monthly average pan-Antarctic snow freeboard found using this method is similar to what was observed from NASA's Ice, Cloud, and land Elevation Satellite (ICESat), where the difference (standard deviation) between October 2011–2017 CryoSat-2 mean snow freeboard and spring 2003–2007 mean freeboard from ICESat is 1.92 cm (9.23 cm). While our results suggest that this physical model and waveform fitting method can be used to retrieve snow freeboard from CryoSat-2, allowing for the potential to join laser and radar altimetry data records in the Antarctic, larger (∼30 cm) regional differences from ICESat and along-track differences from ATM do exist, suggesting the need for future improvements to the method. Snow–ice interface elevation retrieval is also explored as a potential to obtain snow depth measurements. However, it is found that this retrieval method often tracks a strong scattering layer within the snow layer instead of the actual snow–ice interface, leading to an overestimation of ice freeboard and an underestimation of snow depth in much of the Southern Ocean but with promising results in areas such as the East Antarctic sector.

Fons, Steven W.↗
NASA NTRS1992

Variations in lithospheric thickness on Venus

Recent analyses of Magellan data have indicated many regions exhibiting topograhic flexure. On Venus, flexure is associated predominantly with coronae and the chasmata with Aphrodite Terra. Modeling of these flexural signatures allows the elastic and mechanical thickness of the lithosphere to be estimated. In areas where the lithosphere is flexed beyond its elastic limit the saturation moment provides information on the strength of the lithosphere. Modeling of 12 flexural features on Venus has indicated lithospheric thicknesses comparable with terrestrial values. This has important implications for the venusian heat budget. Flexure of a thin elastic plate due simultaneously to a line load on a continuous plate and a bending moment applied to the end of a broken plate is considered. The mean radius and regional topographic gradient are also included in the model. Features with a large radius of curvature were selected so that a two-dimensional approximation could be used. Comparisons with an axisymmetric model were made for some features to check the validity of the two-dimensional assumption. The best-fit elastic thickness was found for each profile crossing a given flexural feature. In addition, the surface stress and bending moment at the first zero crossing of each profile were also calculated. Flexural amplitudes and elastic thicknesses obtained for 12 features vary significantly. Three examples of the model fitting procedures are discussed.

Johnson, C. L.↗
NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗
NASA NTRS1996

X-Ray Spectral Variability in NGC 7469

We present analyses of two Ginga observations and two observations from the ROSAT database of NGC 7469, focusing on the spectral variability observed on timescales of days and longer. During the 1988 Ginga observation, the hardness ratio (8-21 keV/3.4-5.7 keV) increased significantly as the total flux decreased by 30%. As the spectrum is well fit by the reflection model and since the spectra variability dominates the higher energy band, this could be explained by either a variation in the power law index or in the effective covering fraction of the reflecting material. This ambiguity is inherent in reflection modeling of Ginga spectra from moderate flux Seyfert 1 galaxies. Assuming that the power law index did not change, we find that the reflected flux is consistent with being constant, suggesting that much of the reflecting material may be located more than 3 light-days from the continuum source with the molecular torus being a plausible site. This scenario is also supported by the report of a narrow rather than broad iron K-alpha line in the ASCA data by Guainazzi et al. NGC 7469 was faint during the 1989 Ginga observation, but variability was observed with doubling timescale of 5 hr, and the spectrum was harder. A reflection component could not be constrained, and the change in the spectrum could be explained by an increase in neutral absorption. The brighter of two ROSAT spectra was significantly softer, and in both spectra there was evidence of spectral complexity, as has been previously reported by Turner, George, & Mushotzky and Brandt et al. The spectrum could be fit by a variety of two-component models, including a warm absorber model, an ionized disk model, and a thermal model with single-component blackbody spectrum, but joint fitting of the 1988 average Ginga spectrum and the nonsimultaneous ROSAT spectra favored thermal models, and other models required an anomalously high reflection ratio. This model is supported by the observation of a soft excess component and the lack of ionized absorption edges in the ASCA spectrum by Guainazzi et al. The long-term spectral variability could be explained by relative variability between the power-law and soft excess component normalizations, perhaps implying that hard X-ray reprocessing in thermal material does not dominate on long timescales.

Leighly, Karen↗
NASA NTRS1980

Statistical Analysis of daily mean temperatures

Data of daily mean temperatures recorded at the Kennedy Center during the period of 1957-1977 were analyzed to forecast daily mean temperatures and their thirty-day moving averages for a period of ten to fifteen days in a given month. Since it is found that the standard deviation is linear in the mean, a logarithmic transformation of the data is used for finding an integrated moving average process IMA by the Box-Jenkins aproach. The first differences of the transformed data seem to fit a moving average model with parameter value 2, MA(2). The consideration of seasonality factor makes the fit worse.

Ross, D. C.↗