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

ALMA Detection and Astrobiological Potential of Vinyl Cyanide on Titan

Recent simulations have indicated that vinyl cyanide is the best candidate molecule for the formation of cell membranes/vesicle structures in Titan's hydrocarbon-rich lakes and seas. Although the existence of vinyl cyanide (C2H3CN) on Titan was previously inferred using Cassini mass spectrometry, a definitive detection has been lacking until now. We report the first spectroscopic detection of vinyl cyanide in Titan's atmosphere, obtained using archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), collected from February to May 2014. We detect the three strongest rotational lines of C2H3CN in the frequency range of 230 to 232 GHz, each with greater than 4 sigma confidence. Radiative transfer modeling suggests that most of the C2H3CN emission originates at altitudes of approx. greater than 200 km, in agreement with recent photochemical models. The vertical column densities implied by our best-fitting models lie in the range of 3.7 x 10(exp 13) to 1.4 x 10(exp 14) cm(exp −2). The corresponding production rate of vinyl cyanide and its saturation mole fraction imply the availability of sufficient dissolved material to form approx. 10(exp 7) cell membranes/cu cm in Titan's sea Ligeia Mare.

Palmer, Maureen Y.↗

Uniqueness of a solution of a steady state photochemical problem: Applications to Mars

Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Krasnopolsky, Vladimir A.↗

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↗

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↗

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↗

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↗

Aeroelastic modeling for the FIT team F/A-18 simulation

Some details of the aeroelastic modeling of the F/A-18 aircraft done for the Functional Integration Technology (FIT) team's research in integrated dynamics modeling and how these are combined with the FIT team's integrated dynamics model are described. Also described are mean axis corrections to elastic modes, the addition of nonlinear inertial coupling terms into the equations of motion, and the calculation of internal loads time histories using the integrated dynamics model in a batch simulation program. A video tape made of a loads time history animation was included as a part of the oral presentation. Also discussed is work done in one of the areas of unsteady aerodynamic modeling identified as needing improvement, specifically, in correction factor methodologies for improving the accuracy of stability derivatives calculated with a doublet lattice code.

Zeiler, Thomas A.↗

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

