Search NASA⌕ Search

SEARCH · Search NASA

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,513 records · Page 84

A unified exploration of the chronology of the Galaxy

The Milky Way has distinct structural stellar components linked to its formation and subsequent evolution, but disentangling them is non-trivial. With the recent availability of high-quality data for a large numbers of stars in the Milky Way, it is a natural next step for research in the evolution of the Galaxy to perform automated explorations with unsupervised methods of the structures hidden in the combination of large-scale spectroscopic, astrometric, and asteroseismic data sets. We determine precise stellar properties for 21 076 red giants, mainly spanning 2–15 kpc in Galactocentric radii, making it the largest sample of red giants with measured asteroseismic ages available to date. We explore the nature of different stellar structures in the Galactic disc by using Gaussian mixture models as an unsupervised clustering method to find substructure in the combined chemical, kinematic, and age subspace. The best-fitting mixture model yields four distinct physical Galactic components in the stellar disc: the thin disc, the kinematically heated thin disc, the thick disc, and the stellar halo. We find hints of an age asymmetry between the Northern and Southern hemisphere, and we measure the vertical and radial age gradient of the Galactic disc using the asteroseismic ages extended to further distances than previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗

ACCRETION FLOW DYNAMICS OF MAXI J1836-194 DURING ITS 2011 OUTBURST FROM TCAF SOLUTION

The Galactic transient X-ray binary MAXI J1836-194 was discovered on 2011 August 29. Here we make a detailed study of the spectral and timing properties of its 2011 outburst using archival data from the RXTE Proportional Counter Array instrument. The evolution of accretion flow dynamics of the source during the outburst through spectral analysis with Chakrabarti–Titarchuk’s two-component advective flow (TCAF) solution as a local table model in XSPEC. We also fitted spectra with combined disk blackbody and power-law models and compared it with the TCAF model fitted results. The source is found to be in hard and hard-intermediate spectral states only during the entire phase of this outburst. No soft or soft-intermediate spectral states are observed. This could be due to the fact that this object belongs to a special class of sources (e.g., MAXI J1659-152, Swift J1753.5-0127, etc.) that have very short orbital periods and that the companion is profusely mass-losing or the disk is immersed inside an excretion disk. In these cases, flows in the accretion disk are primarily dominated by low viscous sub-Keplerian flow and the Keplerian rate is not high enough to initiate softer states. Low-frequency quasi-periodic oscillations (QPOs) are observed sporadically although as in normal outbursts of transient black holes, monotonic evolutions of QPO frequency during both rising and declining phases are observed. From the TCAF fits, we find the mass of the black hole in the range of 7.5–11 M{sub ⊙}, and from time differences between peaks of the Keplerian and sub-Keplerian accretion rates we obtain a viscous timescale for this particular outburst, ∼10 days.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of the SSRCS engine with hydrazine as a fuel, phase 2

Test firing of the Space Shuttle Reaction Control Thruster (SSRCT) was conducted to determine the characteristic velocity and chamber throat temperatures when the fuel is changed from monomethylhydrazine (MMH) to hydrazine. Thruster performance with hydrazine was essentially as predicted. Characteristic velocity at a mixture ratio of 1.4 (equal volumetric flow) for hydrazine was 5180 feet/second compared with 5110 feet/second at a mixture of 1.6 for MMH. Specific impulse with 22:1 nozzle is calculated to be 280 lbf-sec/lbm. Thermal performance, as measured by chamber throat temperature, was dramatically different (colder) than predicted. Throat temperatures of 2330 F were predicted from a cooling model which assumed a reactive liquid cooling film. Throat temperatures of 800 - 1000 F were measured. An attempt to reconcile the differences between film cooling predictions and the measurements, indicated that an evaporation (nonreactive fluid) model rather than a decomposition (reactive fluid) model best fits the test data at 0/F = 1.4. At higher mixture ratios, a transition to temperatures typical of the reactive fluid model were observed.

Minton, S. J.↗

Bimodal star formation - Constraints from galaxy colors at high redshift

The possibility that at early epochs the light from elliptical galaxies is dominated by stars with an initial mass function (IMF) which is deficient in low-mass stars, relative to the solar neighborhood is investigated. V-R colors for the optical counterparts of 3CR radio sources offer the most severe constraints on the models. Reasonable fits are obtained to both the blue, high-redshift colors and the redder, low-redshift colors with a model galaxy which forms with initially equal star formation rates in each of two IMF modes: one lacking low-mass stars, and one with stars of all masses. The net effect is that the time-integrated IMF has twice as many high-mass stars as the solar neighborhood IMF, relative to low mass stars. A conventional solar neighborhood IMF does not simultaneously account for both the range in colors at high redshift and the redness of nearby ellipticals, with any single star formation epoch. Models with a standard IMF require half the stellar population to be formed in a burst at low redshift z of about 1.

