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 1,099 records · Page 61

THE 2011 JUNE 23 STELLAR OCCULTATION BY PLUTO: AIRBORNE AND GROUND OBSERVATIONS

On 2011 June 23, stellar occultations by both Pluto (this work) and Charon (future analysis) were observed from numerous ground stations as well as the Stratospheric Observatory for Infrared Astronomy (SOFIA). This first airborne occultation observation since 1995 with the Kuiper Airborne Observatory resulted in the best occultation chords recorded for the event, in three visible wavelength bands. The data obtained from SOFIA are combined with chords obtained from the ground at the IRTF, the U.S. Naval Observatory Flagstaff Station, and Leeward Community College to give the detailed state of the Pluto-Charon system at the time of the event with a focus on Pluto's atmosphere. The data show a return to the distinct upper and lower atmospheric regions with a knee or kink in the light curve separating them as was observed in 1988, rather than the smoothly transitioning bowl-shaped light curves of recent years. The upper atmosphere is analyzed by fitting a model to all of the light curves, resulting in a half-light radius of 1288 {+-} 1 km. The lower atmosphere is analyzed using two different methods to provide results under the differing assumptions of particulate haze and a strong thermal gradient as causes for the lower atmospheric diminution of flux. These results are compared with those from past occultations to provide a picture of Pluto's evolving atmosphere. Regardless of which lower atmospheric structure is assumed, results indicate that this part of the atmosphere evolves on short timescales with results changing the light curve structures between 1988 and 2006, and then reverting these changes in 2011 though at significantly higher pressures. Throughout these changes, the upper atmosphere remains remarkably stable in structure, again except for the overall pressure changes. No evidence of onset of atmospheric collapse predicted by frost migration models is seen, and the atmosphere appears to be remaining at a stable pressure level, suggesting it should persist at this full level through New Horizon's flyby in 2015.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust in the small Magellanic Cloud. 2: Dust models from interstellar polarization and extinction data

We study the dust in the Small Magellanic Cloud using our polarization and extinction data (Paper 1) and existing dust models. The data suggest that the monotonic SMC extinction curve is related to values of lambda(sub max), the wavelength of maximum polarization, which are on the average smaller than the mean for the Galaxy. On the other hand, AZV 456, a star with an extinction similar to that for the Galaxy, shows a value of lambda(sub max) similar to the mean for the Galaxy. We discuss simultaneous dust model fits to extinction and polarization. Fits to the wavelength dependent polarization data are possible for stars with small lambda(sub max). In general, they imply dust size distributions which are narrower and have smaller mean sizes compared to typical size distributions for the Galaxy. However, stars with lambda(sub max) close to the Galactic norm, which also have a narrower polarization curve, cannot be fit adequately. This holds true for all of the dust models considered. The best fits to the extinction curves are obtained with a power law size distribution by assuming that the cylindrical and spherical silicate grains have a volume distribution which is continuous from the smaller spheres to the larger cylinders. The size distribution for the cylinders is taken from the fit to the polarization. The 'typical', monotonic SMC extinction curve can be fit well with graphite and silicate grains if a small fraction of the SMC carbon is locked up in the grain. However, amorphous carbon and silicate grains also fit the data well. AZV456, which has an extinction curve similar to that for the Galaxy, has a UV bump which is too blue to be fit by spherical graphite grains.

Rodrigues, C. V.↗

DISCOVERY OF GAMMA-RAY PULSATIONS FROM THE TRANSITIONAL REDBACK PSR J1227-4853

