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 217 records · Page 12

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗

Atomic structure of an FeCrMoCBY metallic glass revealed by high energy x-ray diffraction

Amorphous bulk metallic glasses with the composition Fe 48 Cr 15 Mo 14 C 15 B 6 Y 2 have been of interest due to their special mechanical and electronic properties, including corrosion resistance, high yield-strength, large elasticity, catalytic performance, and soft ferromagnetism. Here, we apply a reverse Monte Carlo technique to unravel the atomic structure of these glasses. The pair-distribution functions for various atomic pairs are computed based on the high-energy x-ray diffraction data we have taken from an amorphous sample. Monte Carlo cycles are used to move the atomic positions until the model reproduces the experimental pair-distribution function. The resulting fitted model is consistent with our ab initio simulations of the metallic glass. Our study contributes to the understanding of functional properties of Fe-based bulk metallic glasses driven by disorder effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanistic Modeling of Creep and Aging for Time Extrapolation in LPBF 316 SS - LANL

A polycrystal constitutive model was extended to quantify the effects heterogeneity in the microstructure, porosity, and internal stresses on the mechanical response (creep, tensile) of 316H steel processed by additive manufacturing. The starting point is a recently developed and fully calibrated advanced microstructuresensitive mechanistic constitutive model fit to literature data pertaining to conventionally processed steels. The model is embedded within the full field fast Fourier transform framework (LApx) and several hundred simulations were performed. Synthetic microstructures with different grain aspect ratios were generated using the Neper package to study the role of grain morphology on the mechanical responses. The effect of initial cell interior dislocation density, cell-wall density, precipitate content, porosity, and internal stresses on creep rates are quantified for different polycrystalline microstructures. The study allowed to prospectively quantify the role of distinct features of the microstructure on the mechanical response under tension and creep loading scenarios. Critically it is found that minor anisotropies in the tensile response of samples, emanating from the aspect ratio of grains correlated into significant differences in creep rate as a function of loading directions. The study also allows to identify initial microstructure porosities below which the materials response only moderately degraded by the initial microstructure.

36 MATERIALS SCIENCE↗

Multi-wavelength observations of GRB 111228A and implications for the fireball and its environment

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 s to about 50 hr after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data from the Swift and Fermi missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features a low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow light curves with the standard afterglow model via the Markov Chain Monte Carlo technique shows that ϵ{sub e}=(6.9±0.3)×10{sup −2}, ϵ{sub B}=(7.73±0.62)×10{sup −6},E{sub K}=(6.32±0.86)×10{sup 53} erg, n=0.100±0.014 cm{sup −3}. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. A chromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The E{sub p} of its time-integrated prompt gamma-ray spectrum is ∼26 KeV. Using the initial Lorentz factor (Γ{sub 0}=476{sub −237}{sup +225}) derived from our afterglow model fit, it is found that GRB 111228A satisfies the L{sub iso}−E{sub p,z}−Γ{sub 0} relation and bridges the typical GRBs and low luminosity GRBs in this relation.

79 ASTRONOMY AND ASTROPHYSICS↗

RELICS: A Strong Lens Model for SPT-CLJ0615–5746, a z = 0.972 Cluster

In this work, we present a lens model for the cluster SPT-CLJ0615-5746, which is the highest-redshift ( z = 0.972) system in the Reionization of Lensing Clusters Survey, making it the highest-redshift cluster for which a full, strong lens model is published. We identify three systems of multiply imaged lensed galaxies, two of which we spectroscopically confirm at z = 1.358 and z = 4.013, which we use as constraints for the model. We find a foreground structure at z ~ 0.4, which we include as a second cluster-sized halo in one of our models; however, two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within r = 26 7 (the region where the strong lensing constraints exist) to be M ⊙ . If we extrapolate out to r 500 , our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev–Zel'dovich effect ( M ~ 10 15 M ⊙ ). This cluster is lensing a previously reported z ~ 10 galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

79 ASTRONOMY AND ASTROPHYSICS↗

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Non-local gravity and comparison with observational datasets

We study the cosmological predictions of two recently proposed non-local modifications of General Relativity. Both models have the same number of parameters as ΛCDM, with a mass parameter m replacing the cosmological constant. We implement the cosmological perturbations of the non-local models into a modification of the CLASS Boltzmann code, and we make a full comparison to CMB, BAO and supernova data. We find that the non-local models fit these datasets very well, at the same level as ΛCDM. Among the vast literature on modified gravity models, this is, to our knowledge, the only example which fits data as well as ΛCDM without requiring any additional parameter. For both non-local models parameter estimation using Planck +JLA+BAO data gives a value of H{sub 0} slightly higher than in ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗

NEW CONSTRAINTS ON THE BLACK HOLE LOW/HARD STATE INNER ACCRETION FLOW WITH NuSTAR

