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,261 records · Page 70

LICK INDICES AND SPECTRAL ENERGY DISTRIBUTION ANALYSIS BASED ON AN M31 STAR CLUSTER SAMPLE: COMPARISONS OF METHODS AND MODELS

Application of fitting techniques to obtain physical parameters—such as ages, metallicities, and α-element to iron ratios—of stellar populations is an important approach to understanding the nature of both galaxies and globular clusters (GCs). In fact, fitting methods based on different underlying models may yield different results and with varying precision. In this paper, we have selected 22 confirmed M31 GCs for which we do not have access to previously known spectroscopic metallicities. Most are located at approximately one degree (in projection) from the galactic center. We performed spectroscopic observations with the 6.5 m MMT telescope, equipped with its Red Channel Spectrograph. Lick/IDS absorption-line indices, radial velocities, ages, and metallicities were derived based on the EZ-Ages stellar population parameter calculator. We also applied full spectral fitting with the ULySS code to constrain the parameters of our sample star clusters. In addition, we performed χ{sub min}{sup 2} fitting of the clusters’ Lick/IDS indices with different models, including the Bruzual and Charlot models (adopting Chabrier or Salpeter stellar initial mass functions and 1994 or 2000 Padova stellar evolutionary tracks), the galev, and the Thomas et al. models. For comparison, we collected their UVBRIJK photometry from the Revised Bologna Catalogue (v.5) to obtain and fit the GCs’ spectral energy distributions (SEDs). Finally, we performed fits using a combination of Lick/IDS indices and SEDs. The latter results are more reliable and the associated error bars become significantly smaller than those resulting from either our Lick/IDS indices-only or our SED-only fits.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral Photometric Properties of the Moon

We modeled the solar phase curves of the moon at a series of wavelengths using the full disk telescopic observations [1]. We endeavored to keep the database self-contained, that is, to use the values derived for the solar magnitude and phase curves of the disk-integrated [1]. These observations were made in a suite of 10 narrowband filters between 0.315 microns and 1.06 microns, and in the broad band Johnson UBV filters, as part of a larger program to obtain photoelectric photometry of the larger planets. Two aspects of the lunar observations are unique. First, the observations cover phase angles from 6deg through 120deg. More importantly, the observers used a special 20-mm diameter f/15 fused quartz lens constructed solely for this purpose. The lens reduced the whole lunar image in the focal plane to a size comparable to the planets observed as part of the same program. This image was fed directly into the photometer. Thus, these observations constitute the only existing set of phase curves of the entire lunar disk over a range of wavelengths. Table 1 lists the values of the Hapke model parameters which fit the data. Figure 1 is an example of the model fits to the data.

Dominque, D.↗

Deeper snow increases the net soil organic carbon accrual rate in moist acidic tussock tundra: 210Pb evidence from Arctic Alaska

The net change in the carbon inventory of arctic tundra remains uncertain as global warming leads to shifts in arctic water and carbon cycles. To better understand the response of arctic tundra carbon to changes in winter precipitation amount, we investigated soil depth profiles of carbon concentration and radionuclide activities ( 7 Be, 137 Cs, 210 Pb, and 241 Am) in the active layer of a twenty-two-year winter snow depth manipulation experiment in moist acidic tussock tundra at Toolik Lake, Alaska. Depth correlations of cumulative carbon dry mass (g cm -2 ) vs. unsupported 210 Pb activity (mBq g -1 ) were examined using a modified constant rate of supply (CRS) model. Results were best fit by two-slope CRS models indicating an apparent step temporal increase in the accumulation rate of soil organic carbon. Most of the best-fit model chronologies indicated that the increase in carbon accumulation rate apparently began and persisted after snow fence construction in 1994. The inhomogeneous nature of permafrost soils and their relatively low net carbon accumulation rates make it challenging to establish robust chronologic records. Nonetheless, the data obtained in this study support a decadal-scale increase in net soil organic carbon accumulation rate in the active layer of arctic moist acidic tussock tundra under conditions of increased winter precipitation.

54 ENVIRONMENTAL SCIENCES↗

Pore growth in U-Mo/Al dispersion fuel

U-Mo/Al dispersion fuel is currently under development in the DOE’s Material Management and Minimization program to convert HEU-fueled research reactors to LEU-fueled reactors. In some demanding conditions in high-power and high-performance reactors, large pores form in the interaction layers between the U-Mo fuel particles and the Al matrix, which pose a potential to cause fuel failure. In this study, comprehension of the formation and growth of these pores was explored. As a product, a model to predict pore growth and porosity increase was developed. Well-characterized in-pile data from reduced-size plates were used to fit the model parameters. A data set of full-sized plates, independent and distinctively different from those used to fit the model parameters, was used to examine the accuracy of the model.

in-pile data↗

Satellite ephemerides for the Voyager Uranus encounter

