Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,171 records · Page 65

Tropopause Characteristics Based on Long-Term ARM Radiosonde Data: A Fine-Scale Comparison at the Extratropical SGP Site and Arctic NSA Site

The variations in the characteristics of the tropopause are sensitive indicators for the climate system and climate change. By using Atmospheric Radiation Measurement (ARM) radiosonde data that were recorded at the extratropical Southern Great Plains (SGP) and Arctic North Slope of Alaska (NSA) sites over an 18-year period (January 2003 to December 2020), this study performs a fine-scale comparison of the climatological tropopause features between these two sites that are characterized by different climates. The static stability increases rapidly above the tropopause at both sites, indicating the widespread existence of a tropopause inversion layer. The structures of both the tropopause inversion layer and the stability transition layer are more obvious at NSA than at SGP, and the seasonal variation trends of the tropopause inversion layer and stability transition layer are distinctly different between the two sites. A fitting method was used to derive the fitted tropopause height and tropopause sharpness (λ). Although this fitting method may determine a secondary tropopause rather than the primary tropopause when multiple tropopause heights are identified on one radiosonde profile, the fitted tropopause heights generally agree well with the observed tropopause heights. Broad tropopause sharpness values (λ > 2 km) occur more frequently at SGP than at NSA, resulting in a greater average tropopause sharpness at SGP (1.0 km) than at NSA (0.6 km). Significant positive trends are exhibited by the tropopause heights over the two sites, with rates of increase of 23.7 ± 6.5 m yr ⁻1 at SGP and 28.0 ± 4.0 m yr ⁻1 at NSA during the study period.

54 ENVIRONMENTAL SCIENCES↗

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Statistical Characterization of Temperature and Pressure Vertical Profiles for the Analysis of Laser Heterodyne Radiometry Data

The statistical analysis of historic pressure and temperature profiles from radiosonde launches for use in the fitting of molecular oxygen line shapes is presented. As the O 2 mixing ratio is nearly constant throughout the lower atmosphere, only variations in pressure and temperature profiles will affect the fit of observed O 2 features in Laser Heterodyne Radiometry (LHR) spectra. Radiosonde temperature and pressure data are extracted from the Integrated Global Radiosonde Archive (IGRA) for a given station, date, and launch time. Data may be extracted for a single launch, for the same date over several years, and/or within a window centered on a target date. The temperature and pressure profiles are further characterized by the statistical variation in coefficients of polynomial fits in altitude. The properties of the probability distributions for each coefficient are used to constrain fits of O 2 line shapes through Nelder–Mead optimization. The refined temperature and pressure profiles are then used in the retrieval of vertically resolved mixing ratios for greenhouse gases (GHGs) measured in the same instrument. In continuous collections, each vertical profile determination may be treated as a Bayesian prior to inform subsequent measurements and provide an estimate of uncertainties.

47 OTHER INSTRUMENTATION↗

Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood

Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.

Astrostatistics (1882)↗

The Paschen Jump as a Diagnostic of the Diffuse Nebular Continuum Emission in Active Galactic Nuclei*

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope Space Telescope Imaging Spectrograph spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multicomponent spectral models over the range 6800–9700 Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from ~10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe ii emission, high-order Balmer line emission, and other components.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar Population Properties in the Stellar Streams around SPRC047

Abstract We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z= 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6μm data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6μm IRAC image, thereby reducing the effect of the point-spread function aliasedemissionfrom the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars.

Astronomy & Astrophysics↗

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H 2 abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ as the formation time of H 2 becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected ƒ H 2 measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy, the overall molecular gas mass is reproduced to better than a factor of ≲1.5 across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous ƒ H 2 models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

Chemical Cartography with APOGEE: Mapping Disk Populations with a 2-process Model and Residual Abundances

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] versus [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core-collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01–0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow for the comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R = 3–13 kpc), we find nearly identical abundance patterns in the outer disk (R = 15–17 kpc), 0.05–0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4–1 dex) of multiple elements in ω Cen. The residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of an adhesively bonded beryllium propulsion structure for the Mariner Mars 1971 spacecraft

