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 379 records · Page 21

Evaluation of an Image-Derived Input Function for Kinetic Modeling of Nicotinic Acetylcholine Receptor-Binding PET Ligands in Mice

Positron emission tomography (PET) radioligands that bind with high-affinity to α4β2-type nicotinic receptors (α4β2Rs) allow for in vivo investigations of the mechanisms underlying nicotine addiction and smoking cessation. Here, we investigate the use of an image-derived arterial input function and the cerebellum for kinetic analysis of radioligand binding in mice. Two radioligands were explored: 2-[ 18 F]FA85380 (2-FA), displaying similar pKa and binding affinity to the smoking cessation drug varenicline (Chantix), and [ 18 F]Nifene, displaying similar pKa and binding affinity to nicotine. Time–activity curves of the left ventricle of the heart displayed similar distribution across wild type mice, mice lacking the β2-subunit for ligand binding, and acute nicotine-treated mice, whereas reference tissue binding displayed high variation between groups. Binding potential estimated from a two-tissue compartment model fit of the data with the image-derived input function were higher than estimates from reference tissue-based estimations. Rate constants of radioligand dissociation were very slow for 2-FA and very fast for Nifene. We conclude that using an image-derived input function for kinetic modeling of nicotinic PET ligands provides suitable results compared to reference tissue-based methods and that the chemical properties of 2-FA and Nifene are suitable to study receptor response to nicotine addiction and smoking cessation therapies.

2-FA85380↗

Silicon vidicon imaging of Jupiter - 4100- to 8300-A absolute reflectivities and limb darkening of spatially resolved regions

Jupiter was observed in six continuum wavelength channels in the region 4100-8300 A, using a silicon vidicon imaging photometer. Spectral reflectivities and high spatial resolution limb-darkening curves for several belts and zones have been extracted from the data. Simple model fits to the data yield information regarding spectral and spatial variations in single-scattering albedos and shape of particle single-scattering phase functions. Belts appear to be more backscattering than zones, particularly in the blue. The data are in moderate agreement with limb-darkening predicted by models derived from the center-to-limb variation in equivalent width of the H2 4-0 S(1) quadrupole line (Cochran, 1976) in the South Tropical Zone, but strongly disagree with the results of such models for the North Equatorial Belt.

Diner, D. J.↗

Mapping tall shrub biomass in Alaska at landscape scale using structure-from-motion photogrammetry and lidar

Warming in arctic and boreal regions is increasing shrub cover and biomass. In southcentral Alaska, willow (Salix spp.) and alder (Alnus spp.) shrubs grow taller than many tree species and account for a substantial proportion of aboveground biomass, yet they are not individually measured as part of the operational Forest Inventory and Analysis (FIA) Program. The goal of this research was to test methods for landscape-scale mapping of tall shrub biomass in upper montane and subalpine environments using FIA-type plot measurements (n = 51) and predictor variables from imagery-based structure-from-motion (SfM) and airborne lidar. Specifically, we compared biomass models constructed from imagery acquired by unmanned aerial vehicle (UAV; ~1.7 cm pixels), imagery from the NASA Goddard's Lidar, Hyperspectral, and Thermal Airborne Imager (G-LiHT; ~3.1 cm pixels), and concomitant G-LiHT small-footprint lidar. Tall shrub biomass was most accurately predicted at 5 m resolution (R^2 = 0.81, RMSE = 1.09 kg/sq. m) using G-LiHT SfM color and structure variables. Lidar-only models had lower precision (R^2 = 0.74, RMSE = 1.26 kg/sq. m), possibly due to reduced model information content from variable multicollinearity or lower data density. Separate models for upper montane zones with trees and shrubs and subalpine zones with only shrubs were always chosen over single models based on minimization of Akaike's Information Criterion, indicating the need for variable sets robust to overhanging tree canopy. Decreasing point density from UAV (5000–8000 pts./sq. m) to the G-LiHT SfM point cloud (500–2000 pts./sq. m) had little impact on model fit, suggesting that high-resolution airborne imagery can extend SfM approaches well beyond line-of-sight restrictions for UAV platforms. Overall, our results confirmed that SfM from high-resolution imagery is a viable approach to estimate shrub biomass in the boreal region, especially when an existing lidar terrain model and local field calibration data are available to quantify uncertainty in the SfM point cloud and landscape-scale estimates of shrub biomass.

Michael Alonzo↗

Unified and Consistent Structure Growth Measurements from Joint ACT, SPT, and Planck CMB Lensing

