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,711 records · Page 95

DISCOVERY OF AN ULTRAMASSIVE PULSATING WHITE DWARF

We announce the discovery of the most massive pulsating hydrogen-atmosphere white dwarf (WD) ever discovered, GD 518. Model atmosphere fits to the optical spectrum of this star show it is a 12, 030 {+-} 210 K WD with a log g =9.08 {+-} 0.06, which corresponds to a mass of 1.20 {+-} 0.03 M{sub Sun }. Stellar evolution models indicate that the progenitor of such a high-mass WD endured a stable carbon-burning phase, producing an oxygen-neon-core WD. The discovery of pulsations in GD 518 thus offers the first opportunity to probe the interior of a WD with a possible oxygen-neon core. Such a massive WD should also be significantly crystallized at this temperature. The star exhibits multi-periodic luminosity variations at timescales ranging from roughly 425 to 595 s and amplitudes up to 0.7%, consistent in period and amplitude with the observed variability of typical ZZ Ceti stars, which exhibit non-radial g-mode pulsations driven by a hydrogen partial ionization zone. Successfully unraveling both the total mass and core composition of GD 518 provides a unique opportunity to investigate intermediate-mass stellar evolution, and can possibly place an upper limit to the mass of a carbon-oxygen-core WD, which in turn constrains Type Ia supernovae progenitor systems.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗

The First Multiple Layer Doppler Imaging of an Active Binary

Preliminary results were presented at the Cool Stars, Stellar Systems, and the Sun in Cambridge, MA in July 1997. A copy of the two published papers (in press) is attached. The project has met or exceeded our expectations. The rapid readout data have given us an excellent data set to model in detail the flare behavior. The large number of spectral features observed between the EUVE and HST data have allowed us to compute a mean model atmosphere and compare the results to another well studied system (HR 1099 - Cycle 3 HST observation previously published). I developed a model (anisotropic macroturbulence) that fits the CIV and MgII better than previously achieved. These results have been constrained by the EUVE data. In early studies, 2 gaussians were applied to the profile. The interpretation of these features was unclear. The anisotropic macroturbulence model fits the data better than previously possible and gives a physically reasonable interpretation: there appears to be an asymmetrical distribution between the radial and tangential velocity fields. This is similar to case of the Sun.

Dempsey, Robert C.↗

The 1% concordance Hubble constant

The determination of the Hubble constant has been a central goal in observational astrophysics for nearly a hundred years. Extraordinary progress has occurred in recent years on two fronts: the cosmic distance ladder measurements at low redshift and cosmic microwave background (CMB) measurements at high redshift. The CMB is used to predict the current expansion rate through a best-fit cosmological model. Complementary progress has been made with baryon acoustic oscillation (BAO) measurements at relatively low redshifts. While BAO data do not independently determine a Hubble constant, they are important for constraints on possible solutions and checks on cosmic consistency. A precise determination of the Hubble constant is of great value, but it is more important to compare the high and low redshift measurements to test our cosmological model. Significant tension would suggest either uncertainties not accounted for in the experimental estimates or the discovery of new physics beyond the standard model of cosmology. In this paper we examine in detail the tension between the CMB, BAO, and cosmic distance ladder data sets. We find that these measurements are consistent within reasonable statistical expectations and we combine them to determine a best-fit Hubble constant of 69.6 ± 0.7 km s{sup –1} Mpc{sup –1}. This value is based upon WMAP9+SPT+ACT+6dFGS+BOSS/DR11+H {sub 0}/Riess; we explore alternate data combinations in the text. The combined data constrain the Hubble constant to 1%, with no compelling evidence for new physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Capella: Structure and Abundances

This grant covers the analysis of EUVE spectra of the cool star binary system Capella. This project has also required the analysis of simultaneous Advanced Satellite for Cosmology and Astrophysics (ASCA) data. The ASCA spectrum of Capella could not be fit with standard models; by imposing models based on strong lines observed with EUVE, a problem wavelength region was identified. Correcting the problem required calculations of atomic collision strengths of higher principal quantum number than had ever been calculated. With these new models applied to the ASCA spectrum, better fits were obtained. Findings are that: (1) ASCA and EUVE spectra are both dominated by a region at 6 x 10(exp 6) K. (2) The high energy cut-off of the ASCA spectrum is consistent with emission from the highest ionization stages of EUVE, namely Fe XXIV. (3) EUVE requires a continuous emission measure distribution with more than two temperatures. (4) The ASCA spectra are of such high statistical significance that systematic uncertainties dominate, including atomic physics issues and calibration issues. (5) While the ASCA spectral fits achieve lower Chi(exp 2 with two-temperature fits, the EUVE-derived emission measure distribution models are also consistent with the spectra. (6) The Fe/H ratio obtained from the ASCA fit is within 20 % of the Fe/H abundance obtained from the summed spectra of Capella over 5 EUVE pointings, as well as the 1996 EUVE data. This result confirms our claims that quasi-continua composed of weak emission lines in the short wavelength spectrometer of EUVE are not major contributors to the measured Capella continuum. Other abundance ratios are also determined from the ASCA data, using models derived with EUVE. Si, Si, and Mg appear to be close to solar photospheric values, while the ratio of Ne/Fe is three to four times lower than solar photospheric values. Whether there is a general First Ionization Potential (FIP) effect or a specific neon anomaly cannot be determined from these data. (8) EUVE line ratios indicating low optical depth place limits on the extent to which the X-ray lines are optically, thick. While the brightest lines may be marginally thick, the ASCA spectrum is not of sufficient resolution to determine the optical depth in these lines.

Brickhouse, Nancy S.↗