Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

Measurements of the Hubble Constant with a Two-rung Distance Ladder: Two Out of Three Ain’t Bad

Abstract The three-rung distance ladder, which calibrates Type Ia supernovae (SNe Ia) through stellar distances linked to geometric measurements, provides the highest precision direct measurement of the Hubble constant. In light of the Hubble tension, it is important to test the individual components of the distance ladder. For this purpose, we report a measurement of the Hubble constant from 35 extragalactic Cepheid hosts measured by the SH0ES team, using their distances and redshifts at cz ≤ 3300 km s −1 , instead of any more distant SNe Ia, to measure the Hubble flow. The Cepheid distances are calibrated geometrically in the Milky Way, NGC 4258, and the Large Magellanic Cloud. Peculiar velocities are a significant source of systematic uncertainty at z ∼ 0.01, and we present a formalism for both mitigating and quantifying their effects, making use of external reconstructions of the density and velocity fields in the nearby universe. We identify a significant source of uncertainty originating from different assumptions about the selection criteria of this sample, whether distance or redshift limited, as it was assembled over three decades. Modeling these assumptions yields central values ranging from H 0 = 71.7 to 76.4 km s −1 Mpc −1 . Combining the four best-fitting selection models yields H 0 = 72.9 − 2.2 + 2.4 as a fiducial result, at 2.4 σ tension with Planck. While SNe Ia are essential for a precise measurement of H 0 , unknown systematics in these supernovae are unlikely to be the source of the Hubble tension.

79 ASTRONOMY AND ASTROPHYSICS↗

Anisotropies of the cosmic microwave background in nonstandard cold dark matter models

Small angular scale cosmic microwave anisotropies in flat, vacuum-dominated, cold dark matter cosmological models which fit large-scale structure observations and are consistent with a high value for the Hubble constant are reexamined. New predictions for CDM models in which the large-scale power is boosted via a high baryon content and low H(0) are presented. Both classes of models are consistent with current limits: an improvement in sensitivity by a factor of about 3 for experiments which probe angular scales between 7 arcmin and 1 deg is required, in the absence of very early reionization, to test boosted CDM models for large-scale structure formation.

Vittorio, Nicola↗

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗

Bulge Growth and Quenching Since Z=2.5 in Candels/3D-HST

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift informationprovided by 3D-HST, we investigate the relation between structure and stellar populations fora mass-selected sample of 6764 galaxies above 1010 M, spanning the redshift range 0.5 z 2.5.For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e.,bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar massmaps reconstructed from resolved stellar population modeling. We confirm that the increased bulgeprominence among quiescent galaxies, as reported previously based on rest-optical observations, remainsin place when considering the distributions of stellar mass. Moreover, we observe an increaseof the typical Sersic index and bulge-to-total ratio (with median BT reaching 40-50) among starforminggalaxies above 1011 M. Given that quenching for these most massive systems is likely tobe imminent, our findings suggest that significant bulge growth precedes a departure from the starformingmain sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge andtotal mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy thanthe total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model bySomerville et al. In the latter, bulges and black holes grow hand in hand through merging andordisk instabilities, and AGN-feedback shuts off star formation. Further observations will be requiredto pin down star formation quenching mechanisms, but our results imply they must be internal to thegalaxies and closely associated with bulge growth.

Quenching↗

High-Resolution 2D Lg and Pg Attenuation Models in the Basin and Range Region with Implications for Frequency-Dependent Q

