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 739 records · Page 41

The magnetic anomaly of the Ivreazone

A magnetic field survey was made in the Ivreazone in 1969/70. The results were: significant anomaly of the vertical intensity is found. It follows the basic main part of the Ivrea-Verbano zone and continues to the south. The width of the anomaly is about 10 km, the maximum measures about +800 gamma. The model interpretation shows that possibly the anomaly belongs to an amphibolitic body, which in connection with the Ivrea-body was found by deep seismic sounding. Therefore, the magnetic anomaly provides further evidence for the conception that the Ivrea-body has to be regarded as a chip of earthmantle material pushed upward by tectonic processes.

Vertical Intensity of Magnetic Field↗

Itokawa's Opposition Surge seen by Hayabusa/AMICA

Using images acquired by the Hayabusa/AMICA instrument, along with Lederer et al.'s (2008) ground-based observations, we re-examine Itokawa's disk-integrated phase curve. The AMICA images provide critical opposition measurements (between 0.7deg - 9.3deg phase at 540 nm). Using Hapke's model (2012), we fit the updated phase curves at 5 different wavelengths. Preliminary modeling results show a range of porosity values commensurate with those in the literature (Ostro et al. 2004, Gundlach and Blum, 2012, Kiuchi and Nakamura 2014) based on an impact-generated grain size distribution function and grain size range evaluations from the AMICA data (Yano et al. 2006). This wide range on a global porosity is indicative of a highly variable porosity across the surface. The derived transport mean free path and the generally forward scattering nature of the global regolith are indicative of scattering centers (such as cracks, bubbles, and inclusions) that are small compared to the observational wavelengths. The derived regolith properties are compared with the imaging and sample analysis results, providing a test of the predictive capabilities of global disk-integrated measurements. This work suggests that the sub-pixel grain information could be extracted from the photometry, especially around opposition.

Dominque, D.↗

Mass loss in the 96 day binary UU Cancri

A series of 16 high-dispersion spectra at H-alpha have been obtained for the tidally distorted K giant-containing long-period binary UU Cnc, in order to both study the K giant's Doppler profiles and search for the effects of accretion onto the second component. While Doppler profiles of the system for a semidetached configuration fit the observations very well, those for existing overcontact light-curve solutions all yield poorer fits. Although the H-alpha line is always stronger than those of the common giants, its equivalent width is consistent with a K4 II classification reflective of the star's 50 solar radius size. Emission wings in H-alpha are possible evidence for an accretion disk.

Eaton, Joel A.↗

The Organic Aerosols of Titan's Atmosphere

One of Titan's many characteristics is the presence of a haze that veils its surface. This haze is composed of heavy organic particles and determining the chemical composition of these particles is a primary objective for future probes that would conduct in situ analysis. Meanwhile, solar occultations provide constraints on the optical characteristics of the haze layer. This paper describes solar occultation observations obtained by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft. These observations strongly constrain the optical characteristics of the haze layer. We detail the different steps involved in the processing of these data and apply them to two occultations that were observed at the South Pole and at the equator in order to investigate the latitudinal dependence of optical properties. The light curves obtained in seven atmospheric windows between 0.933-microns to 5-microns allow us to characterize atmospheric layers from 300 km to the surface. Very good fits of the light curves are obtained using a simple profile of number density of aerosols that is characterized by a scale height. The main difference between the South Pole and the equator is that the value of the scale height increases with altitude at the South Pole whereas it decreases at the equator. The vertically integrated amount of aerosols is similar at the two locations. The curve describing the cross-section versus wavelength is identical at the two locations suggesting that the aerosols have similar characteristics. Finally, we find that the two-way vertical transmission at 5-microns is as large as 80% at both locations.

solar occultation↗

A technique for constructing spectral reflectance curves from Viking lander camera multispectral data

A technique for evaluating the construction of spectral reflectance curves from multispectral data obtained with the Viking lander cameras is presented. The multispectral data is limited to 6 channels in the wave-length range 0.4 to 1.1 microns, and several of the channels suffer from appreciable out-of-band response. The technique represents the estimated reflectance curves as a linear combination of known basic functions with coefficients determined to minimize the error in the representation, and it permits all channels, with and without out-of-band response, to contribute equally valid information. The technique is evaluated for known spectral reflectance curves of 8 materials felt likely to be present on the Martian surface. The technique provides an essentially exact fit if the the reflectance curve has no pronounced maxima and minima. Even if the curve has pronounced maxima and minima, the fit is good and reveals the most dominant features. Since only 6 samples are available some short period features are lost. This loss is almost certainly due to undersampling rather than out-of-band channel response.

Park, S. K.↗

The late behavior of supernova 1987A. I - The light curve. II - Gamma-ray transparency of the ejecta

