Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,243 records · Page 69

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Fermi/GBM Observations of SGRJ0501 + 4516 Bursts

We present our temporal and spectral analyses of 29 bursts from SGRJ0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T(sub 90) durations of the bursts can be fit with a log-normal distribution with a mean value of approx. 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T(sub 90)s estimated in count space (following a log-normal distribution with a mean value of approx. 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that E(sub peak) decreases with energy flux (and fluence) to a minimum of approx. 30 keV at F = 8.7 x 10(exp -6)erg/sq cm/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550 - 5418 and 1806 - 20. The isotropic luminosity, L(sub iso), corresponding to these flux values is roughly similar for all sources (0.4 - l.5 x 10(exp 40) erg/s.

Lin, Lin↗

Algorithm for Compressing Time-Series Data

An algorithm based on Chebyshev polynomials effects lossy compression of time-series data or other one-dimensional data streams (e.g., spectral data) that are arranged in blocks for sequential transmission. The algorithm was developed for use in transmitting data from spacecraft scientific instruments to Earth stations. In spite of its lossy nature, the algorithm preserves the information needed for scientific analysis. The algorithm is computationally simple, yet compresses data streams by factors much greater than two. The algorithm is not restricted to spacecraft or scientific uses: it is applicable to time-series data in general. The algorithm can also be applied to general multidimensional data that have been converted to time-series data, a typical example being image data acquired by raster scanning. However, unlike most prior image-data-compression algorithms, this algorithm neither depends on nor exploits the two-dimensional spatial correlations that are generally present in images. In order to understand the essence of this compression algorithm, it is necessary to understand that the net effect of this algorithm and the associated decompression algorithm is to approximate the original stream of data as a sequence of finite series of Chebyshev polynomials. For the purpose of this algorithm, a block of data or interval of time for which a Chebyshev polynomial series is fitted to the original data is denoted a fitting interval. Chebyshev approximation has two properties that make it particularly effective for compressing serial data streams with minimal loss of scientific information: The errors associated with a Chebyshev approximation are nearly uniformly distributed over the fitting interval (this is known in the art as the "equal error property"); and the maximum deviations of the fitted Chebyshev polynomial from the original data have the smallest possible values (this is known in the art as the "min-max property").

Hawkins, S. Edward, III↗

Analyses of Cometary Silicate Crystals: DDA Spectral Modeling of Forsterite

Comets are the Solar System's deep freezers of gases, ices, and particulates that were present in the outer protoplanetary disk. Where comet nuclei accreted was so cold that CO ice (approximately 50K) and other supervolatile ices like ethane (C2H2) were preserved. However, comets also accreted high temperature minerals: silicate crystals that either condensed (greater than or equal to 1400 K) or that were annealed from amorphous (glassy) silicates (greater than 850-1000 K). By their rarity in the interstellar medium, cometary crystalline silicates are thought to be grains that formed in the inner disk and were then radially transported out to the cold and ice-rich regimes near Neptune. The questions that comets can potentially address are: How fast, how far, and over what duration were crystals that formed in the inner disk transported out to the comet-forming region(s)? In comets, the mass fractions of silicates that are crystalline, f_cryst, translate to benchmarks for protoplanetary disk radial transport models. The infamous comet Hale-Bopp has crystalline fractions of over 55%. The values for cometary crystalline mass fractions, however, are derived assuming that the mineralogy assessed for the submicron to micron-sized portion of the size distribution represents the compositional makeup of all larger grains in the coma. Models for fitting cometary SEDs make this assumption because models can only fit the observed features with submicron to micron-sized discrete crystals. On the other hand, larger (0.1-100 micrometer radii) porous grains composed of amorphous silicates and amorphous carbon can be easily computed with mixed medium theory wherein vacuum mixed into a spherical particle mimics a porous aggregate. If crystalline silicates are mixed in, the models completely fail to match the observations. Moreover, models for a size distribution of discrete crystalline forsterite grains commonly employs the CDE computational method for ellipsoidal platelets (c:a:b=8.14x8.14xl in shape with geometrical factors of x:y:z=1:1:10, Fabian et al. 2001; Harker et al. 2007). Alternatively, models for forsterite employ statistical methods like the Distribution of Hollow Spheres (Min et al. 2008; Oliveira et al. 2011) or Gaussian Random Spheres (GRS) or RGF (Gielen et al. 200S). Pancakes, hollow spheres, or GRS shapes similar to wheat sheaf crystal habit (e.g., Volten et al. 2001; Veihelmann et al. 2006), however, do not have the sharp edges, flat faces, and vertices seen in images of cometary crystals in interplanetary dust particles (IDPs) or in Stardust samples. Cometary forsterite crystals often have equant or tabular crystal habit (J. Bradley). To simulate cometary crystals, we have computed absorption efficiencies of forsterite using the Discrete Dipole Approximation (DDA) DDSCAT code on NAS supercomputers. We compute thermal models that employ a size distribution of discrete irregularly shaped forsterite crystals (nonspherical shapes with faces and vertices) to explore how crystal shape affects the shape and wavelength positions of the forsterite spectral features and to explore whether cometary crystal shapes support either condensation or annealing scenarios (Lindsay et al. 2012a, b). We find forsterite crystal shapes that best-fit comet Hale-Bopp are tetrahedron, bricks or brick platelets, essentially equant or tabular (Lindsay et al. 2012a,b), commensurate with high temperature condensation experiments (Kobatake et al. 2008). We also have computed porous aggregates with crystal monomers and find that the crystal resonances are amplified. i.e., the crystalline fraction is lower in the aggregate than is derived by fitting a linear mix of spectral features from discrete subcomponents, and the crystal resonances 'appear' to be from larger crystals (Wooden et al. 2012). These results may indicate that the crystalline mass fraction in comets with comae dominated by aggregates may be lower than deduced by popular methods that only emoy ensembles of discrete crystals.

Wooden, Diane↗

Assessment of Mars Atmospheric Temperature Retrievals from the Thermal Emission Spectrometer Radiances

Motivated by the needs of Mars data assimilation. particularly quantification of measurement errors and generation of averaging kernels. we have evaluated atmospheric temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) radiances. Multiple sets of retrievals have been considered in this study; (1) retrievals available from the Planetary Data System (PDS), (2) retrievals based on variants of the retrieval algorithm used to generate the PDS retrievals, and (3) retrievals produced using the Mars 1-Dimensional Retrieval (M1R) algorithm based on the Optimal Spectral Sampling (OSS ) forward model. The retrieved temperature profiles are compared to the MGS Radio Science (RS) temperature profiles. For the samples tested, the M1R temperature profiles can be made to agree within 2 K with the RS temperature profiles, but only after tuning the prior and error statistics. Use of a global prior that does not take into account the seasonal dependence leads errors of up 6 K. In polar samples. errors relative to the RS temperature profiles are even larger. In these samples, the PDS temperature profiles also exhibit a poor fit with RS temperatures. This fit is worse than reported in previous studies, indicating that the lack of fit is due to a bias correction to TES radiances implemented after 2004. To explain the differences between the PDS and Ml R temperatures, the algorithms are compared directly, with the OSS forward model inserted into the PDS algorithm. Factors such as the filtering parameter, the use of linear versus nonlinear constrained inversion, and the choice of the forward model, are found to contribute heavily to the differences in the temperature profiles retrieved in the polar regions, resulting in uncertainties of up to 6 K. Even outside the poles, changes in the a priori statistics result in different profile shapes which all fit the radiances within the specified error. The importance of the a priori statistics prevents reliable global retrievals based a single a priori and strongly implies that a robust science analysis must instead rely on retrievals employing localized a priori information, for example from an ensemble based data assimilation system such as the Local Ensemble Transform Kalman Filter (LETKF).

Hoffman, Matthew J.↗

Communicating the Threat of a Tropical Cyclone to the Eastern Range

The 45th Weather Squadron (45 WS) has developed a tool to help visualize the Wind Speed Probability product from the National Hurricane Center (NHC) and to help communicate that information to space launch customers and decision makers at the 45th Space Wing (45 SW) and Kennedy Space Center (KSC) located in east central Florida. This paper reviews previous work and presents the new visualization tool, including initial feedback as well as the pros and cons. The NHC began issuing their Wind Speed Probability product for tropical cyclones publicly in 2006. The 45 WS uses this product to provide a threat assessment to 45 SW and KSC leadership for risk evaluations with an approaching tropical cyclone. Although the wind speed probabilities convey the uncertainty of a tropical cyclone well, communicating this information to customers is a challenge. The 45 WS continually strives to provide the wind speed probability information to customers in a context which clearly communicates the threat of a tropical cyclone. First, an intern from the Florida Institute of Technology (FIT) Atmospheric Sciences department, sponsored by Scitor Corporation, independently evaluated the NHC wind speed probability product. This work was later extended into a M.S. thesis at FIT, partially funded by Scitor Corporation and KSC. A second thesis at FIT further extended the evaluation partially funded by KSC. Using this analysis, the 45 WS categorized the probabilities into five probability interpretation categories: Very Low, Low, Moderate, High, and Very High. These probability interpretation categories convert the forecast probability and forecast interval into easily understood categories that are consistent across all ranges of probabilities and forecast intervals. As a follow-on project, KSC funded a summer intern to evaluate the human factors of the probability interpretation categories, which ultimately refined some of the thresholds. The 45 WS created a visualization tool to express the timing and risk for multiple locations in a single graphic. Preliminary results on an on-going project by FIT will be included in this paper. This project is developing a new method of assigning the probability interpretation categories and updating the evaluation of the performance of the NHC wind speed probability analysis.

Winters, Katherine A.↗

Comparison of the Variations of Sunspot Number, Number of Sunspot Groups, and Sunspot Area, 1875 -2013

Examined are the yearly variations and ratios of sunspot number, the number of sunspot groups, and the total corrected sunspot area for the interval 1875-2013. While yearly sunspot number independently correlates strongly (r = 0.98) with the yearly number of sunspot groups (y = -2 + 11.99x) and the total corrected sunspot area (y = 5 + 0.059x), the strongest correlation (Ry12 = 0.99) is the one based on the bivariate fit of sunspot number against the combined variations of the number of sunspot groups and sunspot area (y = 1 + 5.88x1 + 0.031x2, where y refers to sunspot number, x1 refers to the number of sunspot groups, and x2 refers to the sunspot area). While all cycle minima based on the bivariate fit are concurrent with the observed minimum in sunspot number, cycle maxima are sometimes found to differ. For sunspot cycles 12, 19, 20, and 23, cycle maximum is inferred to have occurred in 1884, 1958, 1970, and 2002, respectively, rather than in 1883, 1957, 1968, and 2000, based on the observed sunspot number. Also, cycle 19's maximum amplitude based on observed sunspot number seems too high in comparison to that found using the bivariate fit. During the 139-year interval 1875-2013, the difference between the observed and predicted sunspot number based on the bivariate fit is <1 standard error of estimate (se) (<6.4) for 111 years, between 1 and <2 se (6.4 to <12.8) for 28 years, and ≥2 se (≥12.8) for only 4 years, these years being 1957 (16.6), 1978 (-15.8), 1980 (23), and 1982 (-16.3). For sunspot cycle 24, the difference between observed and predicted values has been only -0.7 and 3.2 (≤0.5 se).

Wilson, Robert M.↗

A New Global Regression Analysis Method for the Prediction of Wind Tunnel Model Weight Corrections

A new global regression analysis method is discussed that predicts wind tunnel model weight corrections for strain-gage balance loads during a wind tunnel test. The method determines corrections by combining "wind-on" model attitude measurements with least squares estimates of the model weight and center of gravity coordinates that are obtained from "wind-off" data points. The method treats the least squares fit of the model weight separate from the fit of the center of gravity coordinates. Therefore, it performs two fits of "wind- off" data points and uses the least squares estimator of the model weight as an input for the fit of the center of gravity coordinates. Explicit equations for the least squares estimators of the weight and center of gravity coordinates are derived that simplify the implementation of the method in the data system software of a wind tunnel. In addition, recommendations for sets of "wind-off" data points are made that take typical model support system constraints into account. Explicit equations of the confidence intervals on the model weight and center of gravity coordinates and two different error analyses of the model weight prediction are also discussed in the appendices of the paper.

strain-gage balance↗

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗

Ultrasonic Measurement and Monitoring of Loads in Bolts used in Structural Joints

The paper is an overview of work by the author in measuring and monitoring loads in bolts using an ultrasonic extensometer. A number of cases of bolted joints are covered. These include, a clamped joint with clearance fit between the bolt and hole, a clamped joint with the bolt in an interference fit with the hole, a flanged joint which allows the flange and bolt to bend; and a shear joint in a clevis and tang configuration. These applications were initially developed for measuring and monitoring preload in the NASA Space Shuttle Orbiter critical joints but are also applicable for monitoring loads in other critical bolted joints of structures such as transportation bridges and other aerospace structures. The paper explains how to set‐up a model to estimate the load factor and accuracy for the ultrasonic preload application in a clamped joint with clearance fit. The ultrasonic preload application for clamped joint with bolt in an interference fit can also be used to measure diametrical interference between the bolt shank and hole; and interference pressure on the bolt shank. Model and experimental data are given to demonstrate use of ultrasonic measurements in a shear joint. A bolt in a flanged joint experiences both tensile and bending loads. This application involves measurement of bending and tensile preload in a bolt. The ultrasonic beam bends due to the bending load on the bolt. A numerical technique to compute the trace of ultrasonic ray is presented.

Koshti, Ajay↗

Quantification of ln-Flight Physical Changes: Anthropometry and Neutral Body Posture

Currently, NASA does not have sufficient in-flight anthropometric data to assess the impact of changes in body shape and size. For developing future planetary and reduced-gravity suits, NASA needs to quantify the impacts of microgravity on anthropometry and body posture to ensure optimal crew performance, fit, and comfort. To obtain data on these changes, circumference, length, height, breadth, and depth for body segments (chest, waist, bicep, thigh, calf) from astronauts for preflight, in-flight, and post-flight conditions needs to be collected. Once these data have been collected, pre-flight, in-flight, and post-flight anthropometric values will be compared, yielding microgravity factors. The neutral body posture (NBP) will also be measured, to determine body posture (joint angle) changes between subjects throughout the duration of a mission. Data collection, starting with Increments 37/38, is still in progress but has been completed for 6 out of 9 subjects. NASA suit engineers and NASA's Extravehicular Activity (EVA) Project Office have identified that suit fit in microgravity could become an issue. It has been noted that crewmembers often need to adjust their suit sizing once they are in orbit. This adjustment could be due to microgravity effects on anthropometry and postural changes, and is necessary to ensure optimal crew performance, fit, and comfort in space. To date, the only data collected to determine the effects of microgravity on physical human changes were collected during Skylab 4, the Apollo-Soyuz Test Project (ASTP), Space Shuttle mission STS-57, and a recent HRP study on seated height changes due to spinal elongation (Spinal Elongation, Master Task List [MTL] #221). The Skylab 4, ASTP, and the STS-57 studies found that, according to photographs, a distinct NBP exists. The still photographs showed a distinguishable posture with the arms raised and the shoulders abducted; in addition, the knees are flexed, with noticeable hip flexion, and the foot is plantar flexed [1,2]. This combination is considered to be the standard set of body joint angles for an NBP in microgravity. A recent simulated microgravity NBP study [3] showed individual variability and inconsistencies in defining NBP. This variation may be influenced by spinal growth, and other potential anthropometry factors such as spinal curvature, age, and gender. Data on the variation of this posture data is required for all kinds of space device designs (such as suits, habitat, and mobility aids). The method proposed in this study considers the dynamic nature of body movement and will use a measurement technique to continually monitor posture and develop a probability likelihood of the neutral posture and how the NBP postures are affected by anthropometry. Additionally, Skylab studies found that crewmembers experienced a stature growth of up to 3 percent. The data included 3 crewmembers and showed that a biphasic stature growth occurs once the crewmember enters into weightlessness. However, the HRP Spinal Elongation study showed that crewmembers could experience about a 6 percent growth in seated height and a 3 percent stature growth, when exposed to microgravity. The results of that study prove that not all anthropometric measurements have the same microgravity percent growth factor. For EVA and suit engineers to properly update the sizing protocol for microgravity, they need additional anthropometric data from space missions. Hence, this study is aimed at gathering additional in-flight anthropometric measurements, such as length, depth, breadth, and circumference, to determine the changes to body shape and size caused by microgravity effects. It is anticipated that by recording the potential changes to body shape and size, NASA will develop a better suit sizing protocol for the International Space Station and other space missions. In essence, this study will help NASA quantify the impacts of microgravity on anthropometry to ensure optimal crew performance, fit, and comfort. This study will use standard anthropometry data collection techniques, 3D laser scanning, digital still photography, and video data, and perform photogrammetric analyses to determine the changes that occur to the body shape and size, and to NBP, when the human body is exposed to a microgravity environment.

Young, K. S.↗

Mission X: Train Like an Astronaut Challenge

The Mission X: Train Like an Astronaut Challenge was developed in 2011 to encourage proper exercise and nutrition at an early age by teaching young people to live and eat like space explorers. The strong correlation between an unhealthy childhood diet and adolescent fitness, and the onset of chronic diseases as an adult is the catalyst for Mission X. Mission X is dedicated to assisting people on a global scale to live healthier lifestyles and learn about human space exploration. The Mission X: Train Like an Astronaut 2015 (MX15) International Challenge hosted almost 40,000 children on 800 teams, 28 countries affiliated with 12 space agencies. The MX15 website included 17 languages. MX15, the fifth annual international fitness challenges sponsored by the NASA Human Research Program worked with the European Space Agency and other space agencies from around the world. In comparison to MX14, MX15 expanded to include four additional new countries, increased the number of students by approximately 68% and the number of teams by 29%. Chile' and South Korea participated in the new fall Astro Charlie Walk Around the Earth Challenge. Pre-challenge training materials were made more readily available from the website. South Korea completed a prospective assessment of the usability of the MX content for improving health and fitness in 212 preschool children and their families. Mission X is fortunate to have the support of the NASA, ESA and JAXA astronaut corps. In MX15, they participated in the opening and closing events as well as while on-board the International Space Station. Italian Astronaut Samantha Cristoretti participated as the MX15 Astronaut Ambassador for health and fitness providing the opening video and other videos from ISS. United Kingdom Astronaut Tim Peake and US Astronaut Kate Rubins have agreed to be the MX Ambassadors for 2016 and 2017 respectively. The MX15 International Working Group Face-to-Face meeting and Closing Event were held at the Agenzia Spaziale Italiana (ASI) in Rome, Italy. A record number of twenty-eight countries participated. Austria and Norway have offered to host the 2016 and 2017 working group meetings. MX16 planning began with the working group meetings and areas of improvement will include another second early challenge to accommodate countries in the Southern Hemisphere, recommended changes to the MX website, development of a more defined approach to metrics, a change to the format of future MX International Working Group meetings, and proposed new activities to be developed by the MX International Educator Working Group. We look forward to welcoming many new participants in 2016!

Lloyd, Charles W.↗

Vertical Structure and Optical Properties of Titans Aerosols from Radiance Measurements Made Inside and Outside the Atmosphere

Prompted by the detection of stratospheric cloud layers by Cassini's Composite Infrared Spectrometer (CIRS; see Anderson, C.M., Samuelson, R.E. [2011]. Icarus 212, 762-778), we have re-examined the observations made by the Descent Imager/Spectral Radiometer (DISR) in the atmosphere of Titan together with two constraints from measurements made outside the atmosphere. No evidence of thin layers (<1 km) in the DISR image data sets is seen beyond the three previously reported layers at 21 km, 11 km, and 7 km by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M.G. [2009]. Icarus 199, 442-448). On the other hand, there is evidence of a thicker layer centered at about 55 km. A rise in radiance gradients in the Downward-Looking Visible Spectrometer (DLVS) data below 55 km indicates an increase in the volume extinction coefficient near this altitude. To fit the geometric albedo measured from outside the atmosphere the decrease in the single scattering albedo of Titan's aerosols at high altitudes, noted in earlier studies of DISR data, must continue to much higher altitudes. The altitude of Titan's limb as a function of wavelength requires that the scale height of the aerosols decrease with altitude from the 65 km value seen in the DISR observations below 140 km to the 45 km value at higher altitudes. We compared the variation of radiance with nadir angle observed in the DISR images to improve our aerosol model. Our new aerosol model fits the altitude and wavelength variations of the observations at small and intermediate nadir angles but not for large nadir angles, indicating an effect that is not reproduced by our radiative transfer model. The volume extinction profiles are modeled by continuous functions except near the enhancement level near 55 km altitude. The wavelength dependence of the extinction optical depth is similar to earlier results at wavelengths from 500 to 700 nm, but is smaller at shorter wavelengths and larger toward longer wavelengths. A Hapke-like model is used for the ground reflectivity, and the variation of the Hapke single scattering albedo with wavelength is given. Fits to the visible spectrometers looking upward and downward are achieved except in the methane bands longward of 720 nm. This is possibly due to uncertainties in extrapolation of laboratory measurements from 1 km-am paths to much longer paths at lower pressures. It could also be due to changes in the single scattering phase functions at low altitudes, which strongly affect the path length through methane that the photons travel. We demonstrate the effects on the model fits by varying each model parameter individually in order to illustrate the sensitivity of our determination of each model parameter.

Titan↗

NuSTAR AND Swift Observations of the Very High State in GX 339-4: Weighing the Black Hole With X-Rays

We present results from spectral fitting of the very high state of GX339-4 with Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift. We use relativistic reflection modeling to measure the spin of the black hole and inclination of the inner disk and find a spin of a = 0.95+0.08/-0.02 and inclination of 30deg +/- 1deg (statistical errors). These values agree well with previous results from reflection modeling. With the exceptional sensitivity of NuSTAR at the high-energy side of the disk spectrum, we are able to constrain multiple physical parameters simultaneously using continuum fitting. By using the constraints from reflection as input for the continuum fitting method, we invert the conventional fitting procedure to estimate the mass and distance of GX 339-4 using just the X-ray spectrum, finding a mass of 9.0+1.6/-1.2 Stellar Mass and distance of 8.4 +/- 0.9 kpc (statistical errors).

accretion↗

The NuSTAR View of the Non-Thermal Emission from PSR J0437-4715

We present a hard X-ray Nuclear Spectroscopic Telescope Array (NuSTAR) observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period approximately 5.76 ms are observed with a significance of 3.7sigma, at energies up to 20 keV above which the NuSTAR background dominates. We measure a photon index gamma = 1.50 +/- 0.25(90 per cent confidence) for the power-law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsars phase-folded light curve with the pulsars well-defined mass and distance from radio timing observations.

stars; neutro↗

Radiative Heating Indicators and Uncertainty Analysis for Mars 2020

This paper will present analysis to evaluate a functional form and associated parameters to determine the radiative heating at various back-shell locations of the Mars 2020 aeroshell. The radiative heating indicators are used for preliminary vehicle design, and to determine the worst-case trajectories for peak radiative heat flux and load. Historically, the functional form for radiative heating is based on free-stream parameters of density, velocity and a length scale, typically the nose-radius (or equivalent). However, a fit of this form has been shown to not provide significant enough accuracy when compared to simulation results. Therefore, a fit based on post-shock equilibrium calculations of CO and CO2 number density, temperature and pressure has been devised. The results from the TPS 15 01 trajectory were used to develop the fit. The fit was then applied to the MSL best estimated trajectory (BET). Furthermore, the paper will quantify the uncertainty in simulations of the radiative heating for Mars 2020.

Mars 2020↗

GOES-17 Magnetometer On-Orbit Calibration Design and Results

The Geostationary Operational Environmental Satellite-R (GOES-R) is a series of four NOAA weather satellites. Each spacecraft flies a Magnetometer instrument which includes two boom mounted fluxgate sensors located 6.3 and 8.5 meters from the spacecraft. Sensor readings are corrupted by both internal (bias and scale factor) and external (spacecraft emitted fields and misalignments) error sources. In order to estimate the bias, an on-orbit calibration is performed, slewing the spacecraft about multiple axes and then performing a fit of the data to determine zero offset and misalignment terms. The first satellite in the series, GOES-16 was launched in November 2016. Its magnetometer calibration maneuver was found to be inadequate providing insufficient rotation about orthogonal axes. The uncertainty in the estimated biases was found to be in the 3-5 nanotesla range, exceeding the magnetometer accuracy requirement of 1.0 nanotesla. Based on lessons learned, the maneuver was redesigned to provide multiple 360 degree rotations about two orthogonal axes, providing optimal visibility into bias and misalignment. Varying ambient fields over the duration of the calibration maneuvers, corrupts the ability to accurately determine parameters of interest. To address this issue, rather than assuming a constant field for the duration of the maneuver, a non-linear least squares filter was employed that uses a series of spline functions to better fit the ambient field. The number of splines can be varied to minimize the least squares filter residuals. Simulations were performed using 360 days of GOES-16 data to bound the uncertainty in the biases for given residuals of the fit to the ambient field. Following the launch of the second satellite in the series, GOES-17 in March 2018, the new calibration maneuver and enhanced bias estimation algorithm were exercised twice. Based on the residuals of the fit, the uncertainty of the estimated zero offset was reduced by an order of magnitude from 3-5 nanotesla to approximately 0.4 nanotesla per axis. This paper provides details of the maneuver design and the innovative algorithm used to estimate the GOES-17 biases and misalignments.

GOES-17↗

NuStar Hard X-Ray View of Low-Luminosity Active Galactic Nuclei: High-Energy Cutoff and Truncated Thin Disk

We report the analysis of simultaneous XMM-Newton+Nuclear Spectroscopic Telescope Array (NuSTAR) observations of two low-luminosity active galactic nuclei (LLAGNs), NGC3998 and NGC4579. We do not detect any significant variability in either source over the ∼3 day length of the NuSTAR observations. The broadband 0.5–60 keV spectrum of NGC3998 is best fit with a cutoff power law, while the one for NGC4579 is best fit with a combination of a hot thermal plasma model, a power law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with 3σ reflection fraction upper limits of R < 0.3 and R < 0.18 for NGC3998 and NGC4579, respectively; (2) the 6–7 keV line complex in NGC4579 could be fit with either a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line or with three relatively narrow lines, which include contribution from Fe XXVI; (3) the NGC4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy = - E 107+ cut 18 27 keV in NGC3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) the NGC3998 spectrum is consistent with a Comptonization model with either a sphere (τ ≈ 3 ± 1) or slab (τ ≈ 1.2 ± 0.6) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGNs, other LLAGNs, and hot accretion flow models.

Younes, George↗