Search NASA⌕ Search

SEARCH · Search NASA

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,639 records · Page 91

Machine learning-based surrogate modeling for data-driven optimization: a comparison of subset selection for regression techniques

Optimization of simulation-based or data-driven systems is a challenging task, which has attracted significant attention in the recent literature. A very efficient approach for optimizing systems without analytical expressions is through fitting surrogate models. Due to their increased flexibility, nonlinear interpolating functions, such as radial basis functions and Kriging, have been predominantly used as surrogates for data-driven optimization; however, these methods lead to complex nonconvex formulations. Alternatively, commonly used regression-based surrogates lead to simpler formulations, but they are less flexible and inaccurate if the form is not known a priori. In this work, we investigate the efficiency of subset selection regression techniques for developing surrogate functions that balance both accuracy and complexity. Subset selection creates sparse regression models by selecting only a subset of original features, which are linearly combined to generate a diverse set of surrogate models. Five different subset selection techniques are compared with commonly used nonlinear interpolating surrogate functions with respect to optimization solution accuracy, computation time, sampling requirements, and model sparsity. Furthermore, our results indicate that subset selection-based regression functions exhibit promising performance when the dimensionality is low, while interpolation performs better for higher dimensional problems.

42 ENGINEERING↗

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗

High-yield hydrogen production from biomass by in vitro metabolic engineering: Mixed sugars coutilization and kinetic modeling

The use of hydrogen (H 2 ) as a fuel offers enhanced energy conversion efficiency and tremendous potential to decrease greenhouse gas emissions, but producing it in a distributed, carbon-neutral, low-cost manner requires new technologies. Herein we demonstrate the complete conversion of glucose and xylose from plant biomass to H 2 and CO 2 based on an in vitro synthetic enzymatic pathway. Glucose and xylose were simultaneously converted to H 2 with a yield of two H 2 per carbon, the maximum possible yield. Parameters of a nonlinear kinetic model were fitted with experimental data using a genetic algorithm, and a global sensitivity analysis was used to identify the enzymes that have the greatest impact on reaction rate and yield. After optimizing enzyme loadings using this model, volumetric H 2 productivity was increased 3-fold to 32 mmol H 2 ∙L ₋1 ∙h ₋1 . The productivity was further enhanced to 54 mmol H 2 ∙L ₋1 ∙h ₋1 by increasing reaction temperature, substrate, and enzyme concentrations—an increase of 67-fold compared with the initial studies using this method. The production of hydrogen from locally produced biomass is a promising means to achieve global green energy production.

08 HYDROGEN↗

PHASE-RESOLVED X-RAY SPECTRA OF MAGNETARS AND THE CORONAL OUTFLOW MODEL

We test a model recently proposed for the persistent hard X-ray emission from magnetars. In the model, hard X-rays are produced by a decelerating electron-positron flow in the closed magnetosphere. The flow decelerates as it radiates its energy away via resonant scattering of soft X-rays, then it reaches the top of the magnetic loop and annihilates there. We test the model against observations of three magnetars: 4U 0142+61, 1RXS J1708-4009, and 1E 1841-045. We find that the model successfully fits the observed phase-resolved spectra. We derive constraints on the angle between the rotational and magnetic axes of the neutron star, the object inclination to the line of sight, and the size of the active twisted region filled with the plasma flow. Using the fit of the hard X-ray component of the magnetar spectrum, we revisit the remaining soft X-ray component. We find that it can be explained by a modified two-temperature blackbody model. The hotter blackbody is consistent with a hot spot covering 1%-10% of the neutron star surface. Such a hot spot is expected at the base of the magnetospheric e {sup ±} outflow, as some particles created in the e {sup ±} discharge flow back and bombard the stellar surface.

79 ASTRONOMY AND ASTROPHYSICS↗

3D Material Response Analysis of PICA Pyrolysis Experiments

