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 199 records · Page 11

High-contrast Thermal Infrared Spectroscopy with ALES: The 3–4 μm Spectrum of κ Andromedae b

We present the first L-band (2.8–4.1 μm) spectroscopy of κ Andromedae b, a ~20 M Jup companion orbiting at 1" projected separation from its B9-type stellar host. We combine our Large Binocular Telescope (LBT) Arizona Lenslets for Exoplanet Spectroscopy (ALES) integral field spectrograph data with measurements from other instruments to analyze the atmosphere and physical characteristics of κ And b. We report a discrepancy of ~20% (2σ) in the L' flux of κ And b when comparing to previously published values. We add an additional L' constraint using an unpublished imaging data set collected in 2013 using the LBT Interferometer/LMIRCam, the instrument in which the ALES module has been built. The LMIRCam measurement is consistent with the ALES measurement, both suggesting a fainter L-band scaling than previous studies. The data, assuming the flux scaling measured by ALES and LMIRCam imaging, are well fit by an L3-type brown dwarf. Atmospheric model fits to measurements spanning 0.9–4.8 μm reveal some tension with the predictions of evolutionary models, but the proper choice of cloud parameters can provide some relief. In particular, models with clouds extending to very low pressures composed of grains ≤1 μm appear to be necessary. If the brighter L' photometry is accurate, there is a hint that subsolar metallicity may be required.

79 ASTRONOMY AND ASTROPHYSICS↗

A SPECTROSCOPIC CATALOG OF THE BRIGHTEST (J < 9) M DWARFS IN THE NORTHERN SKY{sup ,}

We present a spectroscopic catalog of the 1564 brightest (J < 9) M dwarf candidates in the northern sky, as selected from the SUPERBLINK proper motion catalog. Observations confirm 1408 of the candidates to be late-K and M dwarfs with spectral subtypes K7-M6. From the low ({mu} > 40 mas yr{sup -1}) proper motion limit and high level of completeness of the SUPERBLINK catalog in that magnitude range, we estimate that our spectroscopic census most likely includes >90% of all existing, northern-sky M dwarfs with apparent magnitude J < 9. Only 682 stars in our sample are listed in the Third Catalog of Nearby Stars (CNS3); most others are relative unknowns and have spectroscopic data presented here for the first time. Spectral subtypes are assigned based on spectral index measurements of CaH and TiO molecular bands; a comparison of spectra from the same stars obtained at different observatories, however, reveals that spectral band index measurements are dependent on spectral resolution, spectrophotometric calibration, and other instrumental factors. As a result, we find that a consistent classification scheme requires that spectral indices be calibrated and corrected for each observatory/instrument used. After systematic corrections and a recalibration of the subtype-index relationships for the CaH2, CaH3, TiO5, and TiO6 spectral indices, we find that we can consistently and reliably classify all our stars to a half-subtype precision. The use of corrected spectral indices further requires us to recalibrate the {zeta} parameter, a metallicity indicator based on the ratio of TiO and CaH optical bandheads. However, we find that our {zeta} values are not sensitive enough to diagnose metallicity variations in dwarfs of subtypes M2 and earlier ({+-}0.5 dex accuracy) and are only marginally useful at later M3-M5 subtypes ({+-}0.2 dex accuracy). Fits of our spectra to the Phoenix atmospheric model grid are used to estimate effective temperatures. These suggest the existence of a plateau in the M1-M3 subtype range, in agreement with model fits of infrared spectra but at odds with photometric determinations of T{sub eff}. Existing geometric parallax measurements are extracted from the literature for 624 stars, and are used to determine spectroscopic and photometric distances for all the other stars. Active dwarfs are identified from measurements of H{alpha} equivalent widths, and we find a strong correlation between H{alpha} emission in M dwarfs and detected X-ray emission from ROSAT and/or a large UV excess in the GALEX point source catalog. We combine proper motion data and photometric distances to evaluate the (U, V, W) distribution in velocity space, which is found to correlate tightly with the velocity distribution of G dwarfs in the solar neighborhood. However, active stars show a smaller dispersion in their space velocities, which is consistent with those stars being younger on average. Our catalog will be most useful to guide the selection of the best M dwarf targets for exoplanet searches, in particular those using high-precision radial velocity measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

CONSTRAINT-INDEPENDENT CONSTANT CTOA DETERMINATION FOR DUCTILE STABLE CRACK GROWTH