Here, we develop high–resolution, laterally varying attenuation models for the regional crustal phases of Pg and Lg in the area surrounding the Basin and Range Province in the western United States. The models are part of the characterization effort for the Source Physics Experiment (SPE), a series of chemical explosions at the Nevada National Security Site designed to improve our understanding of explosion source phenomenology. To aid in SPE modeling efforts, we focus on improving our ability to accurately predict amplitudes in a set of narrow frequency bands ranging from 0.5 to 16.0 Hz. To explore constraints at higher frequencies where data become more sparse, we test the robustness of the empirically observed power–law relationship between quality factor Q and frequency (Q=Q 0 f γ ). Our methodology uses a staged approach to consider attenuation, physics–based source terms, site terms, and geometrical spreading contributions to amplitude measurements. Tomographic inversion results indicate that the frequency dependence is a reasonable assumption as attenuation varies laterally for this region through all frequency bands considered. Our 2D Pg and Lg attenuation models correlate with underlying physiographic provinces, with the highest Q located in the Sierra Nevada Mountains and the Colorado plateau. Compared to a best–fitting 1D model for the region, the 2D model provides an 81% variance reduction overall for Lg residuals and a 75% reduction for Pg . These detailed attenuation maps at high frequencies will facilitate further study of local and regional distance P/S amplitude discriminants that are typically used to distinguish between earthquakes and underground explosions.

58 GEOSCIENCES↗

Herschel and Hubble Study of a Lensed Massive Dusty Starbursting Galaxy at z ∼ 3

We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared, and Herschel far-infrared observations of an extremely star-forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE-identified galaxy cluster at z ∼ 1. The lensed galaxy is spectroscopically confirmed at z = 2.685 from detection of CO (1→0) by GBT and from detection of CO (3→2) obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed model of the background dusty lensed submillimeter galaxy (SMG), which allows us to study the source plane properties of the target. The best-fit lens model provides magnifications of μ {sub star} = 2.10 ± 0.11 and μ {sub dust} = 2.02 ± 0.06 for the stellar and dust components, respectively. Multiband data yield a magnification-corrected star formation rate of 1900(±200) M {sub ⊙} yr{sup −1} and a stellar mass of 6.8{sub −2.7}{sup +0.9}×10{sup 11} M{sub ⊙}, consistent with a main sequence of star formation at z ∼ 2.6. The CO observations yield a molecular gas mass of 8.3(±1.0) × 10{sup 10} M {sub ⊙}, similar to the most massive star-forming galaxies, which together with the high star formation efficiency, are responsible for the intense observed star formation rates. The lensed SMG has a very short gas depletion timescale of ∼40 Myr. The high stellar mass and small gas fractions observed indicate that the lensed SMG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling battery life through changes in voltage fit coefficients

A number of 12 ampere-hour cells, nickel cadmium, were cycled under various conditions of temperature and depth of discharge. Using this data, it was confirmed that a five parameter fit equation could be used to model the data within a few millivolts. Both charge and discharge curves can be fit with accuracies in the range of one or two millivolts. The fit coefficients when plotted versus cycles show definite trends and patterns which can be used in an operational sense to predict battery voltage as a function of temperature and depth of discharge.

Fox, D.↗

Report on twisted nematic and supertwisted nematic device characterization program

In this study we measured the optical characteristics of normally white twisted nematic (NWTN) and super twisted nematic (STN ) cells. Though no dynamic computer model was available, the static observations were compared with computer simulated behavior. The measurements were taken as a function of both viewing angle and applied voltage and included in the static case not only luminance but also contrast ratio and chromaticity . We employed the computer model Twist Cell Optics, developed at Kent State in conjunction with this study, and whose optical modeling foundation, Iike the ViDEOS program, is the 4 x 4 matrix method of Berreman. In order to resolve discrepancies between the experimental and modeled data the optical parameters of the individual cell components, where not known, were determined using refractometry, profilometry, and various forms of ellipsometry. The resulting agreement between experiment and model is quite good due primarily to a better understanding of the structure and optics of dichroic sheet polarizers. A description of the model and test cells employed are given in section 2. Section 3 contains the experimental data gathered and section 4 gives examples of the fit between model and experiment. Also included with this report are a pair of papers which resulted from the research and which detail the polarizer properties and some of the cell characterization methods.

Source record↗