We report on an observation of the Galactic black hole candidate GRS 1739–278 during its 2014 outburst, obtained with NuSTAR. The source was captured at the peak of a rising ''low/hard'' state, at a flux of ∼0.3 Crab. A broad, skewed iron line and disk reflection spectrum are revealed. Fits to the sensitive NuSTAR spectra with a number of relativistically blurred disk reflection models yield strong geometrical constraints on the disk and hard X-ray ''corona''. Two models that explicitly assume a ''lamp post'' corona find its base to have a vertical height above the black hole of h=5{sub −2}{sup +7} GM/c{sup 2} and h = 18 ± 4 GM/c {sup 2} (90% confidence errors); models that do not assume a ''lamp post'' return emissivity profiles that are broadly consistent with coronae of this size. Given that X-ray microlensing studies of quasars and reverberation lags in Seyferts find similarly compact coronae, observations may now signal that compact coronae are fundamental across the black hole mass scale. All of the models fit to GRS 1739–278 find that the accretion disk extends very close to the black hole—the least stringent constraint is r{sub in}=5{sub −4}{sup +3} GM/c{sup 2}. Only two of the models deliver meaningful spin constraints, but a = 0.8 ± 0.2 is consistent with all of the fits. Overall, the data provide especially compelling evidence of an association between compact hard X-ray coronae and the base of relativistic radio jets in black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗

Bemnifosbuvir: An HCV NS5B Inhibitor With Multiple Modes of Action

Bemnifosbuvir (BEM) is a potent, pan-genotypic inhibitor targeting the hepatitis C virus (HCV) NS5B polymerase. Its antiviral activity was evaluated in an ascending dose phase I clinical trial involving 30 patients treated once a day for 7 days. After treatment initiation, plasma HCV RNA declined in a biphasic manner with a mean reduction of 2.3 log 10 IU/mL after 24 hours and 4.4 log10 IU/mL by day 7 for the highest dose. Alanine aminotransferase (ALT) also normalized in most patients. We employed a multiscale mathematical model fitted to the HCV RNA and ALT dynamics to quantify the antiviral activity and evaluate the modes of action of BEM. We found that models in which BEM only acted as a typical HCV RNA polymerase inhibitor and reduced the intracellular production of HCV RNA did not fit the data as well as models in which BEM had multiple modes of action, including suppressing viral assembly and secretion and enhancing intracellular HCV RNA degradation. BEM's effectiveness in inhibiting intracellular HCV RNA production increased with dose (150 mg/day: 88.2%, 300 mg/day: 98.8%, 600 mg/day: 99.5%), while inhibition of viral assembly and release was ~95% effective regardless of dose. We observed a dose-dependent enhancement in the degradation of intracellular HCV RNA, with degradation rates 1.5-fold higher in patients receiving 300 mg/day and 2.7-fold higher in those receiving 600 mg/day than in patients receiving 150 mg/day. No significant differences in antiviral activity were detected between HCV genotypes 1b and 3 or between patients with and without compensated cirrhosis.

59 BASIC BIOLOGICAL SCIENCES↗

Wideband dynamic radio spectra of two ultra-cool dwarfs

A number of radio-loud ultra-cool dwarf (UCD) stars exhibit both continuous broadband and highly polarized pulsed radio emission. In order to determine the nature of the emission and the physical characteristics in the source region, we have made multi-epoch, wideband spectral observations of TVLM 0513-46 and 2M 0746+20. We combine these observations with archival radio data to fully characterize both the temporal and spectral properties of the radio emission. The continuum spectral energy distribution can be well modeled using gyrosynchrotron emission from mildly relativistic electrons in a dipolar field. The pulsed emission exhibits a variety of time-variable characteristics, including frequency drifts, frequency cutoffs, and multiple pulses per period. For 2M 0746+20 we determine a pulse period consistent with previously determined values. We modeled locations of pulsed emission using an oblique rotating magnetospheric model with beamed electron-cyclotron maser (ECM) sources. The best-fit models have narrow ECM beaming angles aligned with the local source magnetic field direction, except for one isolated burst from 2M 0746+20. For TVLM 0513-46, the best-fit rotation axis inclination is nearly orthogonal to the line of sight. For 2M 0746+20 we found a good fit using a fixed inclination i=36{sup ∘}, determined from optical observations. For both stars the ECM sources are located near feet of magnetic loops with radial extents 1.2R{sub s}–2.7R{sub s} and surface fields 2.2–2.5 kG. These results support recent suggestions that radio over-luminous UCDs have a global “weak field” non-axisymmetric magnetic topologies.

79 ASTRONOMY AND ASTROPHYSICS↗

sparse_bias