Crack tip opening angle (CTOA) has been used as a reliable fracture toughness parameter for decades to characterize stable ductile crack growth for thin-walled aerospace structures in the low-constraint conditions. Recently, the CTOA parameter was also applied to the pipeline industry, and a CTOA test standard ASTM E3039 was thus developed for testing a critical constant CTOA. Research showed that the constant CTOA can reasonably describe fracture toughness required to arrest a dynamic crack propagation for a modern gas pipeline. However, the CTOA fracture criterion requires constraint-independent CTOA toughness against stable ductile crack growth. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen for CTOA testing. Since a shallow crack is used, DWTT measured CTOA may depend on constraint level at the crack tip. To understand if it is the case, this paper evaluates the critical CTOA for a set of fracture toughness tests on single edge notched bend (SENB) specimens with shallow and deep cracks based on four CTOA estimation models. In which, the Ln(P)-LLD linear fit model is similar to that used by ASTM E3039 in the CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel are considered in the CTOA evaluation. The results show that the four CTOA models can determine a crack size-independent constant CTOA over stable ductile crack growth for the SENB specimens. As a result, CTOA determined by ASTM E3039 is constraint&#x2;independent and transferable to use for an actual crack propagating in a gas pipeline.

Zhu, Xian-Kui↗

Kinetic Study of Coal and Biomass Co-Pyrolysis Using Thermogravimetry

The objectives of this study are to investigate thermal behavior of coal and biomass blends in inert gas environment at low heating rates and to develop a simplified kinetic model using model fitting techniques based on TGA experimental data. Differences in thermal behavior and reactivity in co-pyrolysis of Powder River Basin (PRB) sub-bituminous coal and pelletized southern yellow pine wood sawdust blends at low heating rates are observed. Coal/wood blends have higher reactivity compared to coal alone in the lower temperature due to the high volatile matter content of wood. As heating rates increase, weight loss rates increase. The experiment data obtained from TGA has a better fit with proposed two step first order reactions model compared single first order reaction model.

01 COAL, LIGNITE, AND PEAT↗

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

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 glasses for flowsheet evaluation, testing, and design of the High-Level Waste Facility at the Hanford Waste Treatment and Immobilization Plant. These models and constraints are meant to be used as a placeholder 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, and 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↗

Dark Energy Survey Year 6 Results: Point-Spread Function Modeling

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗

First Principle based Thermal Conductivity Modeling of Metal Alloys with Applications on Uranium Alloys

Thermal conductivity is an important materials property related to heat transport, which is essential to many applications, ranging from thermoelectrics to nuclear reactor materials. High- quality thermal conductivity data is critical to these materials and their associated technologies, including the nuclear fuel materials like traditional oxide fuels and metallic uranium (U) fuels, e.g. U-Zr and U-Mo alloys. And thermal conductivity modeling is widely used to interpolate or extrapolate experimental data, to obtain high-quality thermal conductivity data over wide temperature and composition range, and to understand the impacts of different factors including defects and microstructures. However, a general thermal conductivity model for metal alloys which can work regrading different phase components is still missing. Also, for U alloys, a thermal conductivity model working with different phases is needed. Therefore, in this work, we developed thermal conductivity models for metal alloys, based on ab-initio calculations, semi-classical physics rules, and limited experimental data. The goals of these models are to help obtain high- quality thermal conductivity data within a little experimental input as possible, have models that are extendable to varying microstructures, including different types of irradiation effects, and provide mechanistic understanding of heat transfer in the modeled alloys. In this work, our model solves several challenges in the development. The DFT-BTE approach is applied to decrease the reliance on experimental data and make the model extendable to composition changes and defects. A practical and efficient way to combine the DFT inputs with physics rules is pointed out in this work too. A staged approach, which starts from simple cases of elemental metals, then extends to solid solutions, different compound phases, and multi-phase mixtures, is presented to work with the metal alloys in different phase components. Our models are demonstrated respectively on aU for elemental metal model, on U-Zr and U-Mo alloys in aU temperature range for multi-phase mixture model, on high-temperature U-Zr, U-Nb, and U-Mo alloys in the body-centered cubic phase for concentrated solid solution model, and on irradiated U-Mo alloys for irradiated metal alloy model. All models show great agreement with experimental data. In these demonstrations, our model shows its advantages compared to previous empirical fitting model. Our model requires fewer experimental data, due to the inputs from DFT. The quantitative insights into the different physical factors are provided in our model, as it incorporates the electron and phonon scattering mechanisms. Our model also can be extended to incorporate the contributions of point defects, grain boundaries, and noble gas bubbles, to integrate their effects on heat transfer. This model can serve as both a foundation for understanding the more complex thermal conductivity of realistic U alloy fuels and a useful tool to guide further modeling of thermal conductivity to aid materials and device design.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A New Method for the Evaluation of Majority and Minority Carrier Contact Resistivity of Polysilicon on Oxide Contacts