Primarily interested in improving ablation modeling for use in inverse reconstruction of flight environments on ablative heat shields. Ablation model is essentially a component of the heat flux sensor, so model uncertainties lead to measurement uncertainties. Non-equilibrium processes have been known to be significant in low density ablators for a long time, but increased accuracy requirements of the reconstruction process necessitates incorporating this physical effect. Attempting to develop a pyrolysis model for implementation in material response based on the PICA data produced by Bessire and Minton. Pyrolysis gas species molar yields as a function of temperature and heating rate. Several problems encountered while trying to fit Arrhenius models to the data led to further investigation of the experimental setup.

Oliver, Brandon A.↗

Radio-gamma-ray connection and spectral evolution in 4C +49.22 (S4 1150+49): the Fermi, Swift and Planck view

We report the Large Area Telescope on board the Fermi Gamma-ray Space Telescope detected a strong γ-ray flare on 2011 May 15 from a source identified as 4C +49.22, a flat spectrum radio quasar (FSRQ) also known as S4 1150+49. This blazar, characterized by a prominent radio–optical–X-ray jet, was in a low γ-ray activity state during the first years of Fermi observations. Simultaneous observations during the quiescent, outburst and post-flare γ-ray states were obtained by Swift, Planck and optical–IR–radio telescopes (Instituto Nacional de Astrofísica, Óptica y Electrónica, Catalina Sky Survey, Very Long Baseline Array [VLBA], Metsähovi). The flare is observed from microwave to X-ray bands with correlated variability and the Fermi, Swift and Planck data for this FSRQ show some features more typical of BL Lac objects, like the synchrotron peak in the optical band that outshines the thermal blue-bump emission, and the X-ray spectral softening. Multi-epoch VLBA observations show the ejection of a new component close in time with the GeV γ-ray flare. In conclusion, the radio-to-γ-ray spectral energy distribution is modelled and fitted successfully for the outburst and the post-flare epochs using either a single flaring blob with two emission processes (synchrotron self-Compton (SSC), and external-radiation Compton), and a two-zone model with SSC-only mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

A statistical model for interpreting computerized dynamic posturography data

Computerized dynamic posturography (CDP) is widely used for assessment of altered balance control. CDP trials are quantified using the equilibrium score (ES), which ranges from zero to 100, as a decreasing function of peak sway angle. The problem of how best to model and analyze ESs from a controlled study is considered. The ES often exhibits a skewed distribution in repeated trials, which can lead to incorrect inference when applying standard regression or analysis of variance models. Furthermore, CDP trials are terminated when a patient loses balance. In these situations, the ES is not observable, but is assigned the lowest possible score--zero. As a result, the response variable has a mixed discrete-continuous distribution, further compromising inference obtained by standard statistical methods. Here, we develop alternative methodology for analyzing ESs under a stochastic model extending the ES to a continuous latent random variable that always exists, but is unobserved in the event of a fall. Loss of balance occurs conditionally, with probability depending on the realized latent ES. After fitting the model by a form of quasi-maximum-likelihood, one may perform statistical inference to assess the effects of explanatory variables. An example is provided, using data from the NIH/NIA Baltimore Longitudinal Study on Aging.

NASA Discipline Neuroscience↗

The Atacama Cosmology Telescope: A Measurement of the 600 less than l less than 8000 Cosmic Microwave Background Power Spectrum at 148 GHz

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 deg(sup 2) of the southern sky, in a 4 deg. 2-wide strip centered on declination 53 deg. South. The CMB at arc minute angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect from galaxy dusters, and to emission by radio sources and dusty galaxies. After masking the 108 brightest point sources in our maps, we estimate the power spectrum between 600 less than l less than 8000 using the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. Our absolute calibration is based on observations of Uranus. To verify the calibration and test the fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and recover the WMAP power spectrum from 250 less than l less than 1150. The power beyond the Silk damping tail of the CMB (l approximately 5000) is consistent with models of the emission from point sources. We quantify the contribution of SZ clusters to the power spectrum by fitting to a model normalized to sigma 8 = 0.8. We constrain the model's amplitude A(sub sz) less than 1.63 (95% CL). If interpreted as a measurement of as, this implies sigma (sup SZ) (sub 8) less than 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to a 6-parameter ACDM model plus point sources and the SZ effect is consistent with these results.

