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 235 records · Page 13

Reduction Kinetics of Perovskite Oxides for Selective Hydrogen Combustion in the Context of Olefin Production

Chemical looping represents a novel approach for generating light olefins in which thermal cracking or catalytic dehydrogenation is coupled with selective hydrogen combustion (SHC) by a metal oxide redox catalyst, which enables autothermal operation, increased per-pass conversion, and greater-than-equilibrium yields. Recent studies indicate that Na 2 WO 4 -promoted perovskite oxides are effective redox catalysts with high olefin selectivity. Herein, kinetic parameters, rates, and reaction models for the reduction of unpromoted and Na 2 WO 4 -promoted CaMnO 3 redox catalysts by H 2 , C 2 H 4 , and C 2 H 6 , is reported. Reduction rates of CaMnO 3 under ethylene and ethane are significantly lower than under H 2 . Additionally, model fitting of reduction kinetics show good agreement with reaction order–controlled models for CaMnO 3 reduction and predict greater oxygen site dependence and higher activation energy for CaMnO 3 reduction by C 2 H 4 as compared with H 2 . Avrami–Erofe'ev nucleation and growth models provide the best fit to the reduction of Na 2 WO 4 /CaMnO 3 in H 2 and in C 2 H 4 . After Na 2 WO 4 promotion, the reduction rate of CaMnO 3 is three orders of magnitude lower in ethylene in comparison to hydrogen, consistent with its superior selectivity to hydrogen combustion. The models developed can be applied toward reactor design and optimization in the context of enhanced olefin production via SHC under a cyclic redox scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models of the Hydrodynamic Histories of Post-AGB Stars. I. Multiflow Shaping of OH 231.8+04.2

We present a detailed hydrodynamic model that matches the present structure of the well-observed preplanetary nebula (“pPN”) OH 231.8+04.2 (“OH231”). The purpose of the model is to present a physically justified and coherent picture of its evolutionary history from about 100 years from the start of the formation of its complex outer structures to the present. We have adopted a set of initial conditions that are heavily constrained by high-quality observations of its present structure and kinematics. The shaping of the nebula occurs while the densities of the flows are “light,” i.e., less than the surrounding AGB-wind environment. The simulations show that pairs of essentially coeval clumps and sprays of the same extent and density, but different outflow speeds, sculpted both the pair of thin axial flow “or spine” and the bulbs. The total ejected mass and momentum in the best-fit model are surprisingly large—3 M {sub ⊙} and 2.2 × 10{sup 41} gm cm s{sup −1}, respectively—however, these values are reduced by up to a factor of 10 in other models that fit the data almost as well. Our ultimate goal is to combine the present model results of masses, momenta, flow speeds, and flow geometries for OH231 with those of other models to be published in the future in order to find common attributes of their ejection histories.

79 ASTRONOMY AND ASTROPHYSICS↗

INTEGRAL Observations of the Enigmatic Be Stars (gamma) Cassiopeiae and HD 110432

We present the results of a hard X-ray study of the Be stars gamma Cassiopeiae and HD 110432 based on observations made with the INTEGRAL observatory. These stars are known to be moderately strong, X-ray sources (L(sub x) approx. equal to = 10(sup 32)-10(sup 33) erg per second). These values are at the extreme high end of the known luminosity distribution for active coronal systems, but several orders of magnitude below typical X-ray binaries. The hard X-ray spectra for these systems are quite similar. They can be well fitted by either optically thin thermal plasma models with kT = 12.5 - 14 keV or a cutoff powerlaw + gaussian line model with photon indices in the 1.3 - 1.5 range and a line energy of 6.7 keV. The 20-50 keV light curves show no evidence for flaring and no significant evidence for periodic variability. It has been proposed that the X-ray emission is due to either accretion onto a white dwarf companion or magnetic activity near the surface of the Be star. We discuss in detail the pros and cons of each scenario towards explaining our spectral and temporal results. Given that both thermal and nonthermal models fit the data equally well, we cannot use the spectra to delineate between these two scenarios. Recent observations indicate that gamma Cas has a approx. 1 solar mass companion in a 203.59 day orbit. This is consistent with the white dwarf - Be star binary model but the lack of periodic modulation of the flux on this timescale calls this conclusion into question. On the other hand the lack of flaring activity may rule against the magnetic activity model. We discuss advances in observations and theory that need to be made to resolve the origin of these systems.

Sturner, S. J.↗

Response of CR-39 nuclear track detectors to protons with non-normal incidence