The 1.69 ms spin period of PSR J1227−4853 was recently discovered in radio observations of the low-mass X-ray binary XSS J12270−4859 following the announcement of a possible transition to a rotation-powered millisecond pulsar state, inferred from decreases in optical, X-ray, and gamma-ray flux from the source. We report the detection of significant (5σ) gamma-ray pulsations after the transition, at the known spin period, using ∼1 year of data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. The gamma-ray light curve of PSR J1227−4853 can be fit by one broad peak, which occurs at nearly the same phase as the main peak in the 1.4 GHz radio profile. The partial alignment of light-curve peaks in different wavebands suggests that at least some of the radio emission may originate at high altitude in the pulsar magnetosphere, in extended regions co-located with the gamma-ray emission site. We folded the LAT data at the orbital period, both pre- and post-transition, but find no evidence for significant modulation of the gamma-ray flux. Analysis of the gamma-ray flux over the mission suggests an approximate transition time of 2012 November 30. Continued study of the pulsed emission and monitoring of PSR J1227−4853, and other known redback systems, for subsequent flux changes will increase our knowledge of the pulsar emission mechanism and transitioning systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Lunar occultation of Saturn. I - The diameters of Tethys, Dione, Rhea, Titan, and Iapetus

The diameters of Tethys, Dione, Rhea, Titan and Iapetus were determined from observations of their March 30, 1974, lunar occultations, made with the Mauna Kea 224 and 61 cm telescopes. Light curves were obtained simultaneously in four colors, and the difference between the time of occultation at the two telescopes provided a direct measurement of the slope of the lunar limb, found to be small in all cases. The satellite diameters were determined by least-squares fits of model occultation light curves to the data. Titan's diameter must be at least 5800 km, much larger than the currently accepted value of 5000 km, making it the largest satellite in the solar system. This larger diameter implies a low mean density. For the other four satellites arguments are presented in favor of accepting the occultation diameters corresponding to limb darkened disks. Except for Titan, the lunar occultation diameters generally agree with previous diskmeter and radiometric determinations.

Elliot, J. L.↗

A matched-filter approach to radio variability and transients: searching for orphan afterglows in the VAST Pilot Survey

Radio transient searches using traditional variability metrics struggle to recover sources whose evolution time-scale is significantly longer than the survey cadence. Motivated by the recent observations of slowly evolving radio afterglows at gigahertz frequency, we present the results of a search for radio variables and transients using an alternative matched-filter approach. We designed our matched-filter to recover sources with radio light curves that have a high-significance fit to power-law and smoothly broken power-law functions; light curves following these functions are characteristic of synchrotron transients, including ‘orphan’ gamma-ray burst afterglows, which were the primary targets of our search. Applying this matched-filter approach to data from Variables and Slow Transients Pilot Survey conducted using the Australian SKA Pathfinder, we produced five candidates in our search. Subsequent Australia Telescope Compact Array observations and analysis revealed that: one is likely a synchrotron transient; one is likely a flaring active galactic nucleus, exhibiting a flat-to-steep spectral transition over 4 months; one is associated with a starburst galaxy, with the radio emission originating from either star formation or an underlying slowly evolving transient; and the remaining two are likely extrinsic variables caused by interstellar scintillation. The synchrotron transient, VAST J175036.1–181454, has a multifrequency light curve, peak spectral luminosity, and volumetric rate that is consistent with both an off-axis afterglow and an off-axis tidal disruption event; interpreted as an off-axis afterglow would imply an average inverse beaming factor $\langle f^{-1}_{\text{b}} \rangle = 860^{+1980}_{-710}$, or equivalently, an average jet opening angle of $\langle \theta _{\textrm {j}} \rangle = 3^{+4}_{-1}\,$ deg.

79 ASTRONOMY AND ASTROPHYSICS↗

Curved Finite Elements and Curve Approximation

The approximation of parameterized curves by segments of parabolas that pass through the endpoints of each curve segment arises naturally in all quadratic isoparametric transformations. While not as popular as cubics in curve design problems, the use of parabolas allows the introduction of a geometric measure of the discrepancy between given and approximating curves. The free parameters of the parabola may be used to optimize the fit, and constraints that prevent overspill and curve degeneracy are introduced. This leads to a constrained optimization problem in two varibles that can be solved quickly and reliably by a simple method that takes advantage of the special structure of the problem. For applications in the field of computer-aided design, the given curves are often cubic polynomials, and the coefficient may be calculated in closed form in terms of polynomial coefficients by using a symbolic machine language so that families of curves can be approximated with no further integration. For general curves, numerical quadrature may be used, as in the implementation where the Romberg quadrature is applied. The coefficient functions C sub 1 (gamma) and C sub 2 (gamma) are expanded as polynomials in gamma, so that for given A(s) and B(s) the integrations need only be done once. The method was used to find optimal constrained parabolic approximation to a wide variety of given curves.