Wyse, Rosemary F. G.↗

The Interstellar Polarization Feature Associated with the 2175 Angstrom Extinction Bump

The most successful model at fitting the wavelength dependence of interstellar extinction consists of two populations of bare silicate and graphite grains (Mathis, Rumpl & Nordsieck, (MRN)). The graphite grains are needed to fit the strong 2175 A extinction bump, with the silicates providing most of the smooth extinction seen at other wavelengths. From observations of the IR silicate absorption features, it was inferred that the silicate grains, non-spherical and aligned to the Galactic magnetic field, were responsible for the interstellar polarization. Aligned silicate grains make a very good fit to the smooth featureless UV polarization curves that comprise most of the observed sample. The lines of sight showing a UV polarization feature require something different, perhaps a second population of grains as do the fits to the extinction curve. The polarization wavelength dependence of HD 197770 can be fit by a MRN-like mixture of aligned bare silicate and graphite grains. However, none of the popular grain models including MRN envisaged an aligned bump grain population. It has been suggested that the extinction bump is not due to graphite but rather to a population of polycyclic aromatic hydrocarbons (PAHs). If the observed UV polarization features could be associated with the 2175 A bump then it would strengthen the argument that the bump must be due to grains rather than PAHs which are not likely to be aligned to the Galactic Magnetic field. The reality of ultraviolet polarization features has been confirmed by observations with ASTRO-2. The original detection toward HD 197770 has been confirmed and a new feature has been detected toward HD 147933-4. Both features have centroids lying close to 2175 A the location of the UV extinction bump. Two possibilities are considered for the source of the polarization bump, alignment of graphite grains responsible for the extinction bump or changes in the size distribution of the aligned silicate grains responsible for the continuum polarization. It seems the graphite grains are a more likely source of the polarization.

Clayton, Geoffrey C.↗

Neural network uncertainty assessment using Bayesian statistics: a remote sensing application

Neural network (NN) techniques have proved successful for many regression problems, in particular for remote sensing; however, uncertainty estimates are rarely provided. In this article, a Bayesian technique to evaluate uncertainties of the NN parameters (i.e., synaptic weights) is first presented. In contrast to more traditional approaches based on point estimation of the NN weights, we assess uncertainties on such estimates to monitor the robustness of the NN model. These theoretical developments are illustrated by applying them to the problem of retrieving surface skin temperature, microwave surface emissivities, and integrated water vapor content from a combined analysis of satellite microwave and infrared observations over land. The weight uncertainty estimates are then used to compute analytically the uncertainties in the network outputs (i.e., error bars and correlation structure of these errors). Such quantities are very important for evaluating any application of an NN model. The uncertainties on the NN Jacobians are then considered in the third part of this article. Used for regression fitting, NN models can be used effectively to represent highly nonlinear, multivariate functions. In this situation, most emphasis is put on estimating the output errors, but almost no attention has been given to errors associated with the internal structure of the regression model. The complex structure of dependency inside the NN is the essence of the model, and assessing its quality, coherency, and physical character makes all the difference between a blackbox model with small output errors and a reliable, robust, and physically coherent model. Such dependency structures are described to the first order by the NN Jacobians: they indicate the sensitivity of one output with respect to the inputs of the model for given input data. We use a Monte Carlo integration procedure to estimate the robustness of the NN Jacobians. A regularization strategy based on principal component analysis is proposed to suppress the multicollinearities in order to make these Jacobians robust and physically meaningful.

Neural Networks (Computer)↗

Nustar and Chandra Insight into the Nature of the 3-40 Kev Nuclear Emission in Ngc 253