Diffuse Galactic Continuum Gamma Rays. A Model Compatible with EGRET Data and Cosmic-ray Measurements

We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma-rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission it- self. The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly de- tails of Galactic structure. This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the y-ray spectrum better and to the highest EGRET energies. It gives a good representation of the latitude distribution of the y-ray emission from the plane to the poles, and of the longitude distribution. We show that secondary positrons and electrons make an essential contribution to Galactic diffuse y-ray emission.

Strong, Andrew W.↗

Experimental Determination of the Dynamic Hydraulic Transfer Function for the J-2X Oxidizer Turbopump: Methodology - Part One

An advanced methodology for extracting the hydraulic dynamic pump transfer matrix (Yp) for a cavitating liquid rocket engine turbopump inducer+impeller has been developed. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Laboratory pulsed subscale waterflow test of the J-2X oxygen turbo pump is introduced and our new extraction method applied to the data collected. From accurate measures of pump inlet and discharge perturbational mass flows and pressures, and one-dimensional flow models that represents complete waterflow loop physics, we are able to derive Yp and hence extract the characteristic pump parameters: compliance, pump gain, impedance, mass flow gain. Detailed modeling is necessary to accurately translate instrument plane measurements to the pump inlet and discharge and extract Yp. We present the MSFC Dynamic Lump Parameter Fluid Model Framework and describe critical dynamic component details. We report on fit minimization techniques, cost (fitness) function derivation, and resulting model fits to our experimental data are presented. Comparisons are made to alternate techniques for spatially translating measurement stations to actual pump inlet and discharge.

Zoladz, Tom↗

Boron Coordination in Multicomponent Glasses: Analytical Models and Machine Learning With Uncertainty

Borosilicate glasses are extensively used in a variety of applications from kitchenware to nuclear waste immobilization due to the strong network formed by the Si-O-B bond that makes it resistant to chemical corrosion and gives it a low thermal expansion. Boron, however, exists in both trigonal BO3 and tetrahedral BO4 bonds in glass systems, which impacts the chemical durability and thermal resistance of the glass, amongst other properties. Boron coordination (N4), or the ratio of the amount of BO4 to BO3 within a glass, may aid in predicting these properties but is difficult to derive without experimental data due to the complexity of impacts from varied glass compositions and processing factors. For this reason, compositional models have been developed to predict boron coordination, but the models typically include a limited number of glass components. To help fill this gap in the models, in this work, a diverse multicomponent glass dataset of 809 glasses is compiled from a literature search, and then a number of analytical and machine learning (ML) models are trained on the dataset. Previously developed modified Bernstein and modified Du Stebbins analytical models were fitted to update parameters with the new dataset. Then, partially Bayesian neural networks, Gaussian process regressor, and heteroskedastic deterministic neural networks were evaluated. The ML models examined all have different strategies to overcome the potential for overfitting as a result of a limited training dataset, and return results that account for model uncertainty, which can be valuable for understanding model reliability. For the first time, cooling rate is introduced as an input parameter for ML models, showing consistent improvements in performance and solidifying the importance of including parameters outside of composition alone for N4 prediction. The machine learning models examined here show promise in accurate predictions of boron coordination in borosilicate glasses, all achieving R2 values of 0.91.

boron coordination↗

MAUD Interface Tool Kit (MILK)

Materials Analysis Using Diffraction (MAUD) is an open source Rietveld refinement program which fits diffraction models to diffraction spectra allowing quantitative diffraction analysis. The scripting interface developed here using python facilitates the choice of parameters to refine during the Rietveld fitting process and provides a simple MPI framework for running MAUD instances in parallel. The scripting language facilitates batch, custom, and reproducible Rietveld refinements of large datasets.

Savage, Daniel↗

The use of a helicopter mounted ranging scatterometer for estimation of extinction and scattering properties of forest canopies-II: Experimental results for high-density aspen