We study systematics associated with estimating simple stellar population (SSP) parameters – age, metallicity [M/H], α-enhancement [α/Fe], and initial mass function (IMF) shape – and associated M * /L gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs), and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H β , 〈Fe〉, TiO 2SDSS , and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius R e : (a) ages are younger and [α/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about R e /2; (c) this makes M * /L about 2 × larger in the central regions than beyond R e /2. While the models of Conroy et al. return similar [M/H] and [α/Fe] profiles, the age and (hence) M * /L profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs, and S0s: younger ages and less enhanced [α/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred M * /L gradients significantly.

79 ASTRONOMY AND ASTROPHYSICS↗

Suzaku Observations of Thermal and Non-Thermal X-Ray Emission from the Middle-Aged Supernova Remnant G156.2+5.7

We present results from X-ray analysis of a Galactic middle-aged supernova remnant (SNR) G156.2+5.7 which is bright and largely extended in X-ray wavelengths, showing a clear circular shape (radius approx.50'). Using the Suzaku satellite, we observed this SNR in three pointings; partially covering the northwestern (NW) rim, the eastern (E) rim, and the central portion of this SNR. In the NW rim and the central portion, we confirm that the X-ray spectra consist of soft and hard-tail emission, while in the E rim we find no significant hard-tail emission. The soft emission is well fitted by either a one-component or two-component non-equilibrium ionization (NEI) model. In the NW and E rims, a one-component (the swept-up interstellar medium) NEI model well represents the soft emission. On the other hand, in the central portion, a two-component (the interstellar medium and the metal-rich ejecta) NEI model fits the soft emission better than the one-component NEI model from a statistical point of view. The relative abundances in the ejecta component suggest that G156.2+5.7 is a remnant from a core-collapse SN explosion whose progenitor mass is less than 15 Solar Mass. The origin of the hard-tail emission detected in the NW rim and the central portion of the SNR is highly likely non-thermal synchrotron emission from relativistic electrons. In the NW rim, the relativistic electrons seems to be accelerated by a forward shock with a slow velocity of APPROX.500 km/sec.

Katsuda, Satoru↗

Frozen-Flux Modelling for Epochs 1915 and 1980

The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.

OBrien, Michael S.↗

Relationship between ductility and the porosity of additively manufactured AlSi10Mg

Additive manufacturing via selective laser melting can result in variable levels of internal porosity both between build plates and within components from the same build. In this work, sample porosity levels were compared to tensile properties for 176 samples spanning eight different build plates. Sample porosity was measured both by Archimedes density, which provided an estimation of overall porosity, and by observation of voids in the fracture surface, which provided an estimation of the porosity at the failure plane. The porosity observed at the fracture surface consistently demonstrated higher porosity than that suggested by Archimedes density. The porosity values obtained from both methods were compared against the mechanical results. Sample porosity appears to have some correlation to the ultimate tensile strength, yield strength, and modulus, but the strongest relationship is observed between porosity and ductility. Three different models were used to relate the fracture surface porosity to the ductility. The first method was a simple linear regression analysis, while the other two models have been used to relate porosity to ductility in cast alloys. It is shown that all three models fit the data well over the observed porosity ranges, suggesting that the models taken from casting theory can extend to additively manufactured metals. Finally, it is proposed that the non-destructive Archimedes method could be used to estimate an approximate sample ductility through the use of correlations realized here. Such a relationship could prove useful for design and for a deeper understanding of the impact of pores on tensile behavior.

36 MATERIALS SCIENCE↗

Investigation of low-latitude hydrogen emission in terms of a two-component interstellar gas model

High-resolution 21-cm hydrogen line observations at low galactic latitude are analyzed to determine the large-scale distribution of galactic hydrogen. Distribution parameters are found by model fitting, optical depth effects are computed using a two-component gas model suggested by the observations, and calculations are made for a one-component uniform spin-temperature gas model to show the systematic departures between this model and data obtained by incorrect treatment of the optical depth effects. Synthetic 21-cm line profiles are computed from the two-component model, and the large-scale trends of the observed emission profiles are reproduced together with the magnitude of the small-scale emission irregularities. Values are determined for the thickness of the galactic hydrogen disk between half density points, the total observed neutral hydrogen mass of the galaxy, and the central number density of the intercloud hydrogen atoms. It is shown that typical hydrogen clouds must be between 1 and 13 pc in diameter and that optical thinness exists on large-scale despite the presence of optically thin gas.

Baker, P. L.↗

Calculating Robot-Joint Coordinates From Image Coordinates

Detailed knowledge of robot joints not required. Algorithm generates approximate mathematical models of coordinates of joints of robot as functions of coordinates of points in images of work region viewed by television cameras. Joint coordinates necessary to position and orient end effector calculated by mathematical models fitted to experimentally determined data on positions, orientations, and joint coordinates. Generates models as functions of desired location of end effector of robot. Does not require priori knowledge of kinematic equations of robot.

Source record↗

Software reliability studies

There are many software reliability models which try to predict future performance of software based on data generated by the debugging process. Our research has shown that by improving the quality of the data one can greatly improve the predictions. We are working on methodologies which control some of the randomness inherent in the standard data generation processes in order to improve the accuracy of predictions. Our contribution is twofold in that we describe an experimental methodology using a data structure called the debugging graph and apply this methodology to assess the robustness of existing models. The debugging graph is used to analyze the effects of various fault recovery orders on the predictive accuracy of several well-known software reliability algorithms. We found that, along a particular debugging path in the graph, the predictive performance of different models can vary greatly. Similarly, just because a model 'fits' a given path's data well does not guarantee that the model would perform well on a different path. Further we observed bug interactions and noted their potential effects on the predictive process. We saw that not only do different faults fail at different rates, but that those rates can be affected by the particular debugging stage at which the rates are evaluated. Based on our experiment, we conjecture that the accuracy of a reliability prediction is affected by the fault recovery order as well as by fault interaction.

Hoppa, Mary Ann↗

Downscattering due to Wind Outflows in Compact X-ray Sources: Theory and Interpretation

A number of recent lines of evidence point towards the presence of hot, outflowing plasma from the central regions of compact Galactic and extragalactic X-ray sources. Additionally, it has long been noted that many of these sources exhibit an "excess" continuum component, above approx. 10 keV, usually attributed to Compton Reflection from a static medium. Motivated by these facts, as well as by recent observational constraints on the Compton reflection models - specifically apparently discrepant variability timescales for line and continuum components in some cases - we consider possible of effects of out-flowing plasma on the high-energy continuum spectra of accretion powered compact objects. We present a general formulation for photon downscattering diffusion which includes recoil and Comptonization effects due to divergence of the flow. We then develop an analytical theory for the spectral formation in such systems that allows us to derive formulae for the emergent spectrum. Finally we perform the analytical model fitting on several Galactic X-ray binaries. Objects which have been modeled with high-covering-fraction Compton reflectors, such as GS1353-64 are included in our analysis. In addition, Cyg X-3, is which is widely believed to be characterized by dense circumstellar winds with temperature of order 10(exp 6) K, provides an interesting test case. Data from INTEGRAL and RXTE covering the approx. 3 - 300 keV range are used in our analysis. We further consider the possibility that the widely noted distortion of the power-law continuum above 10 keV may in some cases be explained by these spectral softening effects.

Titarchuk, Lev↗

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗

Detection of Pulsed X-Ray Emission from XMM-Newton Observations of PSR J0538+2817

We report on the XMM-Newton observations of the 143 ms pulsar PSR J0538+2817. We present evidence for the first detections of pulsed X-rays from the source at a frequency that is consistent with the predicted radio frequency. The pulse profile is broad and asymmetric, with a pulse fraction of 18% +/- 3%. We find that the spectrum of the source is well-fitted with a blackbody with T(sup infinity) = (2.12 (sup +0.04) (sub -0.03)) x 10(exp 6) K and N(sub H) = 2.5 x 10(exp 2)/sq cm. The radius determined from the model fit of 1.68 +/- 0.05 km suggests that the emission is from a heated polar cap. A fit to the spectra with an atmospheric model reduces the inferred temperature and, hence, increases the radius of the emitting region; however, the pulsar distance determined from the fit is then smaller than the dispersion distance.

McGowan, K. E.↗