Uranian satellite ephemerides are needed by the Voyager project to support both navigation and acquisition of scientific data. This paper presents the approach being taken to develop the ephemerides and details the initial phase of the development. That phase involves the analytical modeling of the satellites' motion and the adjustment of the model to fit astronomical observations. The paper describes the model and gives the result of a fit to 71 years of observations.

Jacobson, R. A.↗

Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models

A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.

phenological fitting↗

f-Orbital Mixing in the Octahedral f 2 Compounds UX 6 2– [X = F, Br, Cl, I] and PrCl 6 3–

Understanding how interactions between the f orbitals and ligand orbitals in lanthanide and actinide systems affect their physical properties is the central issue in f-element chemistry. A wide variety of approaches including both theoretical and experimental tools have been used to study these relationships. Among the most widely used tools has been crystal field theory (CFT), which bridges theory and experiment in that it is a model based largely on atomic theory that is parametrized using experimental data. Crystal field theory is quite accurate for the lanthanides, due in part to the highly contracted nature of the 4f orbitals. For actinides, crystal field theory is less accurate, potentially due to the treatment of orbital mixing. In CFT, orbital mixing is handled implicitly by allowing the electron repulsion parameters (Slater F k parameters) and the spin–orbit coupling constant to vary. As a result, orbital mixing in CFT is isotropic in that the F k parameters and the spin–orbit coupling constant affect all f orbitals equally. This approximation works well for the lanthanides due to the limited degree of orbital mixing in these complexes. In actinide complexes, the 5f orbitals have greater overlap with the ligand orbitals, and this approximation is less accurate than in the lanthanides. In this work, we report a modification of CFT that includes the effect of orbital mixing on electron repulsion and spin–orbit coupling for each f orbital. The model is applied to the tetravalent uranium hexahalide dianions and PrCl 6 3– for which the energies of many low-lying excited states are known. The new model generally fits the data as well the traditional CFT although with fewer parameters. However, the new model does not fit the data better than the more complex CFT models of Faucher and co-workers. The results of the model show in detail how changes in overlap and orbital energies influence the energies of the bonding and antibonding orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

WMAP Detection of Early Reionization

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systenatic errors. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales the WMAP data show the temperature-polarization correlation expected from adiabatic density perturbations. The data for $\ell > 20$ agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. The existence of correlations on super-horizon scales provides significant support for inflationary cosmologies. We detect excess power on large angular scales compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift $11 < z < 30$ at 95\% confidence. A model-independent fit to reionization optical depth yields results consistent with the best-fit $\Lambda$CDM model, with best fit value $\tau = 0.17 \pm 0.04$ at 68\% confidence, including systematic and foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history. I will discuss the WMAP data and its implications for reionization in the early universe.

Kogut, A.↗

WMAP Detection of Early Reionization

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systematic errors. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales the WMAP data show the temperature-polarization correlation expected from adiabatic density perturbations. The data for $\ell > 20$ agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. The existence of correlations on super-horizon scales provides significant support for inflationary cosmologies. We detect excess power on large angular scales compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift $11 < z < 30$ at 95\% confidence. A model-independent fit to reionization optical depth yields results consistent with the best-fit $\Lambda$CDM model, with best fit value $\tau = 0.17 \pm 0.04$ at 68\% confidence, including systematic and foreground uncertainties. This trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history. I will discuss the WMAP data and its implications for reionization in the early universe.

Kogut, A.↗

Measurement of the Se 80 ( γ , n ) reaction with linearly polarized γ rays

This paper reports the results of the first measurements of the differential cross section of the 80 Se(γ , n) 79 Se reaction with a linearly polarized gamma-ray (γ-ray) beam. The cross section was measured at three incident γ-ray beam energies: 15.6, 15.8, and 16.0 MeV, with a beam energy spread of 3.0% full width at half-maximum (FWHM). The differential cross section for the excitation spectrum in 79Se was measured at two scattering angles in the plane of the beam polarization: θ = 90° and 135°, and at one angle in the plane perpendicular to the plane of polarization: θ = 90°. The total photoneutron cross sections determined from these data are between 0.8 and 1.3 standard deviations smaller than previously published results. Here, the excitation spectra measured in this work were fit with a Hauser-Feshbach model. Better fits to the data were obtained with a constant-temperature formulation of the nuclear level density (NLD) than with a Fermi-gas NLD model. The parameters for the constant-temperature NLD model obtained in this work are consistent with those obtained for medium-mass nuclei in previous studies.

59 ≤ A ≤ 89↗

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗

Modeling and Maximum Likelihood Fitting of Gamma-Ray and Radio Light Curves of Millisecond Pulsars Detected with Fermi

Pulsed gamma rays have been detected with the Fermi Large Area Telescope (LAT) from more than 20 millisecond pulsars (MSPs), some of which were discovered in radio observations of bright, unassociated LAT sources. We have fit the radio and gamma-ray light curves of 19 LAT-detected MSPs in the context of geometric, outermagnetospheric emission models assuming the retarded vacuum dipole magnetic field using a Markov chain Monte Carlo maximum likelihood technique. We find that, in many cases, the models are able to reproduce the observed light curves well and provide constraints on the viewing geometries that are in agreement with those from radio polarization measurements. Additionally, for some MSPs we constrain the altitudes of both the gamma-ray and radio emission regions. The best-fit magnetic inclination angles are found to cover a broader range than those of non-recycled gamma-ray pulsars.