We present results from three nearly simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR) and Chandra monitoring observations between 2012 September 2 and 2012 November 16 of the local star-forming galaxy NGC 253. The 3-40 kiloelectron volt intensity of the inner approximately 20 arcsec (approximately 400 parsec) nuclear region, as measured by NuSTAR, varied by a factor of approximately 2 across the three monitoring observations. The Chandra data reveal that the nuclear region contains three bright X-ray sources, including a luminous (L (sub 2-10 kiloelectron volt) approximately few × 10 (exp 39) erg per s) point source located approximately 1 arcsec from the dynamical center of the galaxy (within the sigma 3 positional uncertainty of the dynamical center); this source drives the overall variability of the nuclear region at energies greater than or approximately equal to 3 kiloelectron volts. We make use of the variability to measure the spectra of this single hard X-ray source when it was in bright states. The spectra are well described by an absorbed (power-law model spectral fit value, N(sub H), approximately equal to 1.6 x 10 (exp 23) per square centimeter) broken power-law model with spectral slopes and break energies that are typical of ultraluminous X-ray sources (ULXs), but not active galactic nuclei (AGNs). A previous Chandra observation in 2003 showed a hard X-ray point source of similar luminosity to the 2012 source that was also near the dynamical center (Phi is approximately equal to 0.4 arcsec); however, this source was offset from the 2012 source position by approximately 1 arcsec. We show that the probability of the 2003 and 2012 hard X-ray sources being unrelated is much greater than 99.99% based on the Chandra spatial localizations. Interestingly, the Chandra spectrum of the 2003 source (3-8 kiloelectron volts) is shallower in slope than that of the 2012 hard X-ray source. Its proximity to the dynamical center and harder Chandra spectrum indicate that the 2003 source is a better AGN candidate than any of the sources detected in our 2012 campaign; however, we were unable to rule out a ULX nature for this source. Future NuSTAR and Chandra monitoring would be well equipped to break the degeneracy between the AGN and ULX nature of the 2003 source, if again caught in a high state.

Chandra↗

Characteristics of DSD Bulk Parameters: Implication to Radar Rain Retrieval

With the use of 213,456 one-minute measured data of droplet-size distribution (DSD) of rain collected during several National Aeronautics and Space Administration (NASA)-sponsored field campaigns, the relationships between rainfall rate R, mass-weighted diameter D(m) and normalized intercept parameter N(w) of the gamma DSD are studied. It is found, based on the simulations of the gamma DSD model, that R, D(m) and N(w) are closely interrelated, and that the ratio of R to N(w) is solely a function of D(m), independent of the shape factor μ of the gamma distribution. Furthermore, the model-produced ratio agrees well with those from the DSD data. When a power-law equation is applied to fit the model data, we have: R=aN(w)D(sup b, sub m) , where a=1.588×10(exp −4) , b=4.706 . Analysis of two-parameter relationships such as R–D(m), N(w)–R and N(w)–D(m) reveals that R and D(m) are moderately correlated while N(w) and D(m) are negatively correlated. N(w) and R, however, are uncorrelated. The gamma DSD model also reveals that variation of R–D(m) relation is caused primarily by N(w). For the application of the Ku- and Ka-band dual-frequency radar for the retrieval of the DSD bulk parameters as well as the specific radar attenuations, the study is carried out to relate the dual-frequency radar reflectivity factors to the DSD and attenuation parameters.

Liang Liao↗

Soil Water Characteristics of Cores from Low- and High-Centered Polygons, Utqiagvik (Barrow), Alaska, 2012

This dataset includes soil water characteristic curves for soil and permafrost in two representative frozen cores collected from a high-center polygon (HCP) and a low-center polygon (LCP) from the Barrow Environmental Observatory. Data include soil water content and soil water potential measured using the simple evaporation method for hydrological and biogeochemical simulations and experimental data analysis. Data can be used to generate a soil moisture characteristic curve, which can be fit to a variety of hydrological functions to infer critical parameters for soil physics. Considering the measured the soil water properties, the van Genuchten model predicted well the HCP, in contrast, the Kosugi model well fitted LCP which had more saturated condition. This dataset includes 8 zipped files, one txt files, and one pdf file.Data were collected to estimate soil water potential and water availability at various temperatures and conditions for hydrological and biogeochemical simulations and experimental data analysis. Data can be used to generate a soil moisture characteristic curve, which can be fit to a variety of hydrological functions to infer critical parameters for soil physics.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

The Syrtis Major volcano, Mars: A multidisciplinary approach to interpreting its magmatic evolution and structural development: EVOLUTION OF THE SYRTIS MAJOR VOLCANO