An analytic expression is derived that describes the backscatter power from a semi-infinite plane parallel homogeneous canopy as a function of distance from an airborne radar. This model is fitted to observed data for a high-density aspen canopy by a modification of a technique developed by Tyapkin (1960). This inversion of the model provides unbiased estimates of the canopy extinction and backscattering parameters. An active radar calibrator located underneath the canopy provides an independent method of determining the volume extinction coefficient. The results reported indicate that the coefficients change throughout the year. A comparison of these coefficients with Eom and Funk's (1984) disk model, using measured canopy properties, shows that at C-band frequency, only a part of the scattering and absorption can be attributed to the canopy leaves.

Pitts, David E.↗

Satellite Sensed Skin Sea Surface Temperature

Quantitative predictions of spatial and temporal changes the global climate rely heavily on the use of computer models. Unfortunately, such models cannot provide the basis for climate prediction because key physical processes are inadequately treated. Consequently, fine tuning procedures are often used to optimize the fit between model output and observational data and the validation of climate models using observations is essential if model based predictions of climate change are to be treated with any degree of confidence. Satellite Sea Surface Temperature (SST) observations provide high spatial and temporal resolution data which is extremely well suited to the initialization, definition of boundary conditions and, validation of climate models. In the case of coupled ocean-atmosphere models, the SST (or more correctly the 'Skin' SST (SSST)) is a fundamental diagnostic variable to consider in the validation process. Daily global SST maps derived from satellite sensors also provide adequate data for the detection of global patterns of change which, unlike any other SST data set, repeatedly extend into the southern hemisphere extra-tropical regions. Such data are essential to the success of the spatial 'fingerprint' technique, which seeks to establish a north-south asymmetry where warming is suppressed in the high latitude Southern Ocean. Some estimates suggest that there is a greater than 80% chance of directly detecting significant change (97.5 % confidence level) after 10-12 years of consistent global observations of mean sea surface temperature. However, these latter statements should be qualified with the assumption that a negligible drift in the observing system exists and that biases between individual instruments required to derive a long term data set are small. Given that current estimates for the magnitude of global warming of 0.015 K yr(sup -1) - 0.025 K yr(sup -1), satellite SST data sets need to be both accurate and stable if such a warming trend is to be confidently detected. Some of these activities are focussed to develop and deploy instrumentation suitable for the collection of precise in situ measurements of the SSST which can be used to improve the accuracy of satellite measurements, while others develop techniques to generate improved global analyses of sea surface temperature using historical data.

Donlon, Craig↗

AN EXTERNAL SHOCK ORIGIN OF GRB 141028A

The prompt emission of the long, smooth, and single-pulsed gamma-ray burst, GRB 141028A, is analyzed under the guise of an external shock model. First, we fit the γ -ray spectrum with a two-component photon model, namely, synchrotron+blackbody, and then fit the recovered evolution of the synchrotron νF{sub ν} peak to an analytic model derived considering the emission of a relativistic blast wave expanding into an external medium. The prediction of the model for the νF{sub ν} peak evolution matches well with the observations. We observe the blast wave transitioning into the deceleration phase. Furthermore, we assume the expansion of the blast wave to be nearly adiabatic, motivated by the low magnetic field deduced from the observations. This allows us to recover within an order of magnitude the flux density at the νF{sub ν} peak, which is remarkable considering the simplicity of the analytic model. Under this scenario we argue that the distinction between prompt and afterglow emission is superfluous as both early-time emission and late-time emission emanate from the same source. While the external shock model is clearly not a universal solution, this analysis opens the possibility that at least some fraction of GRBs can be explained with an external shock origin of their prompt phase.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics‐based iterative reconstruction for dual‐source and flying focal spot computed tomography