The design, testing, and fabrication of the support truss structure for the propulsion system of the Mariner 9 spacecraft are described. Support is provided by an 8.9-kg (19.5-lbm) truss assembly consisting of beryllium tubes adhesively bonded to magnesium end fittings. Beryllium was selected for the tubular struts in the truss because of its exceptionally high stiffness-to-weight ratio. Adhesive bonding, rather than riveting, was utilized to join the struts to the end fittings because of the low toughness (high notch sensitivity) of beryllium. Magnesium, used in the end fittings, resulted in a 50% weight saving over aluminum since geometric factors in the fitting design resulted in low stress areas where magnesium's lower density is a benefit.

Stevens, J. H.↗

Information content of IRIS spectra

Spectra from the satellite instrument IRIS (infra red interferometer spectrometer) were examined to find the number of independent variables needed to describe these broadband high spectral resolution data. The radiated power in the atmospheric window from 771 to 981/cm was the first parameter chosen for fitting observed spectra. At succeeding levels of analysis the residual variability (observed spectrum - best fit spectrum) in an ensemble of observations was partioned into spectral eigenvectors. The eigenvector describing the largest fraction of this variability was examined for a strong spectral signature; the power in the corresponding spectral band was then used as the next fitting parameter. The measured power in nine spectral intervals, when inserted in the spectral fitting functions, was adequate to describe most spectra to within the noise level of IRIS. Considerations of relative signal strength and scales of atmospheric variability suggest a combination sounder (multichannel-broad field of view) scanner (window channel-small field of view) as an efficient observing instrument.

Price, J. C.↗

The 6300 A O/1-D/ airglow and dissociative recombination

Measurements of night-time 6300 A airglow intensities at the Arecibo Observatory have been compared with dissociative recombination calculations based on electron densities derived from simultaneous incoherent backscatter measurements. The agreement indicates that the nightglow can be fully accounted for by dissociative recombination. The comparisons are examined to determine the importance of quenching, heavy ions, ionization above the F-layer peak, and the temperature parameter of the model atmosphere. Comparable fits between the observed and calculated intensities are found for several available model atmospheres. The least-squares fitting process, used to make the comparisons, produces comparable fits over a wide range of combinations of neutral densities and of reaction constants. Yet, the fitting places constraints upon the possible combinations; these constraints indicate that the latest laboratory chemical constants and densities extrapolated to a base altitude are mutually consistent.

Wickwar, V. B.↗

Information content in Iris spectra

Spectra from the satellite instrument Iris (infrared interferometer spectrometer) were examined to find the number of independent variables needed to describe the broad-band high-resolution spectral data. The radiated power in the atmospheric window from 771 to 981 per cm was the first parameter chosen for fitting observed spectra. At succeeding levels of analysis, the residual variability (observed spectrum minus best-fit spectrum) in an ensemble of observations was partitioned into spectral eigenvectors. The eigenvector describing the largest fraction of this variability was examined for a strong spectral signature; the power in the corresponding spectral band was then used as the next fitting parameter. The measured power in nine spectral intervals, when it was inserted in the spectral-fitting functions, was adequate to describe most spectra to within the noise level of Iris. Considerations of relative signal strength and scales of atmospheric variability suggest a combination sounder (multichannel, broad field of view) scanner (window channel, small field of view) as an efficient observing instrument.

Price, J. C.↗

Rotary leveling base platform

A leveling apparatus for the precise adjustment of a scientific instrument is reported. A base member is provided having a hollow cylindrical shape. A table for supporting the instrument rests on the base and has a shaft portion extending below the table. The upper portion of the shaft fits tightly into the hollow portion of the base member whereas the lower portion of the shaft is machined to fit loosely. The lower portion of the shaft is provided with a groove. Adjusting screws are threaded through the hollow cylindrical portion and are adapted to enter the groove. By adjusting the screws, the lower portion of the shaft is moved in a vertical plane since the shaft is loosely fitted into the cylinder. The upper portion of the shaft which is tightly fitted into the upper end of the cylinder causes the cylinder to deform slightly providing a fulcrum point which allows the table to be leveled in response to the adjustment of the adjusting screws.

Delaplaine, R. W.↗

The Construction of Curves and Surfaces Using Numerical Optimization Techniques