This is a python package used to fit an unknown function from data that potential contains systematic biases related to metadata. The model fits the unknown function and uses a Bayesian horseshoe prior model to impose sparsity on the bias terms. This code has been generalized from research code developed for AIACHNE into a package that should have more general application in a wider class of statistical models.

Walton, Noah↗

Predicting environmental mitigation requirements for hydropower projects through the integration of biophysical and socio-political geographies

Uncertainty about environmental mitigation needs at existing and proposed hydropower projects makes it difficult for stakeholders to minimize environmental impacts. Hydropower developers and operators desire tools to better anticipate mitigation requirements, while natural resource managers and regulators need tools to evaluate different mitigation scenarios and order effective mitigation. Here we sought to examine the feasibility of using a suite of multidisciplinary explanatory variables within a spatially explicit modeling framework to fit predictive models for future environmental mitigation requirements at hydropower projects across the conterminous U.S. Using a database comprised of mitigation requirements from more than 300 hydropower project licenses, we were able to successfully fit models for nearly 50 types of environmental mitigation and to apply the predictive models to a set of more than 500 non-powered dams identified as having hydropower potential. The results demonstrate that mitigation requirements have been a result of a range of factors, from biological and hydrological to political and cultural. Furthermore, project developers can use these models to inform cost projections and design considerations, while regulators can use the models to more quickly identify likely environmental issues and potential solutions, hopefully resulting in more timely and more effective decisions on environmental mitigation.

13 HYDRO ENERGY↗

Radiation effects on atherosclerosis in atomic bomb survivors: a cross-sectional study using structural equation modeling

Past reports indicated that total-body irradiation at low to moderate doses could be responsible for cardiovascular disease risks, but the mechanism remains unclear. The purpose of this study was to investigate the association between radiation exposure and atherosclerosis, an underlying pathology of cardiovascular diseases, in the Japanese atomic bomb survivors. We performed a cross-sectional study measuring 14 clinical-physiological atherosclerosis indicators during clinical exams from 2010 to 2014 in 3274 participants of the Adult Health Study cohort. Multivariable analyses were performed by using a structural equation model with latent factors representing underlying atherosclerotic pathologies: (1) arterial stiffness, (2) calcification, and (3) plaque as measured with indicators chosen a priori on the basis of clinical-physiological knowledge. Radiation was linearly associated with calcification (standardized coefficient per Gy 0.15, 95 % confidence interval: CI [0.070, 0.23]) and plaque (0.11, 95 % CI [0.029, 0.20]), small associations that were comparable to about 2 years of aging per Gy of radiation exposure, but not with arterial stiffness (0.036, 95 % CI [-0.025, 0.095]). The model fitted better and had narrower confidence intervals than separate ordinary regression models explaining individual indicators independently. The associations were less evident when the dose range was restricted to a maximum of 2 or 1 Gy. By combining individual clinical-physiological indicators that are correlated because of common, underlying atherosclerotic pathologies, we found a small, but significant association of radiation with atherosclerosis.

59 BASIC BIOLOGICAL SCIENCES↗

qFeature

This package contains statistical routines for extracting features from multivariate time-series data which can then be used for subsequent multivariate statistical analysis to identify patterns and anomalous behavior. It calculates local linear or quadratic regression model fits to moving windows for each series and then summarizes the model coefficients across user-defined time intervals for each series. These methods are domain agnostic-but they have been successfully applied to a variety of domains, including commercial aviation and electric power grid data.

Source record↗

THE AGE OF THE KIC 7177553 SYSTEM

KIC 7177553 is a quadruple system containing two binaries of orbital periods 16.5 and 18 days. All components have comparable masses and are slowly rotating with spectral types of ∼G2V. The longer period binary is eclipsing with component masses and radii M {sub 1} = 1.043 ± 0.014 M {sub ⊙}, R {sub 1} = 0.940 ± 0.005 R {sub ⊙} and M {sub 2} = 0.986 ± 0.015 M {sub ⊙}, R {sub 2} = 0.941 ± 0.005 R {sub ⊙}. The essentially equal radii measurements are inconsistent with the two stars being on the man sequence at the same age using standard nonmagnetic stellar evolution models. Instead a consistent scenario is found if the stars are in their pre-main-sequence phase of evolution and have an age of 32–36 Myr. We have also computed evolutionary models of magnetic stars, but we find that our nonmagnetic models fit the empirical radii and effective temperatures better than the magnetic models.

79 ASTRONOMY AND ASTROPHYSICS↗

qFeature

This package contains statistical routines for extracting features from multivariate time-series data which can then be used for subsequent multivariate statistical analysis to identify patterns and anomalous behavior. It calculates local linear or quadratic regression model fits to moving windows for each series and then summarizes the model coefficients across user-defined time intervals for each series. These methods are domain agnostic-but they have been successfully applied to a variety of domains, including commercial aviation and electric power grid data.

Source record↗