Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Curve 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 973 records · Page 54

A new world survey expression for cosmic ray vertical intensity vs. depth in standard rock

The cosmic ray data on vertical intensity versus depth below 10 to the 5th power g sq cm is fitted to a 5 parameter empirical formula to give an analytical expression for interpretation of muon fluxes in underground measurements. This expression updates earlier published results and complements the more precise curves obtained by numerical integration or Monte Carlo techniques in which the fit is made to an energy spectrum at the top of the atmosphere. The expression is valid in the transitional region where neutrino induced muons begin to be important, as well as at great depths where this component becomes dominant.

Crouch, M.↗

Candidate strongly lensed type Ia supernovae in the Zwicky Transient Facility archive

Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. The candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 ± 3 and 34 ± 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, such as the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints On the Emission Geometries and Spin Evolution Of Gamma-Ray Millisecond Pulsars

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase, we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using amaximum likelihood technique.We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II), or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of the magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

pulsars: general↗

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗

Second ARM Aerosol Chemical Speciation Monitor Users’ Meeting Report

The aerosol chemical speciation monitor (ACSM) was developed to adapt the technology of the aerosol mass spectrometer (AMS) to routine, long-term, standalone monitoring. The calculation of particulate mass concentration from ACSM data requires the measurement of the response of the instrument to aerosol of specific size and composition as well as assumptions about the instrument response based on laboratory measurements and field experience acquired over more than two decades of operation of AMS and a decade of operation of the ACSM. Three parameters in the concentration calculations that are particularly important are the NO3 response factor (RFNO3), relative ionization efficiency (RIE), and the collection efficiency (CE). The values of RFNO3 and RIEs are determined from calibration, however the jump calibration method previously used in calibration can result in errors in the RIE for sulfate. This has been corrected by implementing a continuous calibration method. The default collection efficiency is 0.5. This has been shown to result in mass loadings that do not agree with mass determinations from other instruments because of effects of composition on the vaporization of the particles. The previous work of investigators addressing this issue is discussed. After preliminary work on ACSM and scanning mobility particle sizer (SMPS) data from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) observatory collected in late 2016 and 2017 produced a parameterization of composition dependent collection efficiency very different from the results of previous studies, SMPS data were examined and we determined that there was significant mass that the instrument did not capture because the particles with diameters larger than 465 nm are not counted by this instrument. Data for the ultra-high-sensitivity aerosol spectrometer (UHSAS), SMPS, and ACSM are available for nearly all of 2019. The data from UHSAS and SMPS collected in 2019 were compared. We found that the UHSAS data has particle counts and total volumes significantly less than measured by the condensation particle counter (CPC) and SMPS. The SMPS data were extended by fitting the average volume distribution with a log normal curve and using this relationship to estimate a mass value for the SMPS over extended diameter range. The extended SMPS mass values result in a CDCE parameterization that is in better agreement with the results of other investigators, but is still different from other formulations. The working group recommends that the ACSM data be processed with a collective efficiency (CE)=1, that this be documented clearly in the metadata, and the use of the default CE of 0.5 or a formulation of composition-dependent collection efficiencies (CDCE) chosen by the user should be implemented based on the ammonium nitrate mass fraction. This is clearly necessary for the wintertime SGP ACSM data because of the high nitrate concentrations.

47 OTHER INSTRUMENTATION↗

Observations of the Orion nebula at 100 microns

The Orion nebula spectrum was obtained using an airborne cooled grating spectrometer in the spectral range from 80 to 125 microns at a resolution of 8 microns. The spectrum, corrected for instrument profile and atmospheric water vapor absorption, could be fit by the 60 and 100 K black body curves. The results are in good agreement with those of Erickson et al. (1973) and are consistent with a thermal emission spectrum produced by cool interstellar grains radiating at a wavelength much longer than the grain size.

Ward, D. B.↗

The X-ray source Serpens X-1 - Ariel 5 observations and discussion of models for the spectrum and time variability

Experimental data from Ariel 5 on the spectral shape and time variability of Ser X-1 are presented, and possible explanations for the observations are discussed in terms of current theoretical suggestions for source emission. The observations are summarized in the form of a light curve for 3-7.6-keV photons. The data are fitted with a power law of index -2.3, which yields a hydrogen column density of (1.1 + or - 0.4) x 10 to the 22nd power atoms/sq cm. No persistent periodicity of amplitude greater than 5% of the steady flux is found, but evidence of statistically significant burst activity is obtained. Various emission mechanisms are considered for the time-averaged spectrum and the X-ray bursts. It is suggested that the inverse Compton mechanism is a likely cause for the emission from this source and that the source is radiating near the Eddington limit.

Coe, M. J.↗

The x ray variability of AGN and the anomalous behavior of NGC6814

The power spectra of the X-ray variable AGN are typically scale invariant with no characteristic timescale. The one strong exception to this is NGC6814, where the EXOSAT data showed evidence for a periodic component at 12200 plus or minus 100 seconds. The power spectra of a GINGA lightcurve from this source, found using simulation techniques to account for the uneven sampling, also cannot be well fit by a single power law. A folded light curve analysis of the GINGA data shows a strong peak in the variance, indicative of a periodic component, at 12130.9 plus or minus 0.6 seconds. A second GINGA observation of this source taken one year later is consistent with the periodicity maintaining phase coherence for 7 periods in the range of 12110-12145 seconds. Including the limits from the folding selects a period of 12130.39 plus or minus 0.05 seconds. Phase coherence is not maintained between this and the EXOSAT observations, as the structure of the folded light curves is very different. Thus the periodicity is long lived and stable, but phase coherence is only maintained on timescales of approximately 1 year.

Done, C.↗

The X-ray variability of AGN and the anomalous behavior of NGC 6814

The power spectra of the X-ray variable AGN are typically scale invariant with no characteristic timescale. The one strong exception to this is NGC 6814, where the EXOSAT data showed evidence for a periodic component at 12200 +/- 100 seconds. The power spectra of a GINGA lightcurve from this source, found using simulation techniques to account for the uneven sampling, also cannot be well fit by a single power law. A folded light curve analysis of GINGA data shows strong evidence for this periodic component. A second GINGA observation of this source taken one year later is consistent with the phenomena being completely periodic and phase coherent for 7 periods in the range of 12,110-12,145 seconds. Including the (optimistic) limits from the folding selects a period of 12130.39 +/- 0.05 seconds. Phase coherence is not maintained between this and the EXOSAT observations, as the structure of the folded light curves is very different. Thus the periodicity is long lived and stable, but phase coherence is only maintained on timescales of about 1 year.

Done, C.↗

Metrics of Balance Control for Use in Screening Tests of Vestibular Function

Decrements in balance control have been documented in astronauts after space flight. Reliable measures of balance control are needed for use in postflight field tests at remote landing sites. Diffusion analysis (DA) is a statistical mechanical tool that shows the average difference of the dependent variable on varying time scales. These techniques have been shown to measure differences in open-loop and closed-loop postural control in astronauts and elderly subjects. The goal of this study was to investigate the reliability of these measures of balance control. Eleven subjects were tested using the Clinical Test of Sensory Interaction on Balance: the subject stood with feet together and arms crossed on a stable or compliant surface, with eyes open or closed and with or without head movements in the pitch or yaw plane. Subjects were instrumented with inertial motion sensors attached to their trunk segment. The DA curves for linear acceleration measures were characterized by linear fits measuring open- (Ds) and closed-loop (Dl) control, and their intersection point (X-int, Y-int). Ds and Y-int showed significant differences between the test conditions. Additionally, Ds was correlated with the root mean square (RMS) of the signal, indicating that RMS was dominated by open-loop events (< 0.5 seconds). The Y-int was found to be correlated with the average linear velocity of trunk movements. Thus DA measures could be applied to derive reliable metrics of balance stability during field tests.

Fiedler, Matthew↗

Panchromatic Observations of the Textbook GRB 110205A: Constraining Physical Mechanisms of Prompt Emission and Afterglow

We present a comprehensive analysis of a bright, long duration (T(sub 90) approx. 257 s) GRB 110205A at redshift z = 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Thanks to its long duration, nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray (1 eV - 5 MeV), which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. In particular, by fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/ -ray spectra, it traces the -ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + synchrotron self-Compton (SSC) scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (approx. 1100 s), a bright (R = 14.0) optical emission hump with very steep rise ( alpha approx. 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed.

Zheng, W.↗

AG Channel Measurement and Modeling Results for Over-Sea Conditions

This report describes results from flight tests conducted in an over-sea environment, for the purpose of characterizing the air-to-ground (AG) channel, for future unmanned aircraft system (UAS) communication system analysis and design. These results are for the first of a set of several flight tests conducted in different ground site (GS) environments. An ultimate aim of all these tests is the development of models for the AG channel that can be used in communication system evaluation. In this report we provide measured results for propagation path loss, root-mean square delay spread (RMS-DS), and the correlation coefficient of the primary received signal components on the four antennas (two antennas for C-band, two for L-band). For path loss, the curved-earth two-ray model provides a reasonable fit to the measured data, altered by several dB at the shortest link distances by aircraft antenna pattern effects. This two-ray model also accounts for the majority of measured RMS-DS results of a few tens of nanoseconds, except for the occasional intermittent reflections from surface objects. These intermittent reflections yield RMS-DS values up to several hundred nanoseconds. For portions of the flight path that were over a harbor area highly populated with boats, the channel was found to be more "continuously dispersive," with RMS-DS reaching approximately 250 ns. A separate model will be developed for this over-harbor setting. The correlation coefficient results are still undergoing analysis; preliminary observations are that correlation between signals on the same-band antennas is generally large (>0.6) for the C-band straight flight paths, whereas for the L-band signals and for the oval-shaped flight paths the correlation is generally small (below 0.4). Inter-band correlations are typically very small, and are well modeled as zero-mean Gaussian in distribution, with a standard deviation less than 0.2. Hence the over-sea channel effects in the two bands can be considered uncorrelated, which will allow for good diversity gains in dual-band systems. We describe initial modeling approaches for the over-sea channel; complete models for this and the over-harbor setting will appear in a subsequent report.

Aircraft Communications↗

Earthshine as an Illumination Source at the Moon

Earthshine is the dominant source of natural illumination on the surface of the Moon during lunar night, and at some locations within permanently shadowed regions (PSRs) near the poles that never receive direct sunlight. As such, earthshine has the potential to enable the scientific investigation and exploration of conditions in areas of the Moon that are either temporarily or permanently hidden from the Sun. Earthshine has also been used to refer to Earthlight reflected from the lunar surface, but in this study we use it to refer specifically to Earthlight incident at the Moon. Under certain circumstances, the heat flux from earthshine could also influence the transport and cold-trapping of volatiles present in the very coldest areas within PSRs. In this study, Earth's spectral irradiance, as it would appear at the Moon in the solar reflectance band (0.3–3.0 µm) and at thermal emission wavelengths (3–50 µm), is examined with a suite of model image cubes and whole-disk spectra created using the Virtual Planetary Laboratory (VPL) three-dimensional (latitude, longitude and altitude) modeling capability. At the Moon, the broadband, hemispherical irradiance from Earth at full-phase is approximately 0.15 W m−2 with comparable contributions from solar reflectance and thermal emission; for context, this about 0.01% that of solar irradiance and has an equivalent temperature of around 40 K. Over the simulated timeframe, spanning two lunations, Earth's thermal irradiance shows very little net change (less than a few mW m−2 resulting from cloud variability and the south-to-north motion of the sub-observer latitude on Earth). In the solar band, Earth's diurnally averaged light curve at phase angles g ≤ 60° is well-fit using a Henyey–Greenstein integral phase function. At wavelengths longward of about 0.7 µm, near the well-known vegetation “red edge”, Earth's reflected solar radiance shows significant diurnal modulation as a result of the broad maximum in projected landmass at terrestrial longitudes between 60°W and 0°, as well as from the distribution of clouds. A simple formulation with adjustable coefficients is presented, condensed from the VPL model grid, for estimating Earth's hemispherical irradiance at the Moon as a function of wavelength, phase angle and sub-observer coordinates (terrestrial latitude and longitude). Uncertainties in any one prediction are estimated to be 10–12% at 0.3 µm, rising to >25% near 2.5 µm as a result of the increasing relative brightness and unpredictable influence of clouds. Although coefficient values are derived from a suite of spring equinox models, the approximation appears to be valid for all seasons, to within the stated uncertainties. It is demonstrated that earthshine is sufficiently bright to serve as a natural illumination source for optical measurements on a robotic lander/rover, allowing the identification of water ice mixed with regolith at the percent-level of fractional area.

Surface↗

Future Intensity‐Duration‐Frequency Curves of Extreme Precipitation in the Midwest United States From Convection‐Permitting Modeling

Abstract During the last four decades, global warming has statistically significant intensified extreme precipitation events in the Midwestern United States (defined here as the region covering Illinois, Indiana, Ohio, and Kentucky), leading to increased risks to human life, property, and infrastructure. To enable climate change adaptation and resilience across various economic and social sectors in this region, updated information about future climate changes, specifically at finer spatial scales, is essential. Leveraging a new 150‐year dynamical downscaling data set at convection‐permitting resolution, this study introduces a framework to construct the projected future intensity‐duration‐frequency (IDF) curves of heavy precipitation, which are prominent tools for infrastructure design and water resources management. This framework generates IDF curves at both sub‐daily and multi‐day duration utilizing hourly in situ observations as well as quantile‐based statistical techniques in bias‐correction and return levels selection. The assumption of non‐stationarity in the distribution parameter fitting process is also implemented in this workflow. Compared to historical IDF curves for 1980–2022, future projected IDF curves for 2058–2100 under Representative Concentration Pathway (RCP) 4.5 and RCP 8.5 scenarios indicate an average intensity increase of approximately 15% and 25%, respectively, across 74 stations, considering both annual and seasonal timescales. Future projections suggest that extreme precipitation events may become more severe across six investigated return periods, with longer return periods showing a greater increase. The frequency of future extreme precipitation events in the Midwest region is also projected to double. Furthermore, current results reveal spatial heterogeneity of future trends across stations owing to the high‐resolution input data set. Plain Language Summary This study investigates the evolving nature of extreme precipitation events in the Midwestern United States under a changing climate. By leveraging a high‐resolution dynamical downscaling data set, we construct projected intensity‐duration‐frequency (IDF) curves for future extreme rainfall events. These curves serve as vital tools for infrastructure planning and water resource management. Our analysis reveals a significant increase in both the intensity and frequency of extreme precipitation events in the region. Future projected IDF curves for the late century indicate an average intensity increase of approximately 15%–25% compared to historical values. Moreover, the frequency of such events is expected to double. Spatial heterogeneity in future trends is observed across different stations within the Midwest, highlighting the importance of high‐resolution modeling in capturing localized climate variability. These findings underscore the urgent need for climate adaptation strategies to mitigate the increasing risks associated with extreme precipitation events in the region. Key Points This study introduces a workflow to construct future intensity‐duration‐frequency (IDF) curves over the Midwest United States using a new convection‐permitting modeling data set The current IDF construction workflow reproduces well the historical observed IDF 30 curves in summer months with median relative errors of 2.4% among 74 stations and 6 investigated durations The projected IDF curves show diverse future trends of extreme precipitation across stations, with intensity increases of approximately 15% and 25% under RCP4.5 and RCP8.5 climate scenarios, respectively, and a doubling of frequency on average

Nguyen, Trung↗

Cloud Atlas: Rotational Modulations in the L/T Transition Brown Dwarf Companion HN Peg B

Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 (WFC3) near-infrared G141 taken in six consecutive orbits observations of HNPeg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the 1.1 to 1.7 micron broadband light curve has the amplitude of and period of hour. The modulation amplitude has no detectable wavelength dependence except in the 1.4 micron water absorption band, indicating that the characteristic condensate particle sizes are large (greater than 1 micron). We detect significantly (4.4 sigma) lower modulation amplitude in the 1.4 micron water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3.

brown dwarf↗

Generation of surface-fitted coordinate grids

An accurate numerical solution of the Navier-Stokes equations is strongly dependent on the description of the flow geometry. In the past, rectangular grids led to an ill description of curved boundaries and consequently hindered finite-difference solutions. Now, it is accepted that body-fitted coordinate systems give a much better description of complex flow geometries. The generation of body-fitted curvilinear coordinate grids in two dimensions by the solution of a nonlinear system of elliptic equations has been proposed Thompson. It has been extended to three dimensions and applied in practical configurations.

Garon, A.↗

Solar-cell interconnect design for terrestrial photovoltaic modules

Useful solar cell interconnect reliability design and life prediction algorithms are presented, together with experimental data indicating that the classical strain cycle (fatigue) curve for the interconnect material does not account for the statistical scatter that is required in reliability predictions. This shortcoming is presently addressed by fitting a functional form to experimental cumulative interconnect failure rate data, which thereby yields statistical fatigue curves enabling not only the prediction of cumulative interconnect failures during the design life of an array field, but also the quantitative interpretation of data from accelerated thermal cycling tests. Optimal interconnect cost reliability design algorithms are also derived which may allow the minimization of energy cost over the design life of the array field.

Mon, G. R.↗

Kepler Uniform Modeling of KOIs: MCMC Notes for Data Release 25

This document describes data products related to the reported planetary parameters and uncertainties for the Kepler Objects of Interest (KOIs) based on a Markov-Chain-Monte-Carlo (MCMC) analysis. Reported parameters, uncertainties and data products can be found at the NASA Exoplanet Archive . The codes used for this data analysis are available on the Github website (Rowe 2016). The relevant paper for details of the calculations is Rowe et al. (2015). The main differences between the model fits discussed here and those in the DR24 catalogue are that the DR25 light curves were used in the analysis, our processing of the MAST light curves took into account different data flags, the number of chains calculated was doubled to 200 000, and the parameters which are reported are based on a damped least-squares fit, instead of the median value from the Markov chain or the chain with the lowest 2 as reported in the past.

DR25↗