Baart, M. L.↗

PSR J2030+364I: Radio Discovery and Gamma-ray Study of a Middle-aged Pulsar in the Now Identified Fermi-LAT Source 1FGL J2030.0+3641

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi-LAT sources, we have discovered the middle-aged pulsar J2030+3641, associated with IFGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the entire Fermi mission. With a rotation period of 0.28, spin-down luminosity of 3 x 10(exp 34) erg/s, and characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1 % that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM = 246 pc/cu cm. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive - PSR J2030+364 I would have been found blindly in gamma rays if only > or approx. 0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

Camilo, F.↗

A new form of the calibration curve in radiochromic dosimetry. Properties and results

Purpose: This work describes a new form of the calibration curve for radiochromic dosimetry that depends on one fit parameter. Some results are reported to show that the new curve performs as well as those previously used and, more importantly, significantly reduces the dependence on the lot of films, the film orientation on the scanner, and the time after exposure. Methods: The form of the response curve makes use of the net optical densities ratio against the dose and has been studied by means of the Beer–Lambert law and a simple modeling of the film. The new calibration curve has been applied to EBT3 films exposed at 6 and 15 MV energy beams of linear accelerators and read-out in transmission mode by means of a flatbed color scanner. Its performance has been compared to that of two established forms of the calibration curve, which use the optical density and the net optical density against the dose. Four series of measurements with four lots of EBT3 films were used to evaluate the precision, accuracy, and dependence on the time after exposure, orientation on the scanner and lot of films. Results: The new calibration curve is roughly subject to the same dose uncertainty, about 2% (1 standard deviation), and has the same accuracy, about 1.5% (dose values between 50 and 450 cGy), as the other calibration curves when films of the same lot are used. Moreover, the new calibration curve, albeit obtained from only one lot of film, shows a good agreement with experimental data from all other lots of EBT3 films used, with an accuracy of about 2% and a relative dose precision of 2.4% (1 standard deviation). The agreement also holds for changes of the film orientation and of the time after exposure. Conclusions: The dose accuracy of this new form of the calibration curve is always equal to or better than those obtained from the two types of curves previously used. The use of the net optical densities ratio considerably reduces the dependence on the lot of films, the landscape/portrait orientation, and the time after exposure. This form of the calibration curve could become even more useful with new optical digital devices using monochromatic light.

60 APPLIED LIFE SCIENCES↗

Lunar occultation of Saturn. III - How big is Iapetus

The orbital light curve of Iapetus is analyzed, and the bright material on the trailing face of the satellite is found to have an effective limb-darkening parameter of 0.62 (+0.10, -0.12). This result is combined with observations of a lunar occultation of Iapetus to derive a radius of 718 (+87, -78) km for the satellite. The specific patchy albedo model of Morrison et al. (1975) is shown to fit the observed orbital light curve precisely and to yield an occultation radius of 724 + or - 60 km. It is noted that the values determined for the limb-darkening parameter and for the normal reflectance (in V) of the material at the center of the bright face are consistent with a partial water-frost cover on the bright side of Iapetus. It is concluded that if the albedo model of Morrison et al. reliably represents the surface of Iapetus, the radius of that satellite is 724 + or - 60 km, which is consistent with the radiometric radius of 835 (+50, -75) km derived by Morrison et al. (1975).

Veverka, J.↗

Direct Measurements of Interplanetary Dust Particles in the Vicinity of Earth

The direct measurements made by the Explorer VIII satellite provide the first sound basis for analyzing all available direct measurements of the distribution of interplanetary dust particles. The model average distribution curve established by such an analysis departs significantly from that predicted by the (uncertain) extrapolation of results from meteor observations. A consequence of this difference is that the daily accretion of interplanetary particulate matter by the earth is now considered to be mainly dust particles of the direct measurements range of particle size. Almost all the available direct measurements obtained with microphone systems on rockets, satellites, and spacecraft fit directly on the distribution curve defined by Explorer VIII data. The lack of reliable datum points departing significantly from the model average distribution curve means that available direct measurements show no discernible evidence of an appreciable geocentric concentration of interplanetary dust particles.