Observations of the late (t = 20-1500 days) bolometric light curve and the gamma-lines and X-rays from supernova 1987A are compared to theoretical models. It is found that 0.073 + or - 0.015 solar masses of freshly synthesized Ni-56 must be present to fit the bolometric light curve. The results place limits on the luminosity and presumed period of the newly formed pulsar/neutron star. In the second half of the paper, the problem of computing the luminosities in gamma-ray lines and in X-rays from supernova 1987A is addressed. High-energy observations suggest the development of large-scale clumping and bubbling of radioactive material in the ejecta. A model is proposed with a hydrogen envelope mass of about 7 solar masses, homologous scale expansion velocities of about 3000 km/s, and an approximately uniform mass distribution.

Arnett, W. David↗

The anomalous extinction curve in the direction of Rho Ophiuchi from 950 to 1180 A

We report the interstellar extinction from 950 to 1180 A in the direction of Rho Oph. These data were obtained with a moderate-resolution sounding rocket observation of Ro Oph, and two lightly reddened B stars of identical spectral type. Unlike Voyager observations in this wavelength region, our spectral resolution is sufficient to directly measure the contribution of molecular hydrogen along the line of sight and unambiguously identify the dust contribution. The derived curves compare favorably with the extrapolation of functional fits to the Rho Oph derived from longer wavelengths. In addition, we have fitted a model of dust extinction to the derived curve and find that the best fits indicate that the line of sight exhibits a larger than average grain size and a preferential concentration of graphite relative to silicates. The sharp rise seen at the shortest wavelengths (below 1000 A) cannot be explained by standard silicate-graphite dust models and may be indicative of a nonstandard grain shape or a missing component in the extinction models.

Green, James C.↗

The FU Orionis outburst as a thermal accretion event: Observational constraints for protostellar disk models

The results of the time-dependent disk models developed in Bell & Lin are compared with observed properties of FU Orionis variables. Specific models are fit to the light curves of Fu Ori, V1515 Cyg, and V1057 Cyg. The slow risetime of V1515 Cyg can be matched by a self-regulated outburst model. The rapid risetimes of FU Ori and V1057 Cyg can be fitted with the application of modest perturbations to the disk surface density. Model disks display spectral features characteristic of observed objects. The color evolution of V1057 Cyg is naturally explained if mass flux drops in the inner disk (r less than 1/4 AU) while remaining steady in the outer disk. The decrease in optical line width (rotational velocity) observed during the decay of V1057 Cyg may be accounted for by an outward-propagating ionization front. We predict that before final decay to the quiescent phase, short-wavelength line widths (lambda less than 1.5 microns) will again increase. It is suggested that FU Orionis outbursts primarily occur to systems during the embedded phase with ages less than several times 10(exp 5) yr.

Bell, K. R.↗

VW Hyi - The white dwarf revealed

Nonsimultaneous IUE, optical, and near-IR observations of VW Hyi at quiescence are presented. Using these and UV data from other investigations, a broad feature in the ultraviolet is identified with L-alpha absorption. The presence and width of the line imply that (1) the white dwarf in VW Hyi is directly visible in the UV and (2) the effective temperature of this star is approximately 18,000 + or - 2000 K for log g = 8. The continuum observations, combined with the J and K photometry of Sherrington et al., (1980), can be fit with a combination of this relatively cool white dwarf and a steady-state disk model with an accretion rate of 10 to the -11th solar masses/yr. Additional observations of the hump in the optical light curve can be reasonably fit by a 12,000-K blackbody. Such a source is consistent with the hump being a minor contribution to the system's overall continuum distribution shortward of 2000 A and longward of about 1 micron.

Mateo, M.↗

International planetary patrol observations of Saturn's rings. II - Four color phase curves and their analysis

New phase curves for Saturn's rings at an intermediate tilt angle B of about 17 deg are presented. Quantitative results for each of the A and B rings are reported in terms of the opposition effect, phase coefficient, and best logarithmic fit to the phase curve. There was no significant difference between the shape of the phase curves for the two rings in each of the four colors, and a four-parameter multiple scattering model of the rings was consistent with the observations. In this model, the difference in the phase curves for different colors can be explained by a variation in the single scattering albedo with wavelength. The observations allow the particles to have the same composition in the A and B rings, so that their different photometric behavior is explained by differences in optical depth and volume density in the two rings.

Esposito, L. W.↗

Photometric and Spectral Study of the Saturnian Satellites

Photometric and spectra analysis of data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) has yielded intriguing findings regarding the surface properties of several of the icy Saturnian satellites. Spectral cubes were obtained of these satellites with a wavelength distribution in the IR far more extensive than from any previous observations. Disk-integrated solar phase curves were constructed in several key IR wavelengths that are indicative of key properties of the surface of the body, such as macroscopic roughness, fluffiness (or the porosity of the surface), global albedo and scattering properties of surface particles. Polynomial fits to these phase curves indicate a linear albedo trend of the curvature of the phase functions. Rotational phase functions from Enceladus were found to exhibit a double-peaked sinusoidal curve, which shows larger amplitudes for bands corresponding to water ice and a linear amplitude-albedo trend. These functions indicate regions on the surface of the satellite of more recent geologic activity. In addition, recent images of Enceladus show tectonic features and an absence of impact craters on Southern latitudes which could be indicative of a younger surface. Investigations into the properties of these features using VIMS are underway.

