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 415 records · Page 23

Extracting galactic structure parameters from multivariated density estimation

Multivariate statistical analysis, including includes cluster analysis (unsupervised classification), discriminant analysis (supervised classification) and principle component analysis (dimensionlity reduction method), and nonparameter density estimation have been successfully used to search for meaningful associations in the 5-dimensional space of observables between observed points and the sets of simulated points generated from a synthetic approach of galaxy modelling. These methodologies can be applied as the new tools to obtain information about hidden structure otherwise unrecognizable, and place important constraints on the space distribution of various stellar populations in the Milky Way. In this paper, we concentrate on illustrating how to use nonparameter density estimation to substitute for the true densities in both of the simulating sample and real sample in the five-dimensional space. In order to fit model predicted densities to reality, we derive a set of equations which include n lines (where n is the total number of observed points) and m (where m: the numbers of predefined groups) unknown parameters. A least-square estimation will allow us to determine the density law of different groups and components in the Galaxy. The output from our software, which can be used in many research fields, will also give out the systematic error between the model and the observation by a Bayes rule.

Chen, B.↗

Depth dependence of absorbed dose, dose equivalent and linear energy transfer spectra of galactic and trapped particles in polyethylene and comparison with calculations of models

A matched set of five tissue-equivalent proportional counters (TEPCs), embedded at the centers of 0 (bare), 3, 5, 8 and 12-inch-diameter polyethylene spheres, were flown on the Shuttle flight STS-81 (inclination 51.65 degrees, altitude approximately 400 km). The data obtained were separated into contributions from trapped protons and galactic cosmic radiation (GCR). From the measured linear energy transfer (LET) spectra, the absorbed dose and dose-equivalent rates were calculated. The results were compared to calculations made with the radiation transport model HZETRN/NUCFRG2, using the GCR free-space spectra, orbit-averaged geomagnetic transmission function and Shuttle shielding distributions. The comparison shows that the model fits the dose rates to a root mean square (rms) error of 5%, and dose-equivalent rates to an rms error of 10%. Fairly good agreement between the LET spectra was found; however, differences are seen at both low and high LET. These differences can be understood as due to the combined effects of chord-length variation and detector response function. These results rule out a number of radiation transport/nuclear fragmentation models. Similar comparisons of trapped-proton dose rates were made between calculations made with the proton transport model BRYNTRN using the AP-8 MIN trapped-proton model and Shuttle shielding distributions. The predictions of absorbed dose and dose-equivalent rates are fairly good. However, the prediction of the LET spectra below approximately 30 keV/microm shows the need to improve the AP-8 model. These results have strong implications for shielding requirements for an interplanetary manned mission.

Flight Experiment↗

The resolution of P Cygni's stellar wind

The 6 cm radio emission from P Cygni's stellar wind has been resolved using the Very Large Array. The wind is apparently spherically symmetric, and the radial surface brightness distribution agrees with that predicted by the constant velocity, isothermal wind model. The wind temperature is directly determined from the observations to be within 2000 K of 18,000, which is consistent with an effective temperature of 20,000 K for P Cygni. The technique makes it possible to study the adiabatic cooling of the wind, which should be important for most early-type stars but which cannot be recognized through its effects on the radio spectrum. To avoid systematically underestimating fluxes, it is essential to fit models directly to the (u-v) data when observing stellar winds with the VLA.

White, R. L.↗

The Three-Dimensional Expansion of the Ejecta from Tycho's Supernova Remnant

We present the first 3D measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12 yr baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 tufts of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line-of-sight velocity, we use two different methods: a nonequilibrium ionization model fit to the strong Si and S lines in the 1.22.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the redshift or blueshift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km/s, with a mean of 4430 km/s. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km/s. Some SN Ia explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and we discuss our findings in light of various explosion models, favoring those delayed-detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

ISM: individual objects (Tycho’s SNR) – ISM: s↗

Systematic Matter and Binding-Energy Distributions from a Dispersive Optical Model Analysis