Fowler, J. W.↗

Permeation of CO 2 and N 2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N 2 and CO 2 through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic model development for single step ethanol to butene rich olefin process over Cu-Y/Beta catalysts

Here, this study presents the first intrinsic kinetic model for the single-step conversion of ethanol to butene-rich olefins over bifunctional Cu-Y/Beta catalysts, addressing a critical gap in the design and scale-up of Sustainable Aviation Fuel (SAF) processes. The reaction network comprises ten global steps involving dehydrogenation, aldol condensation, hydrogenation, and dehydration reactions, distributed across Cu and rare-earth (Y) active sites. The model incorporates dual-site functionality (Cu and Y site) and explicitly accounts for key intermediates such as crotonaldehyde and butanal. Reaction rates are formulated using Langmuir–Hinshelwood–Hougen–Watson (LHHW) kinetics. Kinetic parameters are extracted by fitting the model to lab-scale packed-bed reactor data across a wide range of temperatures and space velocities, demonstrating strong agreement in ethanol conversion and product selectivity. The reaction kinetics developed in this work provide a foundational basis for constructing reactor models that enable process optimization and scale-up of ethanol-to-jet fuel technologies.

Cu-Y/Beta catalyst↗

Photometric and spectroscopic studies of massive binaries in the large Magellanic Cloud. II. Three O-type systems in the 30 Dor region

This is the second paper in a series devoted to the study of massive binary systems in the Large Magellanic Cloud (LMC). We mainly aim to provide accurate data that constrains the mass-luminosity relation for the most massive stars but also to address the long lasting problem known as the 'mass discrepancy'. We present here our results for three binaries (LMC 169782, LMC 171520, and [P93] 921) harboring the earliest O-type components—ranging from O4 V to O6.5 V—among our sample of 17 systems. Our photometry provided accurate periods for the studied systems, allowing the spectroscopic observations to be performed at selected phases where the radial velocity separation between binary components is larger. Following the procedure outlined in our first paper of this series, after solving the radial velocity curves for orbital parameters, we used tomographic reconstruction to obtain the individual spectra of each star, from which we determined effective temperatures via a model atmosphere fitting with FASTWIND. This information, combined with the light-curve analysis that was performed with GENSYN, enabled the determination of absolute masses, radii, and bolometric luminosities that are compared with those predicted by modern stellar evolutionary models finding that they agree within the uncertainties. Nevertheless, the comparison seems to confirm the small differences found in the first paper of this series in the sense that the evolutionary masses are slightly larger than the Keplerian ones, with differences averaging ∼10%, or alternatively, the stellar evolutionary models predict luminosities that are somewhat lower than observed. Still, the overall agreement between the current evolutionary models and the empirically determined stellar parameters is remarkable.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectroscopic observations of the detached binary PG 1413 + 015

We present improved estimates of the stellar parameters of the eclipsing, precataclysmic binary system PG 1413 + 015 (GH Vir), which has an orbital period of only 8h16m. Model atmosphere fits a Balmer line profiles yield T(eff) = 48,800 +/- 1200 K and log g = 7.70 +/- 0.11 for the DAO white dwarf primary star, from which a mass of 0.51 +/- 0.04 solar mass is inferred using evolutionary models. An ultraviolet spectrum obtained with the IUE Observatory has a slope consistent with this temperature and the assumption of no interstellar extinction. A red CCD spectrum of the secondary star during the 12-minute total eclipse indicates a spectral type of M3 V-M5 V. Reanalysis of the eclipse light curve leads to an inferred radius of 0.15 solar radius and a mass of 0.10 solar mass for the secondary, the latter being marginally consistent with the spectral type. Reprocessing on the facing side of the secondary produces phase-dependent Balmer line emission and detectable variations in the continuum from 6500-9000 A. The observed levels of reprocessing are consistent with expectations based on the above stellar parameters.