McCracken, C. W.↗

Time-dependent behavioral data from zebrafish reveals novel signatures of chemical toxicity using point of departure analysis

High-content imaging of embryonic zebrafish behavior can be used as a screening approach to rapidly evaluate the relative potential for chemicals to cause toxicity. However, most statistical methods applied to these data so far transform movement values to incidence-based “hits” and calculate lowest effect levels (LELs), which loses individual fish resolution of behavior and defies hazard ranking due to reliance on applied dose levels. Here, we have developed and parallelized a workflow to calculate benchmark dose (BMD) values from dynamic, high-content zebrafish behavior data that scales for high-throughput chemical screening. To appropriately capture the dynamics of zebrafish response from light to dark stimulus, we summarized time-dependent data using both area under the curve and immediate change at the transition point into two novel metrics that characterized abnormal behavior as a function of chemical concentration. The BMD workflow was applied to calculate BMD 10 values of 1,060 ToxCast chemicals for 24 zebrafish endpoints, including behavior, mortality and morphology. The BMD 10 values provided better precision and separation than LELs for clustering chemicals since they were derived from models that best-fit their concentration-response curves. Clustering analysis revealed behavioral signatures as the most sensitive endpoints. High concordance in chemical activity was observed between ToxCast assays and zebrafish behavioral endpoints, however ToxPi analysis indicated that rankings based on in vitro data was not a reliable predictor of in vivo rankings for lower potency chemicals. The methods developed in this work will enable the use of high-content zebrafish behavioral screening data for BMD analysis in toxicological hazard assessment.

59 BASIC BIOLOGICAL SCIENCES↗

Influence of living grass Roots and endophytic fungal hyphae on soil hydraulic properties

Soil hydraulic properties are often estimated based on laboratory data or pedotransfer functions dependent on soil physical properties, which often do not consider potential impacts of soil roots or fungal hyphae. Here, we first review current knowledge of how these soil biotic components affect hydraulic properties, then we conducted laboratory experiments to specifically test if the presence of roots and mycorrhizal fungi had a significant effect on the hydraulic properties of two soils with contrasting textures: Flint sand and Hamblen silt loam. Soil cores were seeded with (Panicum virgatum) and grown in a greenhouse over three separate growth periods. The endophytic fungus Serendipita indica was injected as liquid inoculant into designated mycorrhizal cores. Saturated hydraulic conductivity (K sat ) measurements were made with a constant head permeameter, and soil water retention curves were obtained by the evaporation method, supplemented at the dry end for Hamblen silt loam with water activity meter data. Retention curve parameters were obtained by fitting the van Genuchten equation to the resulting measurements. Mean root volume ratios were higher in the mycorrhizal inoculated treatment than in the uninoculated treatment for both soils. For Flint sand, analysis of variance revealed that K sat was reduced by the presence of roots as compared to bare soil. This was likely due to roots clogging soil pores. Results also indicated the presence of roots changed the shape of the water retention curve for Flint sand by increasing water content at saturation and by reducing the slope of the curve. These changes suggested roots created additional porosity and broadened the pore-size distribution. The presence of mycorrhizal fungi accentuated the root effects. The influence of roots and mycorrhizal fungi on hydraulic properties was less obvious for the Hamblen silt loam, as none of the treatments differed from each other at p < 0.05. The results highlight the necessity to consider the impact of root and fungal structures on models of soil hydraulic properties.

59 BASIC BIOLOGICAL SCIENCES↗

Fantômas unconfined: global QCD fits with Bézier parameterizations

Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.

Bézier curves↗

Investigating the Effect of Cosmic Opacity on Standard Candles