We present the first systematic dispersive-optical-model analysis of 16,18 O, 40,48 Ca, 58,64 Ni, 112,124 Sn, and 208 Pb. In all systems, roughly half the total nuclear binding energy is associated with the most-bound 10% of the total nucleon density. The extracted neutron skins reveal interplay of asymmetry, Coulomb, and shell effects on the skin thickness. Here, our results indicate that optical-model fits on isotopic pairs are effective for constraining nuclear structural quantities otherwise difficult to observe experimentally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ASCA and Radio/RXTE Observations of GX 339-4

We have analyzed three separate archival Advanced Satellite for Cosmology and Astrophysics (ASCA) observations and eight separate Rossi X-ray Timing Explorer (RXTE) observations of the black hole candidate GX 339-4 in its low luminosity, spectrally hard state. Three of the RXTE observations were strictly simultaneous with 843 MHz and 8.3- 9.1 GHz radio observations. All data sets show evidence for an approximately 6.4 keV Fe line with equivalent widths approximately 20-100 eV. 'Reflection models' show a hardening of the RXTE spectra with decreasing X-ray flux; however, these models do not exhibit evidence of a correlation between the photon index of the incident power law flux and the solid angle subtended by the reflector. None of the models fit to the X-ray data, however, simultaneously explain the observed radio properties. We argue that the spatial extent of the observed radio emission is at least 0(10(exp 7 GM/c squared). Timing analysis reveals that all observations save one show evidence of a persistent f(qpo approximately equals 0.3 Hz quasi-periodic oscillations(quasi-periodic oscillations)). The broad band (10-3-102 Hz) power appears to be dominated by two independent processes that can be modeled as very broad Lorentzians with Q approximately less than 1. Similar to Cyg X-1, the hard photon variability is seen to lag the soft photon vaxiability with the lag time increasing with decreasing Fourier frequency. The magnitude of this time lag is seen to be positively correlated with the flux of GX 339-4.

Nowak, M. A.↗

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Evaluation of PV Module Packaging Strategies of Monofacial and Bifacial PERC Using Degradation Pathway Network Modeling

As the PV industry is rapidly expanding, it is important to thoroughly investigate the long-term impact of packaging strategies on the performance of PV modules. In this study, the variants in sets differ on the basis of manufacturer (A/B), encapsulant (EVA/POE), rear encapsulant (UV-cutoff/opaque/transparent), module architecture (GB/DG) and cell type (monofacial/bifacial). The minimodules were exposed for 2520 hours in modified damp heat, with or without full spectrum light. Every 504 hours, stepwise electrical characterization techniques were employed to track changes in minimodules. Degradation pathway modeling using network structural equation modeling was employed to study pairwise relationships between variables and service lifetime prediction in minimodules. Through this study, differences in quality control are identified in minimodules made by different manufacturers. Minimodules with UV-cutoff rear encapsulant show relatively better stability, whereas the ones with opaque rear encapsulant show greater power loss. In addition, GB having UV-cutoff rear encapsulation and GB with POE having opaque rear encapsulation were identified to be stable as they lack a best model fit. Here, the primary power loss mechanism in degrading variants is interconnect corrosion.

14 SOLAR ENERGY↗

Thermal-model-based characterization of heavy-ion-collision systems at chemical freeze-out

We investigate the chemical freeze-out in heavy-ion collisions (HICs) and the impact of the hadronic spectrum on thermal model analyses [1, 2]. Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of HICs. By varying the number of resonances included in Hadron Resonance Gas (HRG) Model calculations, we can shed light on which particles may be produced. Furthermore, we study the influence of the number of states on the so-called two flavor freezeout scenario, in which strange and light particles can freeze-out separately. We consider results for the chemical freeze-out parameters obtained from thermal model fits and from calculating net-particle fluctuations. We will show the effect of using one global temperature to fit all particles and alternatively, allowing particles with and without strange quarks to freeze-out separately.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calibration of Al7075 with Plate-Puncture Predictions

The following details calibration of a material model for Al7075-T6511. This aluminum alloy is commonly used across a host of engineering applications. Owing to its widespread prevalence, there is great benefit in improving simulation predictions for this alloy. In the present effort, a calibration is performed of its elastic-plastic response accounting for both rate and temperature dependence. The calibration is informed by a series of tests that include specimens of different geometries tested at different rates and temperatures. All specimens are derived from the same barstock, 3.5 inches in diameter. The fitted model itself uses an anisotropic, Hill yield surface coupled with a Johnson-Cook hardening model. Failure predictions are had by means of a modified Wilkins failure criterion. Following calibration of the material model, a validation exercise is performed against platepuncture experiments. These experiments include multiple probe shapes, probe diameters, and plate thicknesses. The puncture experiments are replicated in simulation with mesh studies performed to assess uncertainty. Key quantities of interest, notably the absorbed energy up to failure, are compared between simulation and experiment providing a means to assess the suitability of the calibration in puncture simulations.

36 MATERIALS SCIENCE↗

Superluminal motion in the double-lobed quasar 3C 245

VLBI observations of 3C 245 obtained at 10.7 GHz with the Mark III array at five epochs in 1981-1986 are reported. The data-reduction and model-fitting procedures are explained, and the results are presented in tables, graphs, and maps along with previously published data on 11 other superluminal sources. The central component of the quasar is found to have an apparent transverse velocity of (3.1 + or - 1.4) c/h, assuming H0 = 100h km/s Mpc and q0 = 0.5. From the apparent source orientation (based on R, variability, apparent curvature, and projected linear size) it is inferred that the apparent transverse velocity is anticorrelated with the angle to the line of sight, as predicted by relativistic-beaming models of superluminal motion.

Hough, D. H.↗

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗

X-ray Emission from the Pre-planetary Nebula Henize 3-1475

We report the first detection of X-ray emission in a pre-planetary nebula, He 3-1475. Pre-planetary nebulae are rare objects in the short transition stage between the asymptotic giant branch (AGB) and planetary nebula evolutionary phases, and He 3-1475, characterized by a remarkable S-shaped chain of optical knots, is one of the most noteworthy members of this class. Observations with the Advanced CCD Imaging Spectrometer on board the Chandra X-Ray Observatory show the presence of compact emission coincident with the brightest optical knot in this bipolar object, which is displaced from the central star by 2'.7 along the polar axis. Model fits to the X-ray spectrum indicate an X-ray temperature and luminosity, respectively, of (4.3-5.7) x 10(exp 6) K and (4 +/- 1.4) x 10(exp 31) (D/5 kpc)(exp 2) ergs s(exp -1) respectively. Our 3 sigma upper limit on the luminosity of compact X-ray emission from the central star in He 3-1475 is approximately equal to 5 x 10(exp 31) (D/5 kpc)(exp 2) ergs s(exp -1). The detection of X-rays in He 3-1475 is consistent with models in which fast collimated post-AGB outflows are crucial to the shaping of nebulae; we discuss such models in the context of our observations.

circumstellar matter↗

The clustering of DESI-like luminous red galaxies using photometric redshifts

ABSTRACT We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging data set, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 deg2. We have developed new photometric redshift (photo-z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalogue, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogues, and present a new technique for finding best-fitting parameters in this situation. Accompanying this paper, we have released the Photometric Redshifts for the Legacy Surveys catalogue of photo-z’s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 data set. This catalogue provides accurate photometric redshifts for objects with z < 21 over more than 16 000 deg2 of sky.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗

The Swift Burst Alert Telescope Detected Seyfert 1 Galaxies: X-Ray Broadband Properties and Warm Absorbers

We present results from an analysis of the broadband, 0.3-195 keV, X-ray spectra of 48 Seyfert 1-1.5 sources detected in the very hard X-rays with the Swift Burst Alert Telescope (BAT). This sample is selected in an all-sky survey conducted in the 14-195 keV band. Therefore, our sources are largely unbiased toward both obscuration and host galaxy properties. Our detailed and uniform model fits to Suzaku/BAT and XMM-Newton/BAT spectra include the neutral absorption, direct power-law, reflected emission, soft excess, warm absorption, and narrow Fe I K[alpha] emission properties for the entire sample. We significantly detect O VII and O VIII edges in 52% of our sample. The strength of these detections is strongly correlated with the neutral column density measured in the spectrum. Among the strongest detections, X-ray grating and UV observations, where available, indicate outflowing material. The ionized column densities of sources with O VII and O VIII detections are clustered in a narrow range with Nwarm [approx] 1021 cm-2, while sources without strong detections have column densities of ionized gas an order of magnitude lower. Therefore, we note that sources without strong detections likely have warm ionized outflows present but at low column densities that are not easily probed with current X-ray observations. Sources with strong complex absorption have a strong soft excess, which may or may not be due to difficulties in modeling the complex spectra of these sources. Still, the detection of a flat [Gamma] [approx] 1 and a strong soft excess may allow us to infer the presence of strong absorption in low signal-to-noise active galactic nucleus spectra. Additionally, we include a useful correction from the Swift BAT luminosity to bolometric luminosity, based on a comparison of our spectral fitting results with published spectral energy distribution fits from 33 of our sources.

Winter, Lisa M.↗

The European Southern Observatory-MIDAS table file system

The new and substantially upgraded version of the Table File System in MIDAS is presented as a scientific database system. MIDAS applications for performing database operations on tables are discussed, for instance, the exchange of the data to and from the TFS, the selection of objects, the uncertainty joins across tables, and the graphical representation of data. This upgraded version of the TFS is a full implementation of the binary table extension of the FITS format; in addition, it also supports arrays of strings. Different storage strategies for optimal access of very large data sets are implemented and are addressed in detail. As a simple relational database, the TFS may be used for the management of personal data files. This opens the way to intelligent pipeline processing of large amounts of data. One of the key features of the Table File System is to provide also an extensive set of tools for the analysis of the final results of a reduction process. Column operations using standard and special mathematical functions as well as statistical distributions can be carried out; commands for linear regression and model fitting using nonlinear least square methods and user-defined functions are available. Finally, statistical tests of hypothesis and multivariate methods can also operate on tables.

Peron, M.↗

Statistical studies of impulsive events at high latitudes

A statistical study has been made of the high-latitude impulsive events that were observed during the 1985-1986 South Pole Balloon Campaign. The events were selected by searching for unipolar pulses greater than or equal to 10 nT above background in the vertical component of the magnetic field on the ground and/or pedestal or 'W' shaped horizontal electric field perturbations greater than or equal to 10 mV/m in amplitude and accompanied by perturbations in the vertical electric field at balloon altitude. A main event list comprising 112 events was compiled from the 468 hours of data available. Three aspects of the events were examined: the solar wind conditions prior to the event, local time of observation, and intrinsic properties of the events. The local time distribution was obtained from the 112 entry main event list and was found to be nearly uniform across the dayside, with no midday gap. The event rate found using our low-amplitude selection criteria was 0.7 event/hr, comparable to expectations based on in situ studies of the magnetopause. A total of 42 events were found for which data were available from Interplanetary Monitoring Platform (IMP) 8. Of these events, 12 occurred when the Z(sub GSM) component (B(sub Z)) of the interplanetary magnetic field (IMF) was northward and 30 occurred when B(sub Z) was southward or fluctuating. Only three of the B(sub Z) northward cases and only five of the B(sub Z) southward cases were preceded by pressure pulses greater than 0.4 nPa in amplitude. Ten of the events were studied in detail by means of a model-fitting method discussed elsewhere. This method infers values of several parameters, including the total current flowing in a coaxial or monopole system and a two-dimensional dipole system. The intrinsic properties of the events showed that only approximately 10% of the total current contributed to momentum transfer to the high-latitude ionosphere, that the direction of the motion depended more on local time of observation than IMF B(sub y), and that events were usually several hundred kilometers in size. The observed B(sub z) control found in the 42 event list and the prevalence of coaxial current dominated events are inconsistent with the predictions of the pressure pulse model.

Lin, Z. M.↗