Fulbright, Michael S.↗

The Constitutive Modeling of Thin Films with Randon Material Wrinkles

Material wrinkles drastically alter the structural constitutive properties of thin films. Normally linear elastic materials, when wrinkled, become highly nonlinear and initially inelastic. Stiffness' reduced by 99% and negative Poisson's ratios are typically observed. This paper presents an effective continuum constitutive model for the elastic effects of material wrinkles in thin films. The model considers general two-dimensional stress and strain states (simultaneous bi-axial and shear stress/strain) and neglects out of plane bending. The constitutive model is derived from a traditional mechanics analysis of an idealized physical model of random material wrinkles. Model parameters are the directly measurable wrinkle characteristics of amplitude and wavelength. For these reasons, the equations are mechanistic and deterministic. The model is compared with bi-axial tensile test data for wrinkled Kaptong(Registered Trademark) HN and is shown to deterministically predict strain as a function of stress with an average RMS error of 22%. On average, fitting the model to test data yields an RMS error of 1.2%

Murphey, Thomas W.↗

Status Report on Photothermal Radiometry Development

Thermal transport properties such as thermal conductivity and diffusivity, are critical parameters in the development of improved nuclear fuels and components. It is well known that thermal transport properties change temporally and spatially due to dynamic microstructural evolution as a result of neutron irradiation. Laser based photothermal techniques are an ideal candidate for thermal properties measurement of irradiated materials because they are non contact, non destructive and can be performed remotely. In this report, we present the progress in the development of a photothermal radiometry (PTR) system to measure the thermal diffusivity of solid samples. The PTR technique measures the local temperature response from a modulated heat source by detecting the sample radiation heat flux. Initially, a well developed thermal wave analytical model was fit to the data to obtain the thermal diffusivity. However, systematic errors from the diffraction effect and non linearity issues resulted in thermal diffusivity measurements higher than literature values. These effects have been studied, an improved analytical model has been developed, and the physics behind the phenomena are still under investigation. A procedure for optimizing measurements from the benchtop system has been provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Study of the Influence of Xanthate Derivative Structures on Copper Sulfide Mineral Adsorption Under Acidic Conditions

Adsorption of commercial xanthate derivatives on copper sulfide mineral (covellite, CuS) was studied by kinetics and isotherm adsorption experiments. The adsorption of xanthate derivatives was confirmed by FTIR (Fourier transform infrared spectroscopy) and XPS (X-ray photoelectron spectroscopy) results. Experiments were performed with two different xanthate derivatives, C-4410 (O-pentyl S-2-propenyl ester) and C-4940 (isobutyl xanthogen ethyl formate), on individual doses of 0.05 g of powdered covellite. It was found that the equilibrium times at pH 2, 4, and 6 were different for both xanthate derivatives. The shortest times were achieved at pH 2 and 4. The results suggest that C-4110 can be used as collector in a wide range of pH, while C-4940 is limited to lower pH values. Pseudo first- and pseudo second-order kinetics models were thus applied to the experimental data for pH 2. The information obtained from the kinetics models combined with XPS allowed proposing the adsorption mechanism for the covellite-xanthate derivative pair. The adsorption takes place through a non-covalent interaction for C-4410 and chemisorption process for C-4940. The best-fitting isotherm models for C-4410 and C-4940 adsorption were Redlich–Peterson and Freundlich, respectively, which yield a maximum adsorption capacity of 57.07 mg g{sup −1} for C-4410 and 44.62 mg g{sup −1} for C-4940.

36 MATERIALS SCIENCE↗