This paper presents data from experiments with protons at non-normal incidence to CR-39 nuclear track detectors, analyzing the properties of detection efficiency, proton track diameter, track contrast, and track eccentricity. Understanding the CR-39 response to protons incident at an angle is important for designing charged particle detectors for inertial confinement fusion (ICF) applications. This study considers protons with incident energies less than 3 MeV. In this regime, an incident angle of 10° has no effect on CR-39 detection efficiency, and >85% detection efficiency is preserved up through 25° in the range of 1.0 MeV–2.1 MeV. For ICF applications, incident angles above 30° are deemed impractical for detector design due to significant drops in proton detection at all energies. We observe significant reductions in detection efficiency compared to theoretical predictions, particularly at low energies where proton tracks are etched away. The proton track diameter measured by the scan system is observed to decrease with higher incident angles. The track diameters are analyzed with two fitting models, and it is shown that the diameter–energy relation can be fit with the existing models at angles up to 30° . The optical contrast of the tracks tends to increase with the angle, meaning that the tracks are fainter, and a larger increase is observed for higher energies. Eccentricity, a measure of how elongated proton tracks are, increases with the incident angle and drops after the critical angle. The lowest energy tracks remain nearly circular even at higher angles.

47 OTHER INSTRUMENTATION↗

Response of CR-39 nuclear track detectors to protons with non-normal incidence

This paper presents data from experiments with protons at non-normal incidence to CR-39 nuclear track detectors, analyzing the properties of detection efficiency, proton track diameter, track contrast, and track eccentricity. Understanding the CR-39 response to protons incident at an angle is important for designing charged particle detectors for inertial confinement fusion (ICF) applications. This study considers protons with incident energies less than 3 MeV. In this regime, an incident angle of 10° has no effect on CR-39 detection efficiency, and >85% detection efficiency is preserved up through 25° in the range of 1.0 MeV–2.1 MeV. For ICF applications, incident angles above 30° are deemed impractical for detector design due to significant drops in proton detection at all energies. We observe significant reductions in detection efficiency compared to theoretical predictions, particularly at low energies where proton tracks are etched away. The proton track diameter measured by the scan system is observed to decrease with higher incident angles. The track diameters are analyzed with two fitting models, and it is shown that the diameter–energy relation can be fit with the existing models at angles up to 30°. The optical contrast of the tracks tends to increase with the angle, meaning that the tracks are fainter, and a larger increase is observed for higher energies. Eccentricity, a measure of how elongated proton tracks are, increases with the incident angle and drops after the critical angle. The lowest energy tracks remain nearly circular even at higher angles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Missing components in ΛCDM from DESI Y1 baryonic acoustic oscillation measurements: Insights from redshift remapping

We explore transformations of the Friedman-Lemaître-Robertson-Walker (FLRW) metric and cosmological parameters that align with observational data while aiming to gain insights into potential extensions of standard cosmological models. We modified the FLRW metric by introducing a scaling factor, e 2Θ(a) –the cosmological scaling function (CSF), which alters the standard relationship between cosmological redshift and the cosmic scale factor without affecting angular measurements or cosmic microwave background (CMB) anisotropies. Using data from DESI Year 1, Pantheon+ supernovae, and the Planck CMB temperature power spectrum, we constrained both the CSF and cosmological parameters through a Markov chain Monte Carlo approach. Our results indicate that the CSF model fits observational data with a lower Hubble constant (although it is compatible with the value given by Planck 2018 within 1σ) and is predominantly dark matter dominated. Additionally, the CSF model produces temperature and lensing power spectra similar to those predicted by the standard model, though with lower values in the CSF model at large scales. We also checked that when fitting a CSF model without dark energy to the data, we obtain a more negative conformal function. This suggests that the CSF model may offer hints about missing elements and opens up a new avenue for exploring physical interpretations of cosmic acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of sea ice models using an uncertainty-based distance metric for multiple model variables: NEW METRIC FOR SEA ICE MODEL VALIDATION

Here, we implement a variance-based distance metric (D n ) to objectively assess skill of sea ice models when multiple output variables or uncertainties in both model predictions and observations need to be considered. The metric compares observations and model data pairs on common spatial and temporal grids improving upon highly aggregated metrics (e.g., total sea ice extent or volume) by capturing the spatial character of model skill. The D n metric is a gamma-distributed statistic that is more general than the χ 2 statistic commonly used to assess model fit, which requires the assumption that the model is unbiased and can only incorporate observational error in the analysis. The D n statistic does not assume that the model is unbiased, and allows the incorporation of multiple observational data sets for the same variable and simultaneously for different variables, along with different types of variances that can characterize uncertainties in both observations and the model. This approach represents a step to establish a systematic framework for probabilistic validation of sea ice models. The methodology is also useful for model tuning by using the D n metric as a cost function and incorporating model parametric uncertainty as part of a scheme to optimize model functionality. We apply this approach to evaluate different configurations of the standalone Los Alamos sea ice model (CICE) encompassing the parametric uncertainty in the model, and to find new sets of model configurations that produce better agreement than previous configurations between model and observational estimates of sea ice concentration and thickness.