In this contribution we present a novel technique for characterizing the majority and minority carrier contact resistivity in polysilicon on oxide passivating contacts. Traditionally, contact resistivity is measured through the Transfer Length Method (TLM). This method presents difficulties when measuring the minority carrier contact due to the presence of a blocking diode which prevents current flow between the pads. Our novel method is performed by creating a unipolar device which has the same contact on both sides (either symmetric n-type or symmetric p-type) and metalizing the device with a grid on one side, and a full metal contact on the back (the same pattern that will be used in a final cell). We then measure the resistance across the device at increasing illumination. Since the resistivity of the wafer is a function of carrier density, when increasing illumination, the wafer resistivity decreases. We can then fit modeling equations to the measured resistivity and extrapolate the fit to infinite illumination where only the contact resistivity is measured to accurately measure both the minority and majority carrier contact resistivity.

current measurement↗

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The radial profile and flattening of the Milky Way’s stellar halo to 80 kpc from the SEGUE K-giant survey

We characterize the radial density, metallicity, and flattening profile of the Milky Way’s stellar halo, based on the large sample of spectroscopically confirmed giant stars from SDSS/SEGUE-2, spanning galactocentric radii 10 kpc ≤ r{sub GC} ≤ 80 kpc. After excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream), the sample has 1757 K giants, with a typical metallicity precision of 0.2 dex and a mean distance accuracy of 16%. Compared to blue horizontal branch stars or RR Lyrae variables, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modeling of those data, based on ellipsoidal stellar density models, ν{sub *}(R, z), with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of n = 3.1 ± 0.5 with an effective radius r{sub eff}=15±2 kpc and a flattening of q = 0.7 ± 0.02, or comparably well by an equally flattened broken power law, with radial slopes of α{sub in} = 2.1 ± 0.3 and α{sub out} = 3.8 ± 0.1, with a break radius of r{sub break} = 18 ± 1 kpc; this is largely consistent with earlier work. We find a modest but significant metallicity gradient within the “outer” stellar halo, [Fe/H] decreasing outward. If we allow for a variable flattening q=f(r{sub GC}), we find the distribution of halo giants to be considerably more flattened at small radii, q(10 kpc) = 0.55 ± 0.02, compared to q(>30 kpc) = 0.8 ± 0.03. Remarkably, the data are then very well fit by a single power law with index of 4.2 ± 0.1 on the variable r{sub q}≡√(R{sup 2}+(z/q(r)){sup 2}). In this simple and better-fitting model, there is a break in flattening at ∼20 kpc, instead of a break in the radial density function. While different parameterizations of the radial profile vary in their parameters, their implied density gradient, ∂lnν{sub ∗}/∂lnr, is stable along a direction intermediate between major and minor axis; this gradient is crucial in any dynamical modeling that uses halo stars as tracers.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic scale simulations for improved CRUD and fuel performance modeling

A more mechanistic description of fuel performance codes can be achieved by deriving models and parameters from atomistic scale simulations rather than fitting models empirically to experimental data. The same argument applies to modeling deposition of corrosion products on fuel rods (CRUD). Here are some results from publications in 2016 carried out using the CASL allocation at LANL.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field

We have conducted a study of star formation in the outer Galaxy from 65°< l < 265° in the region observed by the GLIMPSE360 program. This Spitzer warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5 μm. We combine the IRAC, Wide-field Infrared Survey Explorer (WISE), and Two Micron All Sky Survey catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 young stellar objects (YSOs) across the field spanning >180° in Galactic longitude. Using the DBSCAN method on the combined catalog, we identify 618 clusters or aggregations of YSOs having five or more members. We identify 10,476 class I, 29,604 class II, and 7325 anemic class II/class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found that 100 of the clusters identified have previously measured distances in the WISE H ii survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with <3σ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters; dividing the clusters by galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M {sub ⊙} for R {sub Gal} < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M {sub ⊙} for R {sub Gal} > 11.5 kpc. The slope of the combined IMF was found to be Γ = 1.92 ± 0.42 above 3 M {sub ⊙}. These values are consistent with each other within the uncertainties and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

