Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 217 records · Page 12

Kepler Data Validation I: Architecture, Diagnostic Tests, and Data Products for Vetting Transiting Planet Candidates

The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for approx. 200,000 individual stellar targets in at least one observing quarter in the primary Kepler Mission. Light curves for target stars are extracted in the Kepler Science Data Processing Pipeline, and are searched for transiting planet signatures. A Threshold Crossing Event is generated in the transit search for targets where the transit detection threshold is exceeded and transit consistency checks are satisfied. These targets are subjected to further scrutiny in the Data Validation (DV) component of the Pipeline. Transiting planet candidates are characterized in DV, and light curves are searched for additional planets after transit signatures are modeled and removed. A suite of diagnostic tests is performed on all candidates to aid in discrimination between genuine transiting planets and instrumental or astrophysical false positives. Data products are generated per target and planet candidate to document and display transiting planet model fit and diagnostic test results. These products are exported to the Exoplanet Archive at the NASA Exoplanet Science Institute, and are available to the community. We describe the DV architecture and diagnostic tests, and provide a brief overview of the data products. Transiting planet modeling and the search for multiple planets on individual targets are described in a companion paper. The final revision of the Kepler Pipeline code base is available to the general public through GitHub. The Kepler Pipeline has also been modified to support the Transiting Exoplanet Survey Satellite (TESS) Mission which is expected to commence in 2018.

data analysis↗

Modeling the structure of H I in the galactic disk

For comparison with a detailed H I survey, we calculated models for the H I intensity distribution, T(B)(l, b, v), for a 18 x 12 deg region in the disk of the Galaxy toward l = 205 deg. In one class of models, neutral hydrogen is in the form of clouds, while in the other class of models it is a continuous medium with holes whose sizes are typical of supernova remnants and superbubbles. The cloud models are a poor representation of the data in that the amplitude fluctuations are dissimilar and the continuity of the H I observations is not reproduced. A continuous H I medium with holes provides a much better match to the data, and the best-fit models suggest that the filling factor for the H I is 20-50 percent.

Bregman, Joel N.↗

The Solar Chromosphere/Corona Interface: FUV-EUV Observations and Modeling of Unresolved Coronal Funnels - I