58 GEOSCIENCES↗

High Removal Capacity of Arsenic from Drinking Water Using Modified Magnetic Polyurethane Foam Nanocomposites

In this study, a novel polyurethane foam (PU) nanocomposite adsorbent based on silane-modified magnetic iron-oxide nanoparticles (Fe{sub 3}O{sub 4}@APTES) is synthesized via a low cost and simple in situ polymerization method for the removal of arsenic ions from aqueous solutions. The chemical structure and surface morphology of the prepared nanoparticles and adsorbent were characterized using Fourier transform infrared spectroscopy, attenuated total reflection, scanning electron microscopy, transmission electron microscopy, and energy-dispersive X-ray spectroscopy. Inductively coupled plasma mass spectrometry was used to measure the arsenic concentration of the treated solutions. Sorption isotherms models were applied to determine the adsorption mechanism and modeling parameters. The removal capacity of the modified PU foam was at its highest during a contact time of four hours which resulted in a removal capacity of 95%. Kinetic studies were conducted to determine the adsorption capacity and the uptake rate of arsenic. A Pseudo-order model was found to be the best fit model for adsorption. The prepared adsorbent can be separated from the solution by using an external magnet field.

36 MATERIALS SCIENCE↗

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗

A broadband X-ray study of the Geminga pulsar with NuSTAR And XMM-Newton

We report on the first hard X-ray detection of the Geminga pulsar above 10 keV using a 150 ks observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) observatory. The double-peaked pulse profile of non-thermal emission seen in the soft X-ray band persists at higher energies. Broadband phase-integrated spectra over the 0.2-20 keV band with NuSTAR and archival XMM-Newton data do not fit to a conventional two-component model of a blackbody plus power law, but instead exhibit spectral hardening above ∼5 keV. We find that two spectral models fit the data well: (1) a blackbody (kT {sub 1} ∼ 42 eV) with a broken power law (Γ{sub 1} ∼ 2.0, Γ{sub 2} ∼ 1.4 and E {sub break} ∼ 3.4 keV) and (2) two blackbody components (kT {sub 1} ∼ 44 eV and kT {sub 2} ∼ 195 eV) with a power-law component (Γ ∼ 1.7). In both cases, the extrapolation of the Rayleigh-Jeans tail of the thermal component is consistent with the UV data, while the non-thermal component overpredicts the near-infrared data, requiring a spectral flattening at E ∼ 0.05-0.5 keV. While strong phase variation of the power-law index is present below ∼5 keV, our phase-resolved spectroscopy with NuSTAR indicates that another hard non-thermal component with Γ ∼ 1.3 emerges above ∼5 keV. The spectral hardening in non-thermal X-ray emission as well as spectral flattening between the optical and X-ray bands argue against the conjecture that a single power law may account for multi-wavelength non-thermal spectra of middle-aged pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud

The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- andintermediate-mass stars is one of the most uncertain phases of stellar evolution and themodels need to be calibrated with the aid of observations. To this purpose, we couplehigh-quality observations of resolved stars in the Small Magellanic Cloud (SMC) withdetailed stellar population synthesis simulations computed with theTRILEGALcode.The strength of our approach relies on the detailed spatially-resolved star formationhistory of the SMC, derived from the deep near-infrared photometry of the VISTAsurvey of the Magellanic Clouds, as well as on the capability to quickly and accuratelyexplore a wide variety of parameters and effects with theCOLIBRIcode for the TP-AGB evolution. Adopting a well-characterized set of observations – star counts andluminosity functions – we set up a calibration cycle along which we iteratively changea few key parameters of the TP-AGB models until we eventually reach a good fit tothe observations. Our work leads to identify two best-fitting models that mainly differin the efficiencies of the third dredge-up and mass loss in TP-AGB stars with initialmasses larger than about 3 M. On the basis of these calibrated models we providea full characterization of the TP-AGB stellar population in the SMC in terms of stel-lar parameters (initial masses, C/O ratios, carbon excess, mass-loss rates). Extensivetables of isochrones including these improved models are publicly available