Newman, Sarah F.↗

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗

Predicting Melt Curves of Energetic Materials Using Molecular Models

Here, in this work, the solid–liquid coexistence curves of classical fully flexible atomistic models of α-RDX and β-HMX were calculated using thermodynamically rigorous methodologies that identify where the free energy difference between the phases is zero. The free energy difference between each phase at a given state point was computed using the pseudosupercritical path (PSCP) method, and Gibbs–Helmholtz integration was used to evaluate the solid–liquid free energy difference as a function of temperature. This procedure was repeated for several pressures to determine points along the coexistence curve, which were then fit to the Simon–Glatzel functional form. While effective, this method is computationally expensive. An alternative approach is to compute the melting point at a single pressure via the PSCP method, and then use the Gibbs–Duhem integration technique to trace out the coexistence curve in a more computationally economical manner. Both approaches were used to determine the coexistence curve of α-RDX. The Gibbs–Duhem integration method was shown to generate a melt curve that is in good agreement with the PSCP-derived melt curve, while only costing ~10 % of the computational resources used for the PSCP method. For α-RDX, the predicted melting temperature increases significantly more for a given increase in pressure when compared to available experimental data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

An analysis of the shapes of ultraviolet extinction curves. II - The far-UV extinction

In this paper the properties of interstellar extinction in the far-ultraviolet region are examined utilizing IUE extinction curves for a primary data sample of 45 reddened Milky Way OB stars. These results are combined with those derived for the 2175 A bump in a previous study. It is found that IUE extinction curves can be represented by linear combinations of a Lorentzian-like 2175 A bump profile, a well-determined FUV curvature term, and an underlying linear component. The parameters of the linear component are strongly correlated and therefore only five free parameters are required to fit all of the curves in our sample. Three parameters describe the 2175 A bump, one parameter describes the linear background, and one parameter describes the strength of the FUV curvature term. The shape of the FUV curvature is found to be identical, to within the observational errors, for all the curves in the sample - which spans virtually the entire range of extinction curve morphologies observed in the Milky Way.

Fitzpatrick, Edward L.↗

Tidal effects in the O-type eclipsing system LY Aurigae

Existing photometry of the O9.5 III eclipsing binary LY Aur is discussed and combined into a visual mean light curve. It is found that the observed light curve is incompatible with previous solutions based on the Russell model, but a respectable fit to the observed light curve can be achieved by assuming von Zeipel gravity darkening provided the system is allowed to be in marginal contact, which involves an unrealistically large change in the mass ratio to obtain the appropriate eclipse depths and light ratio. The alternative possibilities are considered that the primary is tidally distorted more than expected for its relative radius and the observed mass ratio or that von Zeipel gravity darkening is not appropriate for O-type stars and photometric effects exist which magnify the effects of the primary's physical distortions in the light curve. A solution is favored in which the ratio of radii is near 0.75 and the primary is very highly gravity darkened relative to von Zeipel's law.

Eaton, J. A.↗

Insights into Space Solar Cell Durability Using SPICE Simulation Seeded by Current-Voltage Characteristics Parametrized Using the Lambert W Special Function

We developed and validated an automated routine for the fitting of I-V curve data to the single diode model according to an exact analytical solution. Our fitting routine was validated to show good noise immunity and high accuracy usingsimulated values from LTSPICE. We are thus able to automate parameter extraction from a dataset or arbitrary size. This parameterization allows for better simulation of performance of real cells and arrays and provides utility for a host of applications relevant to space arrays. We will use this methodology to determine the array performance of radiated cells over time, and simulate the performance of the arrays with bypass diodes and power electronics.

Solar Cell↗

A Likely Magnetic Activity Cycle for the Exoplanet Host M Dwarf GJ 3512

Current radial velocity data from specialized instruments contain a large amount of information that may pass unnoticed if their analysis is not accurate. The joint use of Bayesian inference tools and frequency analysis has been shown as effective in revealing exoplanets but they have been used less frequently to investigate stellar activity. We intend to use radial velocity data of the exoplanet host star GJ 3512 to investigate its magnetic activity. Our study includes the analysis of the photometric data available. The main objectives of our work are to constrain the orbital parameters of the exoplanets in the system, to determine the current level of activity of the star and to derive an activity cycle length for it. An adaptive importance sampling method was used to determine the parameters of the exoplanets orbit. Generalized Lomb–Scargle periodograms were constructed with both radial velocity curve and photometric data. A careful analysis of the harmonic frequencies was conducted in each periodogram. Our fit to multiple Keplerian orbits constrained the orbital parameters of two giant gas planets orbiting the star GJ 3512. The host star showed an increase of its magnetic activity during the last observing campaign. The accurate fit of the radial velocity curve data to the multi-Keplerian orbit permitted to reveal the star rotation in the residuals of the best fit and estimate an activity cycle length of ~14 yr.

79 ASTRONOMY AND ASTROPHYSICS↗