A coronal funnel model, developed by Rabin (199l), was tested against a calibrated spectroheliogram recorded in 171 - 175 Angstrom bandpass. This image was recorded on board a sounding rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array, II (MSSTA II), MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center and the Lawrence Livermore National Laboratory, is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow bandpass, high resolution images (1 - 3 arc sec) at selected FUV, EUV and soft X-ray wavelengths (44 Angstroms - 1550 Angstroms). Using full-disk images centered at 1550 Angstroms (C IV) and 173 Angstroms (FE IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found, was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171 - 175 Angstrom bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the Chianti spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of FE IX plasma. The funnels are found to have constrictions, Gamma approx. 6 - 20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely, that a global network of funnels must cover approx. 70 - 95 % of the total solar surface area seen in the 171 - 175 Angstrom image, and = 45 % of the disk area seen in the 1550 Angstrom bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et af. (1974) and put forth in more detail by Gabriel (1976). Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by Raymond & Doyle (1981) for temperatures, 250,000 K = T = 650,000 K, based on full-disk observations made on board by SkyLab.

Martinez-Galarce, D. S.↗

The Solar Chromosphere/Corona Interface. I. Far-Ultraviolet to Extreme-Ultraviolet Observations and Modeling of Unresolved Coronal Funnels

A coronal funnel model, developed by D. Rabin, was tested against a calibrated spectroheliogram recorded in the 170-1 75 A bandpass. This image was recorded on board a sounding-rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array II (MSSTA II). MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center, and the Lawrence Livermore National Laboratory' is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow-bandpass, high-resolution images (1 sec.- 3 sec.) at selected far-ultraviolet (FUV), extreme-ultraviolet (EUV), and soft X-ray wavelengths (44-1550 A). Using full disk images centered at 1550 A (C IV) and 173 A (Fe IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171-175 A bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the CHIANTI spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of Fe IX plasma. The funnels are found to have constrictions, Gamma approx. 6-20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely,that a global network of funnels must cover approx. 700 - 95% of the total solar surface area seen in the 171-175 A image, and a 5% of the disk area seen in the 1550 A bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et al. and put forth in more detail by Gabriel. Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by J. C. Raymond & J. G. Doyle for temperatures 250,000 K less than or = T less than or = 650,000 K, based on full-disk observations made on board Skylab.

Martinez-Galarce, Dennis S.↗

The Near-Earth Encounter of Asteroid 308635 (2005 YU55): Thermal IR Observations

The near-Earth approach (0.00217 AU, or 0.845 lunar distances) of the C-type asteroid 308635 (2005 YU55) in November 2011 presented a rare opportunity for detailed observations of a low-albedo NEA in this size range. As part of a multi-telescope campaign to measure visible and infrared spectra and photometry, we obtained mid-infrared (approx. 8 to 22 micron) photometry and spectroscopy of 2005 YU55 using Michelle on the Gemini North telescope on UT November 9 and 10,2011. An extensive radar campaign together with optical light-curves established the rotation state of YU55. In addition, the radar imaging resulted in a shape model for the asteroid, detection of numerous boulders on its surface, and a preliminary estimate of its equatorial diameter at 380 +/- 20 m. In a preliminary analysis, applying the radar and lightcurve-derived parameters to a rough-surface thermophysical model fit to the Gemini/Michelle thermal emission photometry results in a thermal inertia range of approximately 500 to 1500 J/sq m/0.5s/K, with the low-thermal-inertia solution corresponding to the small end of the radar size range and vice versa. Updates to these results will be presented and modeling of the thermal contribution to the measured near-infrared spectra from Palomar/Triplespec and IRTF/SpeX will also be discussed.

Lim, Lucy F.↗

Architectures of Exoplanetary Systems: Towards a Multi-planet Model for Reproducing the Planet Radius Valley and Intra-system Size Similarity

Multi-planet systems provide crucial insights into the architectures and correlations within planetary systems, which in turn offer clues into their formation and evolution histories. Through the thousands of exoplanet candidates that it discovered, NASA's Kepler mission revealed that systems with multiple transiting planets are common. We describe a new model for inferring the exoplanet population observed by Kepler, which combines a semi-parametric description of planetary systems with a prescription for envelope mass-loss due to photoevaporation. This new model is a"hybrid" of the maximum AMD (angular momentum deficit) model in which the orbital excitations of the planets are drawn by distributing the system's critical AMD (He et al. 2020), and a photoevaporation model for sculpting the observed radius valley (Neil & Rogers 2020). This is an update to the SysSim suite of exoplanet population models derived from forward modeling the Kepler detection pipeline. The hybrid model is constrained by the observed statistics of multi-transiting systems, including additional constraints from the distribution of planet radii and their intra-system size patterns. This model is capable of producing a planet radius valley as well as a preference for planet size similarity. However, we find that the Kepler catalog exhibits features that are even more significant (i.e., it has a deeper observed radius valley and stronger size similarity) than the best-fit model. We use this model to infer the intrinsic planet radius valley as well as the distribution of primordial planet masses and core masses.

kepler↗

Silicon vidicon imaging of Jupiter - 4100- to 8300-A absolute reflectivities and limb darkening of spatially resolved regions

Jupiter was observed in six continuum wavelength channels in the region 4100-8300 A, using a silicon vidicon imaging photometer. Spectral reflectivities and high spatial resolution limb-darkening curves for several belts and zones have been extracted from the data. Simple model fits to the data yield information regarding spectral and spatial variations in single-scattering albedos and shape of particle single-scattering phase functions. Belts appear to be more backscattering than zones, particularly in the blue. The data are in moderate agreement with limb-darkening predicted by models derived from the center-to-limb variation in equivalent width of the H2 4-0 S(1) quadrupole line (Cochran, 1976) in the South Tropical Zone, but strongly disagree with the results of such models for the North Equatorial Belt.

Diner, D. J.↗

Mapping tall shrub biomass in Alaska at landscape scale using structure-from-motion photogrammetry and lidar

Warming in arctic and boreal regions is increasing shrub cover and biomass. In southcentral Alaska, willow (Salix spp.) and alder (Alnus spp.) shrubs grow taller than many tree species and account for a substantial proportion of aboveground biomass, yet they are not individually measured as part of the operational Forest Inventory and Analysis (FIA) Program. The goal of this research was to test methods for landscape-scale mapping of tall shrub biomass in upper montane and subalpine environments using FIA-type plot measurements (n = 51) and predictor variables from imagery-based structure-from-motion (SfM) and airborne lidar. Specifically, we compared biomass models constructed from imagery acquired by unmanned aerial vehicle (UAV; ~1.7 cm pixels), imagery from the NASA Goddard's Lidar, Hyperspectral, and Thermal Airborne Imager (G-LiHT; ~3.1 cm pixels), and concomitant G-LiHT small-footprint lidar. Tall shrub biomass was most accurately predicted at 5 m resolution (R^2 = 0.81, RMSE = 1.09 kg/sq. m) using G-LiHT SfM color and structure variables. Lidar-only models had lower precision (R^2 = 0.74, RMSE = 1.26 kg/sq. m), possibly due to reduced model information content from variable multicollinearity or lower data density. Separate models for upper montane zones with trees and shrubs and subalpine zones with only shrubs were always chosen over single models based on minimization of Akaike's Information Criterion, indicating the need for variable sets robust to overhanging tree canopy. Decreasing point density from UAV (5000–8000 pts./sq. m) to the G-LiHT SfM point cloud (500–2000 pts./sq. m) had little impact on model fit, suggesting that high-resolution airborne imagery can extend SfM approaches well beyond line-of-sight restrictions for UAV platforms. Overall, our results confirmed that SfM from high-resolution imagery is a viable approach to estimate shrub biomass in the boreal region, especially when an existing lidar terrain model and local field calibration data are available to quantify uncertainty in the SfM point cloud and landscape-scale estimates of shrub biomass.

Michael Alonzo↗

A moving-window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving-window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm, tests for parameter validity, and the on-line design process.

Montgomery, R. C.↗

A moving window parameter adaptive control system for the F8-DFBW aircraft

This paper describes the moving window parameter adaptive control system developed for the NASA F8-DFBW aircraft. The control system employs a parameter identification process that, iteratively, adjusts parameters of a model of the aircraft motions in a batch-processing manner so that responses generated from the model fit the outputs of sensors stored in a finite record referred to as the moving window. Tests are made on the validity of the parameter estimates before using the parameters in an on-line design process. The on-line design process is an algebraic mapping of the parameters of the model into primary control system feedback and feedforward gains. The mapping was selected to satisfy specific flying quality characteristics over the range of parameter variations expected. Results are presented from simulation studies on the identification algorithm made during the development of the system.

Dunn, H. J.↗

Evaluation of a canopy reflectance model for LAI estimation through its inversion

A technique for estimating leaf area indices (LAIs) based on bidirectional canopy reflectance (CR) data is applied to three plant canopies: a naturally growing healthy soybean canopy; and a clumped and tufted orchardgrass canopy, respectively. The CR data were collected using a PARABOLA instrument which is capable of acquiring complete sky-and-ground looking hemispheres in 11 seconds. The model fit and LAI estimates were good for the soybean and clumped orchardgrass canopies, but poor for the tufted orchardgrass canopy when the maximum zenith angle was less than 50 percent. It is shown that the biophysical parameter estimation based on CR measurements applied well to homogeneous herbaceous vegetation types, while better CR models are needed to adequately represent discontinuous plant canopies.

Goel, N. S.↗

An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Subsolar Water Abundance

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 μm, in addition to a Spitzer 4.5 μm secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free-retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally ≲ 0.1× solar (0.001–0.7× solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical-equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they, too, prefer low metallicities ([M/H] ≲ −2, consistent with the free-retrieval results). However, all of the retrievals should be interpreted with some caution because they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find that the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

Knicole D. Colón↗

Resolving topographic detail on Venus by modeling complex Magellan altimetry echoes

Magellan's altimeter is providing some of the finest resolution topography of Venus achieved to date. Nevertheless, efforts continue to improve the topographic resolution whenever possible. One effort to this end is stereoscopic imaging, which provides topography at scales similar to that of the synthetic aperture radar (SAR). However, this technique requires two SAR images of the same site to be obtained and limits the utility of this method. In this paper, we present another method to resolve topographic features at scales smaller than that of an altimeter footprint, which is more globally applicable than the stereoscopic approach. Each pulse which is transmitted by Magellan's altimeter scatters from the planet and echoes to the receiver, delayed based on the distance between the spacecraft and each surface element. As resolved in time, each element of an altimetry echo represents the sum of all points on the surface which are equidistant from the spacecraft. Thus, individual returns, as a function of time, create an echo profile which may be used to derive properties of the surface, such as the scattering law or, in this case, the topography within the footprint. The Magellan project has derived some of this information by fitting model templates to radar echo profiles. The templates are calculated based on Hagfor's Law, which assumes a smooth, gently undulating surface. In most regions these templates provide a reasonable fit to the observed echo profile; however, in some cases the surface departs from these simple assumptions and more complex profiles are observed. Specifically, we note that sub-footprint topographic relief apparently has a strong effect on the shape of the echo profile. To demonstrate the effects of sub-resolution relief on echo profiles, we have calculated the echo shapes from a wide range of simple topographic models. At this point, our topographic models have emphasized surfaces where only two dominant elevations are contained within a footprint, such as graben, ridges, crater rims, and central features in impact craters.

Lovell, Amy J.↗

The Suzaku Observation of NGC 3516: Complex Absorption and the Broad and Narrow Fe K Lines

We present results from a 150 ksec Suzaku observation of the Seyfert 1 NGC 3516 in October 2005. The source was in a relatively highly absorbed state. Our best-fit model is consistent with partial covering by a lowly-ionized absorber with a column density near 5x10(exp 22) cm(exp -2) and with a covering fraction 96-100 percent. Narrow K-shell absorption features due to He- and H-like Fe confirm the presence of a high-ionization absorbing component as well. A broad Fe K(alpha) diskline is required in all fits, even after the complex absorption is taken into account; an additional partial-covering component is an inadequate substitute for the continuum curvature associated with the broad line. The narrow Fe Ka line at 6.4 keV is resolved, yielding a velocity width commensurate with the optical Broad Line Region. The strength of the Compton reflection hump suggests a contribution mainly from the broad Fe line origin. We include in our model soft band emission lines from He- and H-like ions and radiative recombination lines, consistent with photo-ionization, though a small contribution from collisional ionization is possible.

Markowitz, Alex↗

Refinement and validation of two digital Microwave Landing System (MLS) theoretical models

Two digital microwave landing system theoretical models are considered which are generic models for the Doppler and scanning-beam frequency reference versions of the MLS. These models represent errors resulting from both system noise and discrete multipath. The data used for the validation effort were obtained from the Texas Instrument conventional scanning beam and the Hazeltine Doppler feasibility hardware versions of the MLS. Topics discussed include tape read software, time history plots, computation of power spectral density, smoothed power spectra, best-fit models, different equations for digital simulation, and discrete multipath errors.

Duff, W. G.↗

The Corona of the Young Solar Analog EK Draconis

First coronal microwave and new soft X-ray observations of the very active, near-Zero-Age Main-Sequence (ZAMS) dGOe star EK Dra = HD 129333 show that this analog of the young Sun is more luminous in both emissions than most single M-dwarf flare stars. Variations in the 8.4 GHz flux include modulation with the optically determined rotation period of 2.7 days. This result points to a non-uniform filling of the corona with energetic electrons due to an incomplete coverage of the surface with active regions and a source volume that is not concentric with the star. The radio luminosity varying between log L(sub R) = 13.6 and 14.6 (L(sub R) in erg/s/Hz) shows evidence for unpolarized gyrosynchrotron flares, while strongly polarized flares were absent during the observations. This star is the first young, truly solar-like main sequence G star discovered in microwaves. Having just arrived on the main sequence, it conclusively proves that young, solar-like G stars can maintain very high levels of radio emission after their T Tau phase. The X-ray observations were obtained from the ROSAT All-Sky Survey (RASS). The average X-ray luminosity amounts to log L(sub x) = 29.9 (L(sub x) in erg/s). A Raymond-Smith type plasma model fit yields two plasma components at temperatures of 1.9 and 10 MK, with volume emission measures of 1.2 and 2.5 x 10 (exp 52)/cu cm, respectively. The X-ray light curve is significantly variable, with the photon count rate from the cooler plasma being strongly modulated by the rotation period; the emission from the hotter plasma is only weakly variable. Modeling of the source distribution in the stellar corona yields electron densities of the order of 4 x 10(exp 10)/cu cm or higher for the cool plasma component. It indicates that a considerable portion of EK Dra's high X-ray luminosity is due to high-density plasma rather than large emission volume. Parameters for an X-ray flare indicate an electron density of 1.75 x 10(exp 11)/cu cm and a source height of (1-2) x 10(exp 10) cm, compatible with a few times the scale height of the cooler plasma component.

Gudel, M.↗

A computationally efficient statistically downscaled 100 m resolution Greenland product from the regional climate model MAR

The Greenland Ice Sheet (GrIS) has been contributing directly to sea level rise, and this contribution is projected to accelerate over the next decades. A crucial tool for studying the evolution of surface mass loss (e.g., surface mass balance, SMB) consists of regional climate models (RCMs), which can provide current estimates and future projections of sea level rise associated with such losses. However, one of the main limitations of RCMs is the relatively coarse horizontal spatial resolution at which outputs are currently generated. Here, we report results concerning the statistical downscaling of the SMB modeled by the Modèle Atmosphérique Régional (MAR) RCM from the original spatial resolution of 6 km to 100 m building on the relationship between elevation and mass losses in Greenland. To this goal, we developed a geospatial framework that allows the parallelization of the downscaling process, a crucial aspect to increase the computational efficiency of the algorithm. Using the results obtained in the case of the SMB, surface and air temperature are assessed through the comparison of the modeled outputs with in situ and satellite measurement. The downscaled products show a considerable improvement in the case of the downscaled product with respect to the original coarse output, with the coefficient of determination (R 2 ) increasing from 0.868 for the original MAR output to 0.935 for the SMB downscaled product. Moreover, the value of the slope and intercept of the linear regression fitting modeled and measured SMB values shifts from 0.865 for the original MAR to 1.015 for the downscaled product in the case of the slope and from the value −235 mm w.e. yr -1 (original) to −57 mm w.e. yr -1 (downscaled) in the case of the intercept, considerably improving upon results previously published in the literature.

Greenland Ice Sheet↗

Detection of the SO2 atmosphere on Io with the Hubble Space Telescope

Observations of the trailing hemisphere of Io made with the Faint Object Spectrograph of the Hubble Space Telescope (HST) in March 1992 have resulted in the first detection of atmospheric SO2 absorption bands in the ultraviolet. These observations represent only the third positive means of detection of what is widely believed to be Io's primary atmospheric constituent. Below approximately 2130 A the geometric albedo of the satellite is dominated by SO2 gas absorption band signatures, which have been analyzed using models that include the effects of optical thickness, temperature, and spatial distribution. The disk-intergrated HST data cannot resolve the spatial distribution, but it is possible to define basic properties and set constraints on the atmosphere at the time of the observations. Hemispheric atmospheres with average column density N = 6 - 10 x 10(exp 15)/sq sm and T(gas) = 110 - 500 K fit the data, with preference for temperatures of approximately 200 - 250 K. Better fits are found as the atmosphere is spatially confined, with a limit of approximately 8% hemispheric areal coverage and N approximately equal to 3 x 10(exp 17)/sq cm with colder 110 - 250 K temepratures. A dense (N greater than or equal to 10(exp 16)/sq cm), localized component of SO2 gas, such as that possibly associated with active volcanoes, can generate the observed spectral constrast only when the atmosphere is cold (110 K) and an extended component such as Pele is included. The combination of a dense, localized atmosphere with a tenuous component (N less than 10(exp 16)/sq cm, either patchy or extended) also fits the data. In all cases the best fit models imply a disk-averaged column density larger than exospheric but approximately 10 - 30 times less than the previous upper limit from near-UV observations.

Ballester, G. E.↗