Very weak crustal magnetic fields over the Syrtis Major volcanic complex imply almost total thermal demagnetization via magmatic intrusions over a large area less than ~4Ga. We fit a model of these intrusions and the resulting thermal demagnetization to maps of crustal magnetic field strength at 185km altitude. The best fits are most consistent with a "dog bone"-shaped region of intrusive material, elongated approximately north-south, with an area of ~350,000km 2 and an inferred volume of ~4-19×10 6 km 3 . Such a large volume is best explained by a long-lived mantle plume beneath the Syrtis edifice. A free-air gravity anomaly high over the Syrtis Major caldera is consistent with dense mafic residue remaining at depth following crystal fractionation that produced the silicic magmas seen at the surface. The elongation of this region is consistent with ascent and north-south emplacement of magma enabled by structures parallel to and associated with the preexisting Isidis impact basin.

58 GEOSCIENCES↗

Global reconstructed daily surge levels from the 20th Century Reanalysis (1871–2010)

Studying the effect of global patterns of wind and pressure gradients on the sea level variation (storm surge) is a key issue in understanding the recent climate change effect on the dynamical state of the ocean. The analysis of the spatial and temporal variability of storm surges from observations is a difficult task to accomplish since observations are not homogeneous in time, scarce in space, and moreover, their temporal coverage is limited. A recent global surge database developed by AVISO (DAC, Dynamic Atmospheric Correction) fulfilled the lack of data in terms of spatial coverage, but not regarding time extent, since it only includes the last two decades (1992–2014). In this work, we use the 20th Century Reanalysis V2 (20CR), which spans the years 1871 to 2010, to statistically reconstruct daily maximum surge levels at a global scale. A multivariate linear regression model is fitted between daily mean ERA-interim sea level pressure fields and daily maximum surge levels from DAC. Following, the statistical model is used to reconstruct daily surges using mean sea level pressure fields from 20CR. The verification of the statistical model shows good agreements between DAC levels and the reconstructed surge levels from the 20CR. The validation of the reconstructed surge with tide gauges, distributed throughout the domain, shows good accuracy both in terms of high correlations and small errors. A time series comparison is also depicted at specific tide gauges for the beginning of the 20th century, showing a high concordance. Therefore, this work provides to the scientific community, a daily database of maximum surge levels; which correspond to an extension of the DAC database, from 1871 to 2010. This database can be used to improve the knowledge on historical storm surge conditions, allowing the study of their temporal and spatial variability.

Geology↗

Pseudorapidity distributions of charged particles in pp($\bar p$), p(d)A and AA collisions using Tsallis thermodynamics

The pseudorapidity distributions of charged particles measured in p + p($\overline{\text{p}}$) collisions for energies ranging from $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=23.6$ GeV to 13 TeV and A + A collisions at RHIC and LHC are investigated in the fireball model with Tsallis thermodynamics. We assume that the rapidity axis is populated with fireballs following q-Gaussian distribution and the charged particles follow the Tsallis distribution in the fireball. We also extend the fireball model to asymmetric collision systems, i.e. d + Au collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=200$ GeV and p + Pb collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=5.02$ TeV, by taking into account the asymmetric geometry configuration. The model can fit the experimental data well for all the collision systems and centralities investigated. The collision energy and centrality dependence of the model parameters for the symmetric (asymmetric) collision system, i.e. the central position y 0 (y 0a , y 0A ) and its width σ (σ a , σ A ) of the fireball distribution, are also investigated and discussed. Furthermore, our results suggest that the fireball model with Tsallis thermodynamics can be used as a universal framework for the pseudorapidity distributions of charged particles in high energy collisions at RHIC and LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The magnetic field of Mercury

The paper examines the magnetic field observations and their analyses relating to the determination of the Mercury magnetic field. Methods of analyzing data included: (1) comparison of bow shock and magnetopause relative positions at Mercury to the earth, (2) direct spherical harmonic analysis, (3) magnetosphere modeling by an image dipole, and (4) scaling of a mathematical model for the terrestrial magnetosphere. Dipole moments were determined using partial quadrupole and octupole terms to improve the least-square fit of models to observations; analyses by method (2) yield a convergent series of dipole moments values considered to best represent the intrinsic planetary field. Finally, it is suggested that the origin of the magnetic field of Mercury cannot be uniquely determined, but the sources of convective energy may be radiogenic decay and heat release, gravitational settling, and differentiation of processional torques.

Ness, N. F.↗

Multiple-angle observations of reflectance anisotropy from an airborne linear array sensor