SPITZER TRANSITS OF THE SUPER-EARTH GJ1214b AND IMPLICATIONS FOR ITS ATMOSPHERE

We observed the transiting super-Earth exoplanet GJ1214b using warm Spitzer at 4.5 {mu}m wavelength during a 20 day quasi-continuous sequence in 2011 May. The goals of our long observation were to accurately define the infrared transit radius of this nearby super-Earth, to search for the secondary eclipse, and to search for other transiting planets in the habitable zone of GJ1214. We here report results from the transit monitoring of GJ1214b, including a reanalysis of previous transit observations by Desert et al. In total, we analyze 14 transits of GJ1214b at 4.5 {mu}m, 3 transits at 3.6 {mu}m, and 7 new ground-based transits in the I+z band. Our new Spitzer data by themselves eliminate cloudless solar composition atmospheres for GJ1214b, and methane-rich models from Howe and Burrows. Using our new Spitzer measurements to anchor the observed transit radii of GJ1214b at long wavelengths, and adding new measurements in I+z, we evaluate models from Benneke and Seager and Howe and Burrows using a {chi}{sup 2} analysis. We find that the best-fit model exhibits an increase in transit radius at short wavelengths due to Rayleigh scattering. Pure water atmospheres are also possible. However, a flat line (no atmosphere detected) remains among the best of the statistically acceptable models, and better than pure water atmospheres. We explore the effect of systematic differences among results from different observational groups, and we find that the Howe and Burrows tholin-haze model remains the best fit, even when systematic differences among observers are considered.

79 ASTRONOMY AND ASTROPHYSICS↗

climextremes v0.2.1

Software package for R and Python for analysis of extremes in climate data. Functions for fitting GEV and POT (via point process fitting) models for extremes in climate data, providing return values, return probabilities, and return periods for stationary and nonstationary models. Also provides differences in return values and differences in log return probabilities for contrasts of covariate values. Functions for estimating risk ratios for event attribution analyses, including uncertainty. Under the hood, many of the functions use functions from 'extRemes', including for fitting the statistical models.

Krishnan, Harinarayan↗

In situ synchrotron powder diffraction study of the setting reaction kinetics of magnesium-potassium phosphate cements

This work reports a kinetic study of the formation of magnesium-potassium phosphate cements accomplished using in-situ synchrotron powder diffraction. The reaction: MgO + KH{sub 2}PO{sub 4} + 5H{sub 2}O → MgKPO{sub 4} · 6H{sub 2}O was followed in situ in the attempt of contributing to explain the overall mechanism and assess the influence of periclase (MgO) grain size and calcination temperature (1400-1600 °C) on the reaction kinetics. Numerical kinetic parameters for the setting reaction have been provided for the first time. The best fit to the kinetic data was obtained using a weighted nonlinear model fitting method with two kinetic equations, representing two consecutive, partially overlapping processes. MgO decomposition could be described by a first order (F1) model followed by a Jander diffusion (D3) controlled model. Crystallization of the product of reaction was modelled using an Avrami model (A{sub n}) followed by a first order (F1) chemical reaction. A reaction mechanism accounting for such results has been proposed.

36 MATERIALS SCIENCE↗

Measurements and Simulations of Neutron Emission Versus Deuterium Filling Pressure in Plasma Focus Device PF-24

A series of experiments were carried out using the middle energy (dense) plasma focus device PF-24 with deuterium as a working gas under pressure in a range between 2 and 5 mbar for 17 kV of charging voltage. The relationship between these operating deuterium pressures in reference to the total neutron yield (Y{sub n}) was estimated. The 5-phase Lee code was used to simulate the measured discharge current and neutron yield (Y{sub n}) using a phenomenological beam-target neutron generating mechanism, which was incorporated in the model. Comparison of the Y{sub n} versus pressure using fitted model parameters was made at each point of pressure. The good agreement between measured and computed Y{sub n} values was achieved for discharges with lower neutron emission. The measured Y{sub n} (below 2.6 × 10{sup 9} n/discharge—the median value) has been reproducible by the model for 73% of simulated discharges, while above the median value its prediction were incorrect. The kinetic plasma parameters which were measured and computed using the Lee code for different pressures are: the time to a current sheath collapse (t{sub c}), the average axial current sheath velocity (v{sub z}) and the so called velocity factor (RF). Good agreement was found in the whole range of deuterium pressures between the computed and measured results for these kinematic quantities. Presented findings in this work suggest that the character of neutron emission is more complex than it would seem from classical interpretation of neutron production based on a beam-target model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