Quantifying the effect of multiple load flexibility strategies on commercial building electricity demand and services via surrogate modeling

The expansion of commercial building demand response as a demand-side management resource for the electric grid necessitates new decision support resources for customers seeking to assess the benefit–risk tradeoffs of possible strategies for energy flexible building operations. To address this need, we, in this study, develop surrogate models that predict the impacts of several load flexibility strategies on commercial building electricity demand and indoor temperature, focusing on offices and retail buildings at multiple scales. The surrogate models are fit to a synthetic database generated via whole building simulations, which establish the relationships between the key operational features of a given strategy and potential changes in building demand and temperature across a variety of contexts. The surrogate models are translated to a Bayesian framework to allow straightforward communication of uncertainty and parameter updating given new evidence. We find strong predictive performance across the suite of models, underscoring the usefulness of the approach in guiding decisions about implementing load flexibility strategies under a particular set of operational and environmental conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

PET Imaging of Fatty Acid Amide Hydrolase with [ 18 F]DOPP in Nonhuman Primates

Here, the fatty acid amide hydrolase (FAAH) regulates endocannabinoid signaling. [ 11 C]CURB, an irreversibly binding FAAH inhibitor, has been developed for clinical research imaging with PET. However, no fluorine-18 labeled radiotracer for FAAH has yet advanced to human studies. [ 18 F]DOPP ([ 18 F]3-(4,5-dihydrooxazol-2-yl)phenyl (5-fluoropentyl)carbamate) has been identified as a promising 18 F-labeled analogue based on rodent studies. The goal of this work is to evaluate [ 18 F]DOPP in nonhuman primates to support its clinical translation. High specific activity [ 18 F]DOPP (5–6 Ci·μmol –1 ) was administered intravenously (iv) to three baboons (2M/1F, 3–4 years old). The distribution and pharmacokinetics were quantified following a 2 h dynamic imaging session using a simultaneous PET/MR scanner. Pretreatment with the FAAH-selective inhibitor, URB597, was carried out at 200 or 300 μg/kg iv, 10 min prior to [ 18 F]DOPP administration. Rapid arterial blood sampling for the first 3 min was followed by interval sampling with metabolite analysis to provide a parent radiotracer plasma input function that indicated ~95% baseline metabolism at 60 min and a reduced rate of metabolism after pretreatment with URB597. Regional distribution data were analyzed with 1-, 2-, and 3-tissue compartment models (TCMs), with and without irreversible trapping since [ 18 F]DOPP covalently links to the active site of FAAH. Consistent with previous findings for [ 11 C]CURB, the 2TCM with irreversible binding was found to provide the best fit for modeling the data in all regions. The composite parameter λk 3 was therefore used to evaluate whole brain (WB) and regional binding of [ 18 F]DOPP. Pretreatment studies showed inhibition of λk 3 across all brain regions (WB baseline: 0.112 mL/cm3/min; 300 μg/kg URB597: 0.058 mL/cm 3 /min), suggesting that [ 18 F]DOPP binding is specific for FAAH, consistent with previous rodent data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Continuum model of tensile fracture of metal melts and its application to a problem of high-current electron irradiation of metals

A continuum model of the metal melt fracture is formulated on the basis of the continuum mechanics and theory of metastable liquid. A character of temperature and strain rate dependences of the tensile strength that is predicted by the continuum model is verified, and parameters of the model are fitted with the use of the results of the molecular dynamics simulations for ultra-high strain rates (≥1–10/ns). A comparison with experimental data from literature is also presented for Al and Ni melts. Using the continuum model, the dynamic tensile strength of initially uniform melts of Al, Cu, Ni, Fe, Ti, and Pb within a wide range of strain rates (from 1–10/ms to 100/ns) and temperatures (from melting temperature up to 70–80% of critical temperature) is calculated. The model is applied to numerical investigation of a problem of the high-current electron irradiation of Al, Cu, and Fe targets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