Johnson, T. J.↗

Understanding the extreme luminosity of DES14X2fna

We present DES14X2fna, a high-luminosity, fast-declining Type IIb supernova (SN IIb) at redshift z = 0.0453, detected by the Dark Energy Survey (DES). DES14X2fna is an unusual member of its class, with a light curve showing a broad, luminous peak reaching M r ≃ –19.3 mag 20 d after explosion. This object does not show a linear decline tail in the light curve until ≃60 d after explosion, after which it declines very rapidly (4.30 ± 0.10 mag 100 d –1 in the r band). By fitting semi-analytic models to the photometry of DES14X2fna, we find that its light curve cannot be explained by a standard 56 Ni decay model as this is unable to fit the peak and fast tail decline observed. Inclusion of either interaction with surrounding circumstellar material or a rapidly-rotating neutron star (magnetar) significantly increases the quality of the model fit. Furthermore, we also investigate the possibility for an object similar to DES14X2fna to act as a contaminant in photometric samples of SNe Ia for cosmology, finding that a similar simulated object is misclassified by a recurrent neural network (RNN)-based photometric classifier as an SN Ia in ~1.1–2.4 per cent of cases in DES, depending on the probability threshold used for a positive classification.

79 ASTRONOMY AND ASTROPHYSICS↗

Lithospheric bending at subduction zones based on depth soundings and satellite gravity

A global study of trench flexure was performed by simultaneously modeling 117 bathymetric profiles (original depth soundings) and satellite-derived gravity profiles. A thin, elastic plate flexure model was fit to each bathymetry/gravity profile by minimization of the L(sub 1) norm. The six model parameters were regional depth, regional gravity, trench axis location, flexural wavelength, flexural amplitude, and lithospheric density. A regional tilt parameter was not required after correcting for age-related trend using a new high-resolution age map. Estimates of the density parameter confirm that most outer rises are uncompensated. We find that flexural wavelength is not an accurate estimate of plate thickness because of the high curvatures observed at a majority of trenches. As in previous studies, we find that the gravity data favor a longer-wavelength flexure than the bathymetry data. A joint topography-gravity modeling scheme and fit criteria are used to limit acceptable parameter values to models for which topography and gravity yield consistent results. Even after the elastic thicknesses are converted to mechanical thicknesses using the yield strength envelope model, residual scatter obscures the systematic increase of mechanical thickness with age; perhaps this reflects the combination of uncertainties inherent in estimating flexural wavelength, such as extreme inelastic bending and accumulated thermoelastic stress. The bending moment needed to support the trench and outer rise topography increases by a factor of 10 as lithospheric age increases from 20 to 150 Ma; this reflects the increase in saturation bending moment that the lithosphere can maintain. Using a stiff, dry-olivine rheology, we find that the lithosphere of the GDH1 thermal model (Stein and Stein, 1992) is too hot and thin to maintain the observed bending moments. Moreover, the regional depth seaward of the oldest trenches (approximately 150 Ma) exceeds the GDH1 model depths by about 400 m.

Levitt, Daniel A.↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗

Determining neutron star properties by fitting oblate-star waveform models to X-ray burst oscillations

We describe sophisticated new Bayesian analysis methods that make it possible to estimate quickly the masses and radii of rapidly rotating, oblate neutron stars by fitting oblate-star waveform models to energy-resolved observations of the waveforms of X-ray burst oscillations produced by such stars. We find that a 25% variation of the temperature of the hot spot in the north-south direction does not significantly bias estimates of the mass M and equatorial radius R{sub eq} derived by fitting a model that assumes a uniform-temperature spot. Our results show that fits of oblate-star waveform models to waveform data can simultaneously determine M and R{sub eq} with 1σ uncertainties ≲7% if (1) the star's rotation rate is ≳ 600 Hz; (2) the spot center and observer's sightline are both within 30{sup ∘} of the star's rotational equator; and (3) ∼10{sup 7} counts are collected during burst oscillations that have a fractional rms amplitude ≳10%. A fractional amplitude of ∼10% is realistic, and the accepted NICER and proposed LOFT and AXTAR space missions could collect 10{sup 7} counts from a single star by combining data from many X-ray bursts from the star. Uncertainties ≲7% are small enough to improve substantially our understanding of cold, ultradense matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Long-term sinking of nuclear waste canisters in salt formations by low-stress creep at high temperature

Abstract Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of $$-4.4\times 10^{-7}$$ - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus $$-2.1\times 10^{-2}$$ - 2.1 × 10 - 2 mm/year and $$-3.1\times 10^{-2}$$ - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