Numerical optimization techniques are playing an increasing role in curve and surface construction. Often difficult problems in curve and surface construction, especially when some aspect of shape control is involved, can be phrased as a constrained optimization problem. Four such classes of problems are explored: parametric curve fitting with non-linear shape constraints; explicit surface fitting with linear shape constraints; surface fitting to scattered data giving rise to ill-posed problems; finally, variable knot problems. In each of these problems there is a nonlinear aspect: either the shape of the curve or surface is important for manufacturing or engineering reasons or the shape affects the convergence of numerical algorithms which use the curve or surface or the placement of knots affects the accuracy of the fits. In all cases the class of functions used is that of parametric spline curves and tensor or direct product spline surfaces. The reason for choosing this class is that splines provide flexible models that are easily evaluated and stored. Furthermore, the B-spline representation of splines leads to convenient expressions for shape control over regions.

Ferguson, D. R.↗

Saskatoon M. F. Radar System: (52 Deg N, 107 Deg W), Canada

The system runs continuously, producing 1-h profiles approx. 75 to 110 km, 1978 to 1983: in particular for the 30-d (Nov. 9 to Dec. 7), 10-d (Nov. 16 to 26), 4-d core periods (Nov. 19 to 23). The means (EW, NS) and fourier components (24-, 12-h) fitted for the 4-d are shown. s.d. are shown, as derived from four 24-h fits in the interval. Comparisons with the 10-, 30-d intervals confirm confidence in these data and the lack of major variability in November. The 1-h profiles are means of 12 profiles so intrinsic errors are minimal. The 4-, 10-d fourier fits are compared and also the 30-d means from 48-h fits; in these amplitudes are N2+E2 as N/E approx. 1. Means differ slightly due to planetary wave activity; 24-h tides are variable in Novembers; and the 12-h variation is small less than 80 km. The final figure is a high resolution spectra for the November month: the 12-, 24-h dominate; the 8-h appears less than 100 km; and up to 90 km there is a 4-d oscillation.

Manson, A. H.↗

The X-ray spectrm of the Coma Cluster of galaxies

The HEAO 1 A-2 and OSO 8 observations of the Coma Cluster of galaxies in the 2-60 keV energy band are presented. The spectrum is not well fitted by an isothermal bremsstrahlung model with an Fe line. To model the spectrum more accurately, we have performed fits of polytropic models and find that the gas is well fitted by an equation of state with polytropic index ranging between 1.38 and 1.6. This analysis reduces the total cluster mass by a factor of 3-5 below the optically derived virial mass. Using the best-fitting parameters, the gas is found to be a significant fraction (25 percent) of the total cluster mass. The Fe K-alpha line blend is strongly detected, from which we infer 0.18-0.36 solar abundance in the X-ray emitting gas. The Fe K-beta line blend is also weakly detected. The spectrum is not consistent with the existence of a hard X-ray tail at the flux level reported by Bazzano el al. (1984). In addition, the HEAO I scanning mode data place strong constraints on possible contaminating sources in the region surrounding Coma.

Henriksen, M. J.↗

Modeling of coronal X-ray emission from active cool stars. I Hyades cluster

X-ray pulse height spectra of the most active cool stars in the Hyades cluster obtained with the Einstein IPC cannot be satisfactorily fitted using isothermal thin plasma emission models. Addition of a second isothermal component provides acceptable fits. However, a more physically meaningful set of coronal parameters is provided by models which consist of an ensemble of loops wih a single maximum temperature, but with the temperature distribution within the loop determined by loop physics. Such models have been successfully fitted to the IPC pulse height spectra. Constraints on loop parameters are discussed for the F-G dwarfs BD + 14 deg 693, BD + 14 deg 690, BD + 15 deg 640, and 71 Tau. Models with a large variation of loop cross section from base to top do not fit the data. A consistent physical description is an ensemble of small high-pressure loops of similar maximum temperature which dominate the coronal X-ray spectrum.

Stern, R. A.↗

A spectral study of the persistent X-ray flux from 4 U/MXB 1636-53

A spectral analysis of the persistent X-ray flux from the X-ray burst source 4 U/MXB 1636-53, which was continuously observed with EXOSAT for about 25 h in May, 1984, is presented. It is found that no standard one-component spectral model provides an acceptable representation of all the data. However, almost all 'complex one-component' and 'two-component' models do accurately fit the data. A blackbody plus a 'multicolor' disk spectrum does not fit these data. In addition, some 'two-component' spectral models can be fitted to the data with more than one set of different parameter values. Based on these results, it is suggested that the various physical scenarios that have been proposed to interpret 'two-component' spectra of low-mass X-ray binaries may be based on mathematical parameterizations which may not provide unique representations of the data. If this is the case, the reported 'two-component' mathematical fits may not have any true physical significance.

Vacca, W. D.↗