Purpose For single‐source helical Computed Tomography (CT), both Filtered‐Back Projection (FBP) and statistical iterative reconstruction have been investigated. However, for dual‐source CT with flying focal spot (DS‐FFS CT), a statistical iterative reconstruction that accurately models the scanner geometry and acquisition physics remains unknown to researchers. Therefore, our purpose is to present a novel physics‐based iterative reconstruction method for DS‐FFS CT and assess its image quality. Methods Our algorithm uses precise physics models to reconstruct from the native cone‐beam geometry and interleaved dual‐source helical trajectory of a DS‐FFS CT. To do so, we construct a noise physics model to represent data acquisition noise and a prior image model to represent image noise and texture. In addition, we design forward system models to compute the locations of deflected focal spots, the dimension, and sensitivity of voxels and detector units, as well as the length of intersection between x‐rays and voxels. The forward system models further represent the coordinated movement between the dual sources by computing their x‐ray coverage gaps and overlaps at an arbitrary helical pitch. With the above models, we reconstruct images by an advanced Consensus Equilibrium (CE) numerical method to compute the maximum a posteriori estimate to a joint optimization problem that simultaneously fits all models. Results We compared our reconstruction with Siemens ADMIRE, which is the clinical standard hybrid iterative reconstruction (IR) method for DS‐FFS CT, in terms of spatial resolution, noise profile, and image artifacts through both phantoms and clinical scan datasets. Experiments show that our reconstruction has a higher spatial resolution, with a Task‐Based Modulation Transfer Function (MTF task ) consistently higher than the clinical standard hybrid IR. In addition, our reconstruction shows a reduced magnitude of image undersampling artifacts than the clinical standard. Conclusions By modeling a precise geometry and avoiding data rebinning or interpolation, our physics‐based reconstruction achieves a higher spatial resolution and fewer image artifacts with smaller magnitude than the clinical standard hybrid IR.

Wang, Xiao↗

Modeling inter-particle magnetic correlations in magnetite nanoparticle assemblies using x-ray magnetic scattering data

Magnetic nanoparticles are increasingly used in nanotechnologies and biomedical applications, such as drug targeting, MRI, bio-separation. Magnetite (Fe 3 O 4 ) nanoparticles stand to be effective in these roles due to the non-toxic nature of magnetite and its ease of manufacture. To be more effective in these applications, a greater understanding of the magnetic behavior of a collection of magnetite nanoparticles is needed. This research seeks to discover the local magnetic ordering of ensembles of magnetite nanoparticles occurring under various external fields. To complete this study, we use x-ray resonant magnetic scattering (XRMS). Here we discuss the modeling of the magnetic scattering data using a one-dimensional chain of nanoparticles with a mix of ferromagnetic, anti-ferromagnetic, and random orders. By fitting the model to the experimental data, we extracted information about the magnetic correlations in the nanoparticle assembly.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Active galactic nuclei, neutrinos, and interacting cosmic rays in NGC 253 and NGC 1068

The galaxies M82, NGC 253, NGC 1068, and NGC 4945 have been detected in γ-rays by Fermi. Previously, we developed and tested a model for cosmic-ray interactions in the starburst galaxy M82. Now, we aim to explore the differences between starburst and active galactic nucleus (AGN) environments by applying our self-consistent model to the starburst galaxy NGC 253 and the Seyfert galaxy NGC 1068. Assuming a constant cosmic-ray acceleration efficiency by supernova remnants with Milky Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations, predict the radio and γ-ray spectra, and compare with published measurements. We find that our models easily fit the observed γ-ray spectrum for NGC 253 while constraining the cosmic-ray source spectral index and acceleration efficiency. However, we encountered difficultly modeling the observed radio data and constraining the speed of the galactic wind and the magnetic field strength, unless the gas mass is less than currently preferred values. Additionally, our starburst model consistently underestimates the observed γ-ray flux and overestimates the radio flux for NGC 1068; these issues would be resolved if the AGN is the primary source of γ-rays. We discuss the implications of these results and make predictions for the neutrino fluxes for both galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