Standard candles can probe the evolution of dark energy over a large redshift range. But the cosmic opacity can degrade the quality of standard candles. In this paper, we use the latest observations, including Type Ia supernovae (SNe Ia) from the “joint light-curve analysis” sample and Hubble parameters, to probe the opacity of the universe. A joint fitting of the SNe Ia light-curve parameters, cosmological parameters, and opacity is used in order to avoid the cosmological dependence of SNe Ia luminosity distances. The latest gamma-ray bursts are used in order to explore the cosmic opacity at high redshifts. The cosmic reionization process is considered at high redshifts. We find that the sample supports an almost transparent universe for flat ΛCDM and XCDM models. Meanwhile, free electrons deplete photons from standard candles through (inverse) Compton scattering, which is known as an important component of opacity. This Compton dimming may play an important role in future supernova surveys. From analysis, we find that about a few per cent of the cosmic opacity is caused by Compton dimming in the two models, which can be corrected.

79 ASTRONOMY AND ASTROPHYSICS↗

Young stellar object variability (YSOVAR): Long timescale variations in the mid-infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 and 4.5 μm) time series photometry of the Orion Nebula Cluster plus smaller footprints in 11 other star-forming cores (AFGL 490, NGC 1333, Mon R2, GGD 12-15, NGC 2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC 1396A, and Ceph C). There are ∼29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the 'standard sample' on which we calculate statistics, consisting of fast cadence data, with epochs roughly twice per day for ∼40 days. We also define a 'standard sample of members' consisting of all the IR-selected members and X-ray-selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data—the Stetson index, a χ{sup 2} fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of six to seven years by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data set; out of members and field objects combined, at most 0.02% may have transient IR excesses.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-Ray Pulsar Light Curves in Offset Polar Cap Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres, used to model high-energy light curves have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profile~. We find that. corn pared to the profile:-; derived from :-;ymmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines. increases significantly relative to the off-peak emission. formed along leading field lines. The enhanced contrast produces greatly improved slot gap model fits to Fermi pulsar light curves like Vela, which show very little off-peak emIssIon.

Harding, Alice K.↗

Probability Density Functions of Observed Rainfall in Montana

The question of whether a rain rate probability density function (PDF) can vary uniformly between precipitation events is examined. Image analysis on large samples of radar echoes is possible because of advances in technology. The data provided by such an analysis easily allow development of radar reflectivity factors (and by extension rain rate) distribution. Finding a PDF becomes a matter of finding a function that describes the curve approximating the resulting distributions. Ideally, one PDF would exist for all cases; or many PDF's that have the same functional form with only systematic variations in parameters (such as size or shape) exist. Satisfying either of theses cases will, validate the theoretical basis of the Area Time Integral (ATI). Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89 percent of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit. Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89% of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit.

Larsen, Scott D.↗

SU-F-T-170: Patient Surface Dose Measurements Using Optically Stimulated Luminescence Dosimeters in Scanning Proton Beam Therapy for Breast Cancer

Purpose: To establish patient surface dose dosimetry for scanning proton beam therapy (SPBT) for breast cancer using optically stimulated luminescence dosimeters (OSLD). Methods: OSLDs were calibrated with SPB under the similar conditions as the treatments for breast cancer. A range shifter (RS) of 5 cm water equivalent thickness (WET) was used. The air gap from the surface of the range shifter to the surface of the phantom was 15 cm. A uniform planar dose generated by nominal energy of 118 MeV was delivered. The range of 118 MeV proton beam after the 5cm RS is approximately 5 cm in water, which is the common range for breast treatments. The OSLDs were placed on the surface of high density polyethylene slabs, and a bolus of 1.06 cm WET was used for buildup. A variety of dose levels in the range of 0.5 to 8 Gy were delivered. Under the same condition, an ADCL calibrated parallel plate (PP) chamber was used to measure the reference dose. The correlation between the output signals of OSLDs and the reference doses was established. The calibration of OSLD was verified against the PP chamber measurements for two SPBT breast plans calculated for two patients. Results: the least squares fitting for the OSLD calibration curve was a polynomial function to the order of 2 in the range of 0.5 to 8 Gy (RBE). The differences between the dose measured with OSLDs and PP chamber were within 3% for the two breast proton plans. Conclusion: the calibrated OSLDs under the similar conditions as the treatments can be used for patient surface dose measurements.

60 APPLIED LIFE SCIENCES↗