An airborne pointable imaging multispectral linear array sensor has been developed for the multidirectional observation of surface reflectance anisotropy. The sensor design permits observations up to 45 deg off-nadir in three spectral bands (green, red, and near-infrared). Calibration permits the conversion of sensor data to radiance units with an absolute uncertainty of 6 percent. Observations of five field plots from seven view directions are discussed. Calibration and atmospheric corrections are used to derive hemispherical-directional reflectance factors. A three-term reflectance model is fit to the reflectance factors for each plot to represent the continuous distribution of reflectance factors with view direction. The reflectance model is integrated over all view directions to calculate bihemispherical reflectance factors. The calculated bihemispherical factors differed by 1 to 25 percent from values based on an assumption of isotropic reflectance depending on spectral band and field plot. These calculations demonstrate the technologic and scientific capabilities required for the remote characterization of surface reflectance anisotropy. Remote multidirectional observations are both feasible and needed to fully evaluate land reflectance characteristics.

Irons, James R.↗

On the Statistical Calibration of Physical Models

ABSTRACT We introduce a novel statistical calibration framework for physical models, relying on probabilistic embedding of model discrepancy error within the model. For clarity of illustration, we take the measurement errors out of consideration, calibrating a chemical model of interest with respect to a more detailed model, considered as “truth” for the present purpose. We employ Bayesian statistical methods for such model‐to‐model calibration and demonstrate their capabilities on simple synthetic models, leading to a well‐defined parameter estimation problem that employs approximate Bayesian computation. The method is then demonstrated on two case studies for calibration of kinetic rate parameters for methane air chemistry, where ignition time information from a detailed elementary‐step kinetic model is used to estimate rate coefficients of a simple chemical mechanism. We show that the calibrated model predictions fit the data and that uncertainty in these predictions is consistent in a mean‐square sense with the discrepancy from the detailed model data.

Sargsyan, K.↗

Multiwavelength investigations of PKS 2300–18: S-shaped radio quasar with precessing jets and double-peaked broad emission-line spectrum

ABSTRACT S-shaped radio galaxy jets are prime sources for investigating the dynamic interplay between the central active galactic nucleus (AGN), the jets, and the ambient intergalactic medium. These sources are excellent candidates for studying jet precession, as their S-shaped inversion symmetry strongly indicates underlying precession. We present a multiwavelength analysis of the giant inversion-symmetric S-shaped radio galaxy PKS 2300$-$18, which spans 0.76 Mpc. The host is a quasar at a redshift of 0.128, displaying disturbed optical morphology due to an ongoing merger with a companion galaxy. We conducted a broad-band radio spectral study using multifrequency data ranging from 183 MHz to 6 GHz, incorporating dedicated observations with the upgraded Giant Metrewave Radio Telescope (uGMRT) and Karl G. Jansky Very Large Array (JVLA) alongside archival radio data. A particle injection model was fitted to the spectra of different regions of the source to perform ageing analysis, which was supplemented with a kinematic jet precession model. The ageing analysis revealed a maximum plasma age of $\sim$ 40 Myr, while the jet precession model indicated a precession period of $\sim$ 12 Myr. ROentgen SATellite (ROSAT) data revealed an X-ray halo of Mpc size, and from Chandra the AGN X-ray spectrum was modelled using thermal and power-law components. The optical spectrum displaying double-peaked broad emission lines was modelled, indicating complex broad-line region kinematics at the centre with the possibility of a binary SMBH. We present the results of our multiwavelength analysis of the source, spanning scales from a few light-days to a few Mpc, and discuss its potential evolutionary path.

Misra, Arpita↗

Testing Lepton Flavor Universality with Pion, Kaon, Tau, and Beta Decays

We present an overview of searches for violation of lepton flavor universality with a focus on low energy precision probes using π, K, τ, and nuclear beta decays. We review the current experimental results, summarize the theoretical status within the context of the Standard Model, and discuss future prospects (both experimental and theoretical). We review the implications of these measurements for physics beyond the Standard Model by performing a global model-independent fit to modified W couplings to leptons and four-fermion operators. We also discuss new physics in the context of simplified models and review Standard Model extensions with a focus on those that can explain a possible deviation from unitarity of the Cabibbo–Kobayashi–Maskawa quark mixing matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantile Kriging R Package

A re-implementation of Quantile Kriging. Quantile Kriging was described by Matthew Plumlee & Rui Tuo in their 2014 paper "Building Accurate Emulators for Stochastic Simulations via Quantile Kriging." With computational savings when dealing with replication from the recent paper "Practical heteroskedastic Gaussian process modeling for large simulation experiments”, it is now possible to apply Quantile Kriging to a wider class of problems. In addition to fitting the model, other useful tools are provided such as the ability to automatically perform leave-one-out cross validation.

Quinlan, KevinR.↗