79 ASTRONOMY AND ASTROPHYSICS↗

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Glass Property-Composition Models Update for use in Direct Feed High-Level Waste Flowsheet Development

A set of preliminary glass property models and constraints were developed and augmented by models from literature for use in design of direct-feed high-level waste (DFHLW) glasses for flowsheet evaluation, testing, and design of the Tank Waste Treatment and Immobilization Plant (WTP) high-level waste (HLW) Facility. These models and constraints are meant to be used as a place-holder while glass property-composition data gaps are filled and final plant operating models are developed. This report describes the motivation and intended use of the models, the compilation of data, model fitting and selection, methods to apply the models and constraints in glass design and offers example calculations demonstrating their intended use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Glass Property-Composition Models Update for use in Direct Feed High-Level Waste Flowsheet Development: EWG2.6

A set of preliminary glass property models and constraints were developed and augmented by models from literature for use in design of direct-feed high-level waste (DFHLW) glasses for flowsheet evaluation, testing, and design of the Waste Treatment and Immobilization Plant (WTP) high-level waste (HLW) Facility. These models and constraints are meant to be used as a place-holder while glass property-composition data gaps are filled and final plant operating models are developed. This report describes the motivation and intended use of the models, the compilation of data, model fitting and selection, methods to apply the models and constraints in glass design and offers example calculations demonstrating their intended use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Galilean satellites - Observations of mutual occultations and eclipses in 1973

Seven Galilean satellite mutual events, two occultations and two eclipses of Europa and three eclipses of Io, were observed at three wavelengths (0.35, 0.50, 0.91 micrometers) with a time resolution of 0.1 sec. Preliminary model fits to the light curves are presented. Model satellites with different albedo distributions (uniform disk, bright solar caps, bright quadrant) are used in generating model occultation and eclipse curves to demonstrate the sensitivity of observed light curves to the brightness distribution on the surface of the occulted or eclipsed satellite. At the present the observations yield no conclusive information of the limb-darkening of Io. The best data for Europa indicate that the satellite is limb-darkened at both 0.50 and 0.91 micrometers.

Wasserman, L. H.↗

The phase of the ten-hour modulation in the Jovian magnetosphere /Pioneers 10 and 11/

The paper describes the study of the phase of the 10-hour modulation of energetic electrons seen by Pioneers 10 and 11 in the Jovian magnetosphere. Attention is given to the peaks rather than the valleys of each cycle because the peaks are where physically interesting features occur, such as particle acceleration, current sheets, etc. To identify the peaks, it is required that the instantaneous intensity be higher than the 5-hour running average and the 5-hour running average be greater than the 10-hour running average. These criteria select an interval rather than a point and it is determined that this interval is an appropriate estimate of the experimental uncertainty. When the phases of the peaks are plotted together, they create patterns which are discussed in terms of disk-like, clock-like, and rotating anomaly models of the magnetosphere. Each model fits some of the data, but no model explains all of the data convincingly. It is concluded that there is still no understanding of the configuration of the outer Jovian magnetosphere.

Fillius, W.↗

Interacting binary galaxies. V - NGC 4782/4783 (3C 278): Unbound colliders, not a supermassive pair

CCD imaging data of the elliptical galaxies NGC 4782 and NGC 4783 have been analyzed along with published spectroscopic measurements in order to develop a viable interaction model for this system. The best-fit model suggests that the binary orbit was hyperbolic before the encounter, but that the energy loss during the collision has led to a loosely bound system. The total system mass in the region occupied by the luminous matter is shown to be 1.4 x 10 to the 12th solar masses. The results support the previous theory that the internal velocity dispersions of colliding galaxies tend to increase just past closest approach.

Borne, Kirk D.↗

Chandra spectral measurements of the O supergiant ζ Puppis indicate a surprising increase in the wind mass-loss rate over 18 yr

New long Chandra grating observations of the O supergiant ζ Pup show not only a brightening of the X-ray emission line flux of 13 per cent in the 18 yr since Chandra’s first observing cycle, but also clear evidence – at more than 4σsignificance – of increased wind absorption signatures in its Doppler-broadened line profiles. We demonstrate this with non-parametric analysis of the profiles as well as Gaussian fitting and then use line-profile model fitting to derive a mass-loss rate of 2.47±0.09×10−6Myr−1, which is a 40 per cent increase over the value obtained from the cycle 1 data. The increase in the individual emission line fluxes is greater for short-wavelength lines than long-wavelength lines, as would be expected if a uniform increase in line emission is accompanied by an increase in the wavelength-dependent absorption by the cold wind in which the shock-heated plasma is embedded.

David H. Cohen↗