We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT), and Planck. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing band powers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of A lens recon = 1.025 ± 0.017 with respect to the theory prediction from the best-fit CMB Planck-ACT cosmology. The band powers from all three lensing datasets, analyzed jointly, yield a 1.6% measurement of the parameter combination S 8 CMBL ≡ σ 8 ( Ω m / 0.3 ) 0.25 = 0.82 5 - 0.013 + 0.015 . Including dark energy spectroscopic instrument baryon acoustic oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to σ 8 = 0.829 ± 0.009 (a 1.1% determination). When combining with uncalibrated supernovae from Pantheon+, we present a 4% sound-horizon-independent estimate of H 0 = 66.4 ± 2.5 km s - 1 Mpc - 1 . The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a Λ CDM model fit to the primary CMB data from Planck and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the Λ CDM cosmological model, the combination of primary CMB, BAO, and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

A Simplified Approach Based on Cellular Automata for Describing Direct Reduced Iron Production in Different Reducing Conditions

A quick computation approach based on cellular automata is developed and implemented to describe the reduction of iron ore pellets by a mixture of reducing agents featured by different H 2 /CO ratios. The evolution of oxygen concentration inside the pellet is followed from the beginning to the end of contact between reducing agent and pellet. The variation of thermal state of pellets and gas mixture is computed based on their initial temperature, considering the heat involved and the convective heat exchange between pellet and gas mixture. The use of cellular automata and finite‐difference method to solve the diffusion equation point out the absence of any diffusion coefficient value, allowing to make the model fit the experimental trial, because the problem is that it is not ruled just by diffusion but also by the concentration variation of reducing agent inside the pellet due to porosity increasing during reduction. The updating of the reducing agents concentration implies a sharp decrease of oxygen concentration that the cellular automata model considers. The developed model is able to provide the in‐line control of reduction process and could be used to adjust the chemical concentration and temperature of injected reducing agents.

Metallurgy & Metallurgical Engineering↗

Measurement of 3 Λ H production in Pb–Pb collisions at √ s NN = 5.02 TeV

The first measurement of 3 Λ H and $^{3}_{\overline{Λ}}$$\overline{H}$ differential production with respect to transverse momentum and centrality in Pb–Pb collisions at √ s NN = 5.02 TeV is presented. The 3 Λ H has been reconstructed via its two-charged-body decay channel, i.e., 3 Λ H → 3 He → π - . A Blast-Wave model fit of the p T -differential spectra of all nuclear species measured by the ALICE collaboration suggests that the kinetic freeze-out surface is consistent with that of other nuclei. The ratio between the integrated yields of 3 Λ H and 3 He is compared to predictions from the statistical hadronisation model and the coalescence model, with the latter being favoured by the presented measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Subsolar Water Abundance

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 μm, in addition to a Spitzer 4.5 μm secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free-retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally ≲0.1× solar (0.001–0.7× solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical-equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they, too, prefer low metallicities ([M/H] ≲ -2, consistent with the free-retrieval results). However, all of the retrievals should be interpreted with some caution because they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find that the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b’s dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b’s composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Archetype-based Redshift Estimation for the Dark Energy Spectroscopic Instrument Survey

We present a computationally efficient galaxy archetype-based redshift estimation and spectral classification method for the Dark Energy Survey Instrument (DESI) survey. The DESI survey currently relies on a redshift fitter and spectral classifier using a linear combination of principal component analysis–derived templates, which is very efficient in processing large volumes of DESI spectra within a short time frame. However, this method occasionally yields unphysical model fits for galaxies and fails to adequately absorb calibration errors that may still be occasionally visible in the reduced spectra. Our proposed approach improves upon this existing method by refitting the spectra with carefully generated physical galaxy archetypes combined with additional terms designed to absorb data reduction defects and provide more physical models to the DESI spectra. We test our method on an extensive data set derived from the survey validation (SV) and Year 1 (Y1) data of DESI. Our findings indicate that the new method delivers marginally better redshift success for SV tiles while reducing catastrophic redshift failure by 10%–30%. At the same time, results from millions of targets from the main survey show that our model has relatively higher redshift success and purity rates (0.5%–0.8% higher) for galaxy targets while having similar success for QSOs. These improvements also demonstrate that the main DESI redshift pipeline is generally robust. Additionally, it reduces the false-positive redshift estimation by 5%–40% for sky fibers. We also discuss the generic nature of our method and how it can be extended to other large spectroscopic surveys, along with possible future improvements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A moving-window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving-window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm, tests for parameter validity, and the on-line design process.

Montgomery, R. C.↗

A moving window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch-processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm made during the development of the system.

Dunn, H. J.↗

Evaluation of a canopy reflectance model for LAI estimation through its inversion

A technique for estimating leaf area indices (LAIs) based on bidirectional canopy reflectance (CR) data is applied to three plant canopies: a naturally growing healthy soybean canopy; and a clumped and tufted orchardgrass canopy, respectively. The CR data were collected using a PARABOLA instrument which is capable of acquiring complete sky-and-ground looking hemispheres in 11 seconds. The model fit and LAI estimates were good for the soybean and clumped orchardgrass canopies, but poor for the tufted orchardgrass canopy when the maximum zenith angle was less than 50 percent. It is shown that the biophysical parameter estimation based on CR measurements applied well to homogeneous herbaceous vegetation types, while better CR models are needed to adequately represent discontinuous plant canopies.

Goel, N. S.↗

An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Subsolar Water Abundance

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 μm, in addition to a Spitzer 4.5 μm secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free-retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally ≲ 0.1× solar (0.001–0.7× solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical-equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they, too, prefer low metallicities ([M/H] ≲ −2, consistent with the free-retrieval results). However, all of the retrievals should be interpreted with some caution because they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find that the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

Knicole D. Colón↗

Resolving topographic detail on Venus by modeling complex Magellan altimetry echoes

Magellan's altimeter is providing some of the finest resolution topography of Venus achieved to date. Nevertheless, efforts continue to improve the topographic resolution whenever possible. One effort to this end is stereoscopic imaging, which provides topography at scales similar to that of the synthetic aperture radar (SAR). However, this technique requires two SAR images of the same site to be obtained and limits the utility of this method. In this paper, we present another method to resolve topographic features at scales smaller than that of an altimeter footprint, which is more globally applicable than the stereoscopic approach. Each pulse which is transmitted by Magellan's altimeter scatters from the planet and echoes to the receiver, delayed based on the distance between the spacecraft and each surface element. As resolved in time, each element of an altimetry echo represents the sum of all points on the surface which are equidistant from the spacecraft. Thus, individual returns, as a function of time, create an echo profile which may be used to derive properties of the surface, such as the scattering law or, in this case, the topography within the footprint. The Magellan project has derived some of this information by fitting model templates to radar echo profiles. The templates are calculated based on Hagfor's Law, which assumes a smooth, gently undulating surface. In most regions these templates provide a reasonable fit to the observed echo profile; however, in some cases the surface departs from these simple assumptions and more complex profiles are observed. Specifically, we note that sub-footprint topographic relief apparently has a strong effect on the shape of the echo profile. To demonstrate the effects of sub-resolution relief on echo profiles, we have calculated the echo shapes from a wide range of simple topographic models. At this point, our topographic models have emphasized surfaces where only two dominant elevations are contained within a footprint, such as graben, ridges, crater rims, and central features in impact craters.

Lovell, Amy J.↗

Characterization of the indoor near-field aerosol transmission in a model commercial office building

To evaluate the exposure potential of infectious aerosols containing SARS-CoV-2 in an office building setting, synthetic test aerosols were used to experimental study airborne particle transmission in a multizone small office test building at the Oak Ridge National Laboratory. Nine measurement points in a single zone using active aerosol impactors report that the coefficient of variation of the time-averaged concentration is <10% in two campaigns and < 15% in one campaign, so a nearly well-mixed condition was noted. To understand the effect of HVAC system operation on the dynamic concentration of aerosols in office spaces, an aerosol transport model that includes factors such as outside air (OA) ratio, filtration, return air fraction, transport loss in air ducts, and particle deposition has been developed. The results of model fitting demonstrate strong agreement with experimental data. Overall, our investigation finds the return air fraction effects outweigh other mechanisms for the aerosol recirculation in this study, and the impact of air change rate (ACR) is more important than the small particle deposition for aerosol removal. Because ACR dominates the aerosol transport, the full model can be simplified to just one factor, the ACR, while maintaining an acceptable representation of the experimental data.

60 APPLIED LIFE SCIENCES↗

Unified and consistent structure growth measurements from joint ACT, SPT and \textit{Planck} CMB lensing

We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) and \textit{Planck}. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing bandpowers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of $A_\mathrm{lens}^\mathrm{recon} = 1.025 \pm 0.017$ with respect to the theory prediction from the best-fit CMB \textit{Planck}-ACT cosmology. The bandpowers from all three lensing datasets, analyzed jointly, yield a $1.6\%$ measurement of the parameter combination $S_8^\mathrm{CMBL} \equiv σ_8\,(Ω_m/0.3)^{0.25} = 0.825^{+0.015}_{-0.013}$. Including Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to $σ_8 = 0.829 \pm 0.009$ (a $1.1\%$ determination). When combining with uncalibrated supernovae from \texttt{Pantheon+}, we present a $4\%$ sound-horizon-independent estimate of $H_0=66.4\pm2.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}} $. The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a $Λ$CDM model fit to the primary CMB data from \textit{Planck} and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the $Λ\mathrm{CDM}$ cosmological model, the combination of primary CMB, BAO and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗

The Suzaku Observation of NGC 3516: Complex Absorption and the Broad and Narrow Fe K Lines

We present results from a 150 ksec Suzaku observation of the Seyfert 1 NGC 3516 in October 2005. The source was in a relatively highly absorbed state. Our best-fit model is consistent with partial covering by a lowly-ionized absorber with a column density near 5x10(exp 22) cm(exp -2) and with a covering fraction 96-100 percent. Narrow K-shell absorption features due to He- and H-like Fe confirm the presence of a high-ionization absorbing component as well. A broad Fe K(alpha) diskline is required in all fits, even after the complex absorption is taken into account; an additional partial-covering component is an inadequate substitute for the continuum curvature associated with the broad line. The narrow Fe Ka line at 6.4 keV is resolved, yielding a velocity width commensurate with the optical Broad Line Region. The strength of the Compton reflection hump suggests a contribution mainly from the broad Fe line origin. We include in our model soft band emission lines from He- and H-like ions and radiative recombination lines, consistent with photo-ionization, though a small contribution from collisional ionization is possible.

Markowitz, Alex↗

The role of geography in the complex diffusion of innovations

The urban–rural divide is increasing in modern societies calling for geographical extensions of social influence modelling. Improved understanding of innovation diffusion across locations and through social connections can provide us with new insights into the spread of information, technological progress and economic development. In this work, we analyze the spatial adoption dynamics of iWiW, an Online Social Network (OSN) in Hungary and uncover empirical features about the spatial adoption in social networks. During its entire life cycle from 2002 to 2012, iWiW reached up to 300 million friendship ties of 3 million users. We find that the number of adopters as a function of town population follows a scaling law that reveals a strongly concentrated early adoption in large towns and a less concentrated late adoption. We also discover a strengthening distance decay of spread over the life-cycle indicating high fraction of distant diffusion in early stages but the dominance of local diffusion in late stages. The spreading process is modelled within the Bass diffusion framework that enables us to compare the differential equation version with an agent-based version of the model run on the empirical network. Although both model versions can capture the macro trend of adoption, they have limited capacity to describe the observed trends of urban scaling and distance decay. We find, however that incorporating adoption thresholds, defined by the fraction of social connections that adopt a technology before the individual adopts, improves the network model fit to the urban scaling of early adopters. Controlling for the threshold distribution enables us to eliminate the bias induced by local network structure on predicting local adoption peaks. Finally, we show that geographical features such as distance from the innovation origin and town size influence prediction of adoption peak at local scales in all model specifications.

97 MATHEMATICS AND COMPUTING↗

Modeling suggests that virion production cycles within individual cells is key to understanding acute hepatitis B virus infection kinetics

Hepatitis B virus (HBV) infection kinetics in immunodeficient mice reconstituted with humanized livers from inoculation to steady state is highly dynamic despite the absence of an adaptive immune response. To recapitulate the multiphasic viral kinetic patterns, we developed an agent-based model that includes intracellular virion production cycles reflecting the cyclic nature of each individual virus lifecycle. The model fits the data well predicting an increase in production cycles initially starting with a long production cycle of 1 virion per 20 hours that gradually reaches 1 virion per hour after approximately 3–4 days before virion production increases dramatically to reach to a steady state rate of 4 virions per hour per cell. Together, modeling suggests that it is the cyclic nature of the virus lifecycle combined with an initial slow but increasing rate of HBV production from each cell that plays a role in generating the observed multiphasic HBV kinetic patterns in humanized mice.

59 BASIC BIOLOGICAL SCIENCES↗