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 721 records · Page 40

Constraints on the Recent Rate of Lunar Regolith Accumulation from Diviner Observations

Many large craters on the lunar nearside show radar CPR signatures consistent with the presence of blocky ejecta blankets, to distances pre dicted to be covered by continuous ejecta. However, most of these sur faces show limited enhancements in both derived rock abundance and rock-free regolith temperatures calculated from Diviner nighttime infrar ed observations. This indicates that the surface rocks are covered by a layer of thermally insulating regolith material. By matching the results of one-dimensional thermal models to Diviner nighttime temperat ures, we have constrained the thermophysical properties of the upper regolith, and the thickness of regolith overlying proximal ejecta. We find that for all of the regions surveyed (all in the nearside highla nds), the nighttime cooling curves are best fit by a density profile that varies exponentially with depth, consistent with a linear mixture of rocks and regolith fines, with increasing rock content with depth . Our results show significant spatial variations in the density e-folding depth, H, among young crater ejecta regions, indicating differen ces in the thickness of accumulated regolith. However, away from youn g craters, the average regional "equilibrium" value of H (Heq) is remarkably consistent, and is on the order of 5 cm. As expected, near-rim ejecta associated with young craters show lower values of H, indicating a high rock content in the shallow subsurface; for older craters, the average value of H approaches the regional value of Heq. Calculat ed H values for young craters show a clear correlation with published ages, providing the first observational constraint on the recent rate of lunar regolith accumulation. In addition, this result may help to resolve the apparent discrepancy between ages calculated from small crater counts on melt ponds versus counts on continuous ejecta (e.g., King crater; Ashley et al., 2011, LPSC 42, abstract 2437). This method could, in principle, be extended to other airless bodies (e.g., aste roids), which would in turn constrain the recent impactor flux.

Ghent, R. R.↗

Energies of Various Interactions Between Hydrogen and Helium Atoms and Ions

A compilation is made of those interactions between hydrogen and helium atoms and ions whose potential curves are well established by virtue of being based on experimental data or on essentially exact quantum-mechanical calculations. These potential curves have been fitted with well-known empirical forms which should be useful in calculations of transport properties involving these interactions.

Fallon, Robert J.↗

Percolation Behaviour in the Magnetic Permeability and Electrical Conductivity in Conducting Magnetic - Insulating Non Magnetic Binary Composites

Experimental results of the complex magnetic permeability (µ) and the electrical conductivity (σ) of a granular paramagnetic Gadolinium Gallium Garnet (GGG: 0.3–26 vol%) and Teflon (PTFE) system are presented and discussed in relation to previously published (conductivity) and unpublished (permeability) studies on granular Fe3O4 – talc and Ni – talc wax systems. In these systems plots of the real conductivity (σ'(sub m)) against the volume fraction (φ) lie on characteristic sigmoid curves that when fitted to the Two Exponent Phenomenological Percolation Equation (TEPPE) confirm the existence of “percolation microstructures” with critical volume fractions (φ(sub c)). The plots of the real and imaginary permeability (µ'(sub m)) and (µ"(sub m)) satisfactorily fit to the TEPPE using the φ(sub c) obtained in each case from the “conductivity” measurements. In all three cases the conductivity results gave the exponent t > 2, and the permeability results gave t < 1.

David S McLachlan↗

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗

OGLE-2017-BLG-1186: First Application of Asteroseismology and Gaussian Processes to microlensing

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of ∼9 d. The asteroseismic analysis also provides us source properties including the source angular size (∼27 μas) and distance (∼11.5 kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional χ(exp 2) minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a ∼0.045 Mʘ brown dwarf at distance ∼9.0 kpc, or a ∼0.073 Mʘ ultracool dwarf at distance ∼9.8 kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a ∼ 35 per cent probability for the lens to be a bulge object and ∼ 65 per cent to be a background disc object.

S.-S. Li↗

TOI-674b: An oasis in the desert of exo-Neptunes transiting a nearby M dwarf

Context. The NASA mission TESS is currently doing an all-sky survey from space to detect transiting planets around bright stars. As part of the validation process, the most promising planet candidates need to be confirmed and characterized using follow-up observations. Aims. In this article, our aim is to confirm the planetary nature of the transiting planet candidate TOI-674b using spectroscopic and photometric observations. Methods. We use TESS, Spitzer, ground-based light curves, and HARPS spectrograph radial velocity measurements to establish the physical properties of the transiting exoplanet candidate TOI-674b. We perform a joint fit of the light curves and radial velocity time series to measure the mass, radius, and orbital parameters of the candidate. Results. We confirm and characterize TOI-674b, a low-density super-Neptune transiting a nearby M dwarf. The host star (TIC 158588995, V = 14.2 mag, J = 10.3 mag) is characterized by its M2V spectral type with M⋆ = 0.420 ± 0.010 Mꙩ, R⋆ = 0.420 ± 0.013 Rꙩ, and T(eff) = 3514 ± 57 K; it is located at a distance d = 46.16 ± 0.03 pc. Combining the available transit light curves plus radial velocity measurements and jointly fitting a circular orbit model, we find an orbital period of 1.977143 ± 3 × 10^(−6) days, a planetary radius of 5.25 ± 0.17 Rꚛ, and a mass of 23.6 ± 3.3 Mꚛ implying a mean density of ρp =0.91 ± 0.15 g/cu. cm. A non-circular orbit model fit delivers similar planetary mass and radius values within the uncertainties. Given the measured planetary radius and mass, TOI-674b is one of the largest and most massive super-Neptune class planets discovered around an M-type star to date. It is found in the Neptunian desert, and is a promising candidate for atmospheric characterization using the James Webb Space Telescope.

F. Murgas↗

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

Structure evolution of Bi 4 O 7 under high pressure

X-ray diffraction and Raman spectroscopy measurements were used to study the structure evolution of Bi 4 O 7 during compression. Two isostructural phase transitions were observed under high pressure, with the first transition to phase II occurring at approximately 3.0 (2) GPa and the second to phase III occurring at approximately 13.0 (5) GPa. The pressure–volume ( P – V ) curve was found to be discontinuous at those two pressure transition boundaries. By fitting the P – V curves in three different stages, the bulk moduli were estimated to be approximately 101 (2), 68 (2) and 172 (3) GPa. Such large differences in the bulk moduli confirm two isostructural phase changes in the three stages. An irreversible amorphization process was found to start at approximately 20 (1) GPa, and the sample gradually becomes fully amorphous during compression at 30 (1) GPa.

Chemistry↗

An analysis of the shapes of ultraviolet extinction curves. III - An atlas of ultraviolet extinction curves

The basic data needed to produce UV extinction curves for a sample of 78 stars, including spectral types, visual photometry, parameters used to produce fits to the UV extinction curves, distances to the program stars, and line-of-sight H I column densities, are presented. The observational and numerical techniques used to obtain the results are discussed in detail. The UV extinction curves and their analytic fits are shown graphically.

Fitzpatrick, Edward L.↗

Interstellar absorption lines in the spectrum of Gamma Velorum

Copernicus scans of selected interstellar absorption lines in the UV spectrum of Gamma Vel are analyzed, together with ground-based data, to obtain column densities for various ion states of C, N, O, Na, Mg, Al, Si, P, S, Cl, Ar, Ca, Mn, Fe, and CO. N I and O I are fitted to a single empirical curve of growth with a velocity parameter (b) of 8 km/s; Mg II, Si II, P II, S II, Mn II, and Fe II are fitted to another curve with b between 3 and 9 km/s. Abundance determinations relative to H I show that: (1) C, N, P, S, and Ar are probably close to their solar values; (2) O may be depleted by about a factor of 2; (3) Mg, Al, Si, Cl, Mn, and Fe are depleted by a factor of 4 or more: (4) Al is depleted by at least a factor of 10 in the H II region; and (5) both N V and O VI are present, but not C IV. The N V/O VI ratio implies that the electron temperature in the H II region is about 275,000 K.

Morton, D. C.↗

ZTF SN Ia DR2: Improved SN Ia colors through expanded dimensionality with SALT3+

Context. Type Ia supernovae (SNe Ia) are a key probe in modern cosmology, as they can be used to measure luminosity distances at gigaparsec scales. Models of their light curves are used to project heterogeneous observed data onto a common basis for analysis. Aims. The SALT model currently used for SN Ia cosmology describes SNe as having two sources of variability, accounted for by a color parameter c , and a “stretch” parameter x 1 . We extend the model to include an additional parameter we label x 2 , to investigate the cosmological impact of currently unaddressed light-curve variability. Methods. We constructed a new SALT model, that we dub “SALT3+”. This model was trained by an improved version of the SALTshaker code, using training data combining a selection of the second data release of cosmological SNe Ia from the Zwicky Transient Facility and the existing SALT3 training compilation. Results. We find additional, coherent variability in supernova light curves beyond SALT3. Most of this variation can be described as phase-dependent variation in g − r and r − i color curves, correlated with a boost in the height of the secondary maximum in i -band. These behaviors correlate with spectral differences, particularly in line velocity. We find that fits with the existing SALT3 model tend to address this excess variation with the color parameter, leading to less informative measurements of supernova color. We find that neglecting the new parameter in light-curve fits leads to a trend in Hubble residuals with x 2 of 0.039 ± 0.005 mag, representing a potential systematic uncertainty. However, we find no evidence of a bias in current cosmological measurements. Conclusions. We conclude that extended SN Ia light-curve models promise mild improvement in the accuracy of color measurements, and corresponding cosmological precision. However, models with more parameters are unlikely to substantially affect current cosmological results.

Kenworthy, W. D. (ORCID:0000000251535983)↗

Integrating Light Curve and Atmospheric Modeling of Transiting Exoplanets

Spectral retrieval techniques are currently our best tool to interpret the observed exoplanet atmospheric data. Said techniques retrieve the optimal atmospheric components and parameters by identifying the best fit to an observed transmission/emission spectrum. Over the past decade, our understanding of remote worlds in our galaxy has flourished thanks to the use of increasingly sophisticated spectral retrieval techniques and the collective effort of the community working on exoplanet atmospheric models. A new generation of instruments in space and from the ground is expected to deliver higher quality data in the next decade; it is therefore paramount to upgrade current models and improve their reliability, their completeness, and the numerical speed with which they can be run. In this paper, we address the issue of reliability of the results provided by retrieval models in the presence of systematics of unknown origin. More specifically, we demonstrate that if we fit directly individual light curves at different wavelengths (L-retrieval), instead of fitting transit or eclipse depths, as it is currently done (S-retrieval), the said methodology is more sensitive against astrophysical and instrumental noise. This new approach is tested, in particular, when discrepant simulated observations from Hubble Space Telescope/Wide Field Camera 3 and Spitzer/IRAC are combined. We find that while S-retrievals converge to an incorrect solution without any warning, L-retrievals are able to flag potential discrepancies between the data sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Discrete Gust Model for Launch Vehicle Assessments

Analysis of spacecraft vehicle responses to atmospheric wind gusts during flight is important in the establishment of vehicle design structural requirements and operational capability. Typically, wind gust models can be either a spectral type determined by a random process having a wide range of wavelengths, or a discrete type having a single gust of predetermined magnitude and shape. Classical discrete models used by NASA during the Apollo and Space Shuttle Programs included a 9 m/sec quasi-square-wave gust with variable wavelength from 60 to 300 m. A later study derived discrete gust from a military specification (MIL-SPEC) document that used a "1-cosine" shape. The MIL-SPEC document contains a curve of non-dimensional gust magnitude as a function of non-dimensional gust half-wavelength based on the Dryden spectral model, but fails to list the equation necessary to reproduce the curve. Therefore, previous studies could only estimate a value of gust magnitude from the curve, or attempt to fit a function to it. This paper presents the development of the MIL-SPEC curve, and provides the necessary information to calculate discrete gust magnitudes as a function of both gust half-wavelength and the desired probability level of exceeding a specified gust magnitude.

Leahy, Frank B.↗

Tool For Making Curved Holes

Tool for use in electrical-discharge machining (EDM) guides EDM electrode in making curved holes. Guide rod fits in slot in arm, which moves through arc. Motion drives electrode into workpiece along desired curved path. Electrode burns into workpiece while arm rotates on spindle. Discharge cuts hole of same radius of curvature.

Allard, Robert↗

Optical variability of quasars with 20-yr photometric light curves

ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗

The influence of temperature on area-specific impedance and capacity of Li-ion cells with nickel-containing positive electrodes

The temperature-dependent behaviors of five nickel-containing positive electrodes (NCA, NMC811, NMC622, NMC532, and NMC111) in lithium-ion batteries are investigated using an electrochemical protocol involving rate studies, mild aging (~100 cycles), and hybrid pulse power characterization (HPPC). Tests are conducted using coin-cells with graphite negative electrodes at -20 °C, 0 °C, 20 °C, and 40 °C. Three techniques are compared for determining the area-specific impedance (ASI): i) fits to the rate study average voltages, ii) fitting to the entire voltage curves using a regularization scheme, and iii) HPPC. When fit to an Arrhenius-type equation, all methods yield similar apparent activation energies (±2 kJ/mol) for the impedance, which range from -20 to -31 kJ/mol for the electrodes. Impedance growth increases with temperature but remains at less than 0.2% per cycle for most electrodes and temperatures. NCA and NMC811 are the exceptions, which yield 0.5% and 1.5% increases in ASI per cycle, respectively, at 40 °C. For cells with the same electrodes, the capacities are similar at 20 and 40 °C but reduce at lower temperatures, with up to a 52% reduction at -20 °C and 2C. The fade in energy of the cells during C/3 cycling is attributed to decreasing capacity as opposed to increasing ASI.

25 ENERGY STORAGE↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