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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Prompt Emission of GRB 121217A from Gamma-Rays to the Near-Infrared

The mechanism that causes the prompt-emission episode of gamma-ray bursts (GRBs) is still widely debated despite there being thousands of prompt detections. The favoured internal shock model relates this emission to synchrotron radiation. However, it does not always explain the spectral indices of the shape of the spectrum, which is often fit with empirical functions, such as the Band function. Multi-wavelength observations are therefore required to help investigate the possible underlying mechanisms that causes the prompt emission. We present GRB 121217A, for which we were able to observe its near-infrared (NIR) emission during a secondary prompt-emission episode with the Gamma-Ray burst Optical Near-infrared Detector (GROND) in combination with the Swift and Fermi satellites, which cover an energy range of 5 orders of magnitude (10(exp −3) keV to 100 keV). We determine a photometric redshift of z = 3.1 +/- 0.1 with a line-of-sight with little or no extinction (AV approx. 0 mag) utilising the optical/NIR SED. From the afterglow, we determine a bulk Lorentz factor of Gamma(sub 0) approx. 250 and an emission radius of R < 1018 cm. The prompt-emission broadband spectral energy distribution is well fit with a broken power law with beta1 = −0.3 +/- 0.1 and beta2 = 0.6 +/- 0.1 that has a break at E = 6.6 +/- 0.9 keV, which can be interpreted as the maximum injection frequency. Self-absorption by the electron population below energies of Ea < 6 keV suggest a magnetic field strength of B approx. 10(exp 5) G. However, all the best fit models underpredict the flux observed in the NIR wavelengths, which also only rebrightens by a factor of approx. 2 during the second prompt emission episode, in stark contrast to the X-ray emission, which rebrightens by a factor of approx. 100. This suggests an afterglow component is dominating the emission. We present GRB 121217A, one of the few GRBs that has multi-wavelength observations of the prompt-emission period and shows that it can be understood with a synchrotron radiation model. However, due to the complexity of the GRB's emission, other mechanisms that result in Band-like spectra cannot be ruled out.

prompt↗
NASA NTRS2014

X-Ray Spectral Components Observed in the Afterglow of GRB 130925A

We have identified spectral features in the late-time X-ray afterglow of the unusually long, slow-decaying GRB 130925A using NuSTAR, Swift/X-Ray Telescope, and Chandra. A spectral component in addition to an absorbed power law is required at greater than 4 less than 1 significance, and its spectral shape varies between two observation epochs at 2 x 10 (sup 5) and 10 (sup 6) seconds after the burst. Several models can fit this additional component, each with very different physical implications. A broad, resolved Gaussian absorption feature of several kiloelectronvolts width improves the fit, but it is poorly constrained in the second epoch. An additive blackbody or second power-law component provide better fits. Both are challenging to interpret: the blackbody radius is near the scale of a compact remnant (10 (sup 8) centimeters), while the second power-law component requires an unobserved high-energy cutoff in order to be consistent with the non-detection by Fermi/Large Area Telescope.

SPECTRAL↗
NASA NTRS2015

Sub-GLE Solar Particle Events and the Implications for Lightly-Shielded Systems Flown During an Era of Low Solar Activity

Many of the large space missions must be very rigorous in their designs to reduce risk from radiation damage as much as possible. Some ways of reducing this risk have been to build in multiple redundancies, purchase/develop radiation hardened electronics parts, and plan for worst case radiation environment scenarios. These methods work well for these ambitious missions that can afford the costs associated with these meticulous efforts. However, there have been more small spacecraft and CubeSats with smaller duration missions entering the space arena, which can take some additional risks, but cannot afford to implement all of these risk-reducing methods. Therefore, one way to quantify the radiation exposure risk for these smaller spacecraft would be to investigate the radiation environment pertinent to the mission to better understand these radiation exposures, rather than always designing to the infrequent, worst-case environment. In this study, we have investigated 34 historical solar particle events (1974-2010) that occurred during a time period when the sun spot number (SSN) was less than 30. These events contain Ground Level Events (GLE), sub-GLEs, and sub-sub-GLEs(sup 1-3). GLEs are extremely energetic solar particle events (SPEs) having proton energies often extending into the several GeV range and producing secondary particles in the atmosphere, mostly neutrons, observed with ground station neutron monitors. Sub-GLE events are less energetic, extending into the several hundred MeV range, but without producing detectable levels of secondary atmospheric particles. Sub-sub GLEs are even less energetic with an observable increase in protons at energies greater than 30 MeV, but no observable proton flux above 300 MeV. The spectra for these events were fitted using a double power law fit in particle rigidity, called the Band fit method. The differential spectra were then input into the NASA Langley Research Center HZETRN 2005, which is a high-energy particle transport/dose code, to determine the dose in various thicknesses of aluminum, representing the spacecraft. This paper will detail the absorbed dose results of each of these environments, as well as analyze the data to better understand the doses over small thicknesses that are more relevant to small spacecraft and satellites, such as CubeSats. In addition, we will discuss the implications of these data and provide some recommendations that may be useful to spacecraft designers of these smaller spacecraft/CubeSat-type missions.

Atwell, William↗
NASA NTRS2013

Estimating Ground-Level PM(sub 2.5) Concentrations in the Southeastern United States Using MAIAC AOD Retrievals and a Two-Stage Model

Previous studies showed that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with various health outcomes. Ground in situ measurements of PM(sub 2.5) concentrations are considered to be the gold standard, but are time-consuming and costly. Satellite-retrieved aerosol optical depth (AOD) products have the potential to supplement the ground monitoring networks to provide spatiotemporally-resolved PM(sub 2.5) exposure estimates. However, the coarse resolutions (e.g., 10 km) of the satellite AOD products used in previous studies make it very difficult to estimate urban-scale PM(sub 2.5) characteristics that are crucial to population-based PM(sub 2.5) health effects research. In this paper, a new aerosol product with 1 km spatial resolution derived by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was examined using a two-stage spatial statistical model with meteorological fields (e.g., wind speed) and land use parameters (e.g., forest cover, road length, elevation, and point emissions) as ancillary variables to estimate daily mean PM(sub 2.5) concentrations. The study area is the southeastern U.S., and data for 2003 were collected from various sources. A cross validation approach was implemented for model validation. We obtained R(sup 2) of 0.83, mean prediction error (MPE) of 1.89 micrograms/cu m, and square root of the mean squared prediction errors (RMSPE) of 2.73 micrograms/cu m in model fitting, and R(sup 2) of 0.67, MPE of 2.54 micrograms/cu m, and RMSPE of 3.88 micrograms/cu m in cross validation. Both model fitting and cross validation indicate a good fit between the dependent variable and predictor variables. The results showed that 1 km spatial resolution MAIAC AOD can be used to estimate PM(sub 2.5) concentrations.

aerosol optical depth↗
NASA NTRS2013

Bayesian Revision of Residual Detection Power

This paper addresses some issues with quality assessment and quality assurance in response surface modeling experiments executed in wind tunnels. The role of data volume on quality assurance for response surface models is reviewed. Specific wind tunnel response surface modeling experiments are considered for which apparent discrepancies exist between fit quality expectations based on implemented quality assurance tactics, and the actual fit quality achieved in those experiments. These discrepancies are resolved by using Bayesian inference to account for certain imperfections in the assessment methodology. Estimates of the fraction of out-of-tolerance model predictions based on traditional frequentist methods are revised to account for uncertainty in the residual assessment process. The number of sites in the design space for which residuals are out of tolerance is seen to exceed the number of sites where the model actually fails to fit the data. A method is presented to estimate how much of the design space in inadequately modeled by low-order polynomial approximations to the true but unknown underlying response function.

DeLoach, Richard↗
NASA NTRS2015

Relation Between Residual and Hoop Stresses and Rolling Bearing Fatigue Life

Rolling-element bearings operated at high speed or high vibration may require a tight interference fit between the bore of the bearing and shaft to prevent rotation of the bearing bore around the shaft and fretting damage at the interfaces. Previous work showed that the hoop stresses resulting from tight interference fits can reduce bearing lives by as much as 65 percent. Where tight interference fits are required, case-carburized steel such as AISI 9310 or M50 NiL is often used because the compressive residual stresses inhibit subsurface crack formation and the ductile core inhibits inner-ring fracture. The presence of compressive residual stress and its combination with hoop stress also modifies the Hertz stress-life relation. This paper analyzes the beneficial effect of residual stresses on rolling-element bearing fatigue life in the presence of high hoop stresses for three bearing steels. These additional stresses were superimposed on Hertzian principal stresses to calculate the inner-race maximum shearing stress and the resulting fatigue life of the bearing. The load-life exponent p and Hertz stress-life exponent n increase in the presence of compressive residual stress, which yields increased life, particularly at lower stress levels. The Zaretsky life equation is described and is shown to predict longer bearing lives and greater load- and stress-life exponents, which better predicts observed life of bearings made from vacuum-processed steel.

Fatigue↗
NASA NTRS2016

Evaluation of the Hinge Moment and Normal Force Aerodynamic Loads from a Seamless Adaptive Compliant Trailing Edge Flap in Flight

A seamless adaptive compliant trailing edge (ACTE) flap was demonstrated in flight on a Gulfstream III aircraft at the NASA Armstrong Flight Research Center. The trailing edge flap was deflected between minus 2 deg up and plus 30 deg down in flight. The safety-of-flight parameters for the ACTE flap experiment require that flap-to-wing interface loads be sensed and monitored in real time to ensure that the structural load limits of the wing are not exceeded. The attachment fittings connecting the flap to the aircraft wing rear spar were instrumented with strain gages and calibrated using known loads for measuring hinge moment and normal force loads in flight. The safety-of-flight parameters for the ACTE flap experiment require that flap-to-wing interface loads be sensed and monitored in real time to ensure that the structural load limits of the wing are not exceeded. The attachment fittings connecting the flap to the aircraft wing rear spar were instrumented with strain gages and calibrated using known loads for measuring hinge moment and normal force loads in flight. The interface hardware instrumentation layout and load calibration are discussed. Twenty-one applied calibration test load cases were developed for each individual fitting. The 2-sigma residual errors for the hinge moment was calculated to be 2.4 percent, and for normal force was calculated to be 7.3 percent. The hinge moment and normal force generated by the ACTE flap with a hinge point located at 26-percent wing chord were measured during steady state and symmetric pitch maneuvers. The loads predicted from analysis were compared to the loads observed in flight. The hinge moment loads showed good agreement with the flight loads while the normal force loads calculated from analysis were over-predicted by approximately 20 percent. Normal force and hinge moment loads calculated from the pressure sensors located on the ACTE showed good agreement with the loads calculated from the installed strain gages.

Hinge Moment Loads↗
NASA NTRS2016

Forward Skirt Structural Testing on the Space Launch System (SLS) Program

Structural testing was performed to evaluate heritage forward skirts from the Space Shuttle program for use on the Space Launch System (SLS) program. One forward skirt is located in each solid rocket booster. Heritage forward skirts are aluminum 2219 welded structures. Loads are applied at the forward skirt thrust post and ball assembly. Testing was needed because SLS ascent loads are roughly 40% higher than Space Shuttle loads. Testing objectives were to determine margins of safety, demonstrate reliability, and validate analytical models. Two forward skirts were structurally tested using the test configuration. The test stand applied loads to the thrust post. Four hydraulic actuators were used to apply axial load and two hydraulic actuators were used to apply radial and tangential loads. The first test was referred to as FSTA-1 (Forward Skirt Structural Test Article) and was performed in April/May 2014. The purpose of FSTA-1 was to verify the ultimate capability of the forward skirt subjected to ascent ultimate loads. Testing consisted of two liftoff load cases taken to 100% limit load followed by an ascent load case taken to 110% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The second test was referred to as FSTA-2 and performed in July/August of 2014. The purpose of FSTA-2 was to verify the ultimate capability of the forward skirt subjected to liftoff ultimate loads. Testing consisted of six liftoff load cases taken to 100% limit load followed by the six liftoff cases taken to 140% limit load. Two ascent load cases were then tested to 100% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The forward skirts on FSTA-1 and FSTA-2 successfully carried all applied liftoff and ascent load cases. Both FSTA-1 and FSTA-2 were tested to failure by increasing the ascent loads. Failure occurred in the forward skirt thrust post radius. The forward skirts on FSTA-1 and FSTA-2 had nearly identical failure modes. FSTA-1 failed at 1.72 times limit load and FSTA-2 failed at 1.62 times limit load. This difference is primarily attributed to variation in material properties in the thrust post region. Test data were obtained from strain gages, deflection gages, ARAMIS digital strain measurement, acoustic emissions, and high-speed video. Strain gage data and ARAMIS strain were compared to finite element (FE) analysis predictions. Both the forward skirt and tooling were modeled. This allows the analysis to simulate the loading as close as possible to actual test configuration. FSTA-1 and FSTA-2 were instrumented with over 200 strain gages to ensure all possible failure modes could be captured. However, it turned out that three gages provided critical strain data. One was located in the post bore and two on the post radius. More gages were not specified due to space limitations and the desire to not interfere with the use of the ARAMIS system on the post radius. Measured strains were compared to analysis results for the load cycle to failure. Note that FSTA-1 gages were lost before failure was reached. FSTA-2 gages made it to the failure load but one of the radius gages was lost before testing began. This gage was not replaced because of the time and cost associated with disassembly of the test structure. Correlation to analysis was excellent for FSTA-1. FSTA-2 was not quite as good because there was more residual strain from previous load cycles. FSTA-2 was loaded and unloaded with 12 liftoff cases and two ascent cases before taking the skirt to failure. FSTA-1 only had two liftoff cases and one ascent case before taking the skirt to failure. The ARAMIS system was used to determine strain at the post radius by processing digital images of a speckled paint pattern. Digital cameras recorded images of the speckled paint pattern. ARAMIS strain results for FSTA-2 just prior to failure. Note a high strain location develops near the left side. This high strain compares well to analysis prediction for both FSTA-1 and FSTA-2. The strain at this location was also plotted versus limit load. Both FSTA-1 and FSTA-2 had excellent correlation between ARAMIS and analysis strains. Acoustic emission (AE) sensors were used to monitor for damage formation that may occur during testing (e.g., crack formation and growth or propagation). AE was very important because after disassembly of FSTA-1, a crack was observed in the ball fitting radius. The ball fitting did not crack on FSTA-2. AE data was used to reconstruct when the crack occurred. The AE energy versus time plot for FSTA. The energy increased considerably at 850 seconds (152% limit load), indicating a crack could have formed at this point. The only visual evidence found that could have corresponded to this was the crack that initiated in the ball fitting. The cracks in the forward skirt aluminum structures would likely have been lower energy due to a lower modulus and all that were found after failure correlated to occurring after the initial crack in the post radius. This was verified by high-speed cameras used to record the failure.

Lohrer, J. D.↗
NASA NTRS2016

CryoCart Restoration and Vacuum Pipe Construction

Propulsion systems that utilize hypergolic propellants have been used to power space vehicles since the beginning of the space program. Liquid methane and oxygen propulsion systems have emerged as an alternative and have proven to be more environmentally friendly. The incorporation of liquid methane/liquid oxygen (LOX) into the propulsion system has demonstrated an increase in engine performance, as well as a reduction in the volume, size and complexity of the system. Consequently, reducing the total mass of the vehicle which is a crucial aspect that is considered when planning space missions to both the Moon and Mars [1]. Project Morpheus has made significant advancements in liquid oxygen/liquid methane propulsion system technologies by incorporating a LOX/methane propulsion system to a vertical test bed. The vehicle consisted of a 5,000 lb main engine and four 20 lb remote control system (RCS) engines that utilize liquid methane/LOX as its propellant [1]. The vehicle completed successful flight testing at Kennedy Space Center in 2014 which marked the completion of the Morpheus project. Subsequent projects utilizing Morpheus' vertical test bed have been developed to make further advancements. One of the subsequent projects consisted of the addition of a smaller 2,000 lb main engine and a cold helium heat exchanger which would make it possible for a pressurant tank systems to be send to Mars or the Moon by significantly decreasing the overall mass and volume of the pressurant tank. The hot fire tests of the integrated system with the smaller main engine and cold helium heat exchanger were successful at sea level, but further studies are being conducted to better understand how the vertical test bed will behave under thermal-vacuum conditions. For this reason, the integrated vehicle will be taken to Plum Brook to be tested in a chamber capable of simulating these conditions. To ensure that the vehicle will function properly under vacuum conditions, testing will be first completed at the component level. During this process, the igniter of the main engine and the RCS thrusters will be tested under a vacuum. To complete the testing of the components, the test setup first needed to be finalized. The CryoCart is being used to feed the propellants to the test article. The CryoCart is a movable test set-up that was developed in 2009 to provide a mobile platform for testing oxygen/methane systems with hot-fire capability up to 100 lbf. The CryoCart consists of three different systems: Oxygen, Methane, and liquid Nitrogen. The Oxygen and Methane systems are placed into two different carts while the liquid nitrogen system is mainly located in the methane cart. Over the years, the CryoCart has been utilized for different projects and has undergone deterioration. For this reason, a new phase has been developed to rebuild it to working conditions once again. During my internship, I was aiding in the construction and restoration of the CryoCart. In the initial stages of the process, I updated the fluid and electrical schematics for the oxygen, methane, and test article systems. The original CryoCart consisted of an electrical panel that utilized electromechanical relays and a terminal to drive the igniter power and signal, as well as the main fuel and oxygen valves. This electrical panel connected to the CryoCart through various wire harnesses that could be found exiting from the CryoCart. First, it was determined how these harnesses connected to the electromechanical relays so that they worked correctly. Once the electrical system was understood, an alternative for the electromechanical relays and the Molex connectors used throughout the system was sought since these components can often prove to be unreliable. Solid State relays and MIL connectors were purchased to serve as replacements. Upon arrival of the parts, crimping and wiring was completed to install the new solid state relays and MIL connectors. During the replacement of the relays and connectors, system checks of the electrical system were ran to ensure that the system was working correctly. While completing system checks, the pressure transducers that were not functioning properly were also replaced and any issue with the wiring or signal was addressed. Once the electrical components were replaced, the restoration of the fluid system began. Parts of the tubing in the CryoCart had to be rebuild and often consisted of sizing, cutting, bending, filing, and sanding the tubing to prepare it to be flared. Many components had to be proof-tested to bring their certifications up to date, and several components had to be replaced. Various flex hoses, valves, and fittings were send to the Clean Lab because they were new, dirty, or had gone through proof-testing. Once they arrived from the cleaning lab they had to be put back to the system and leak checks and functional tests were conducted. In the Nitrogen system, the copper tubing located in the Oxygen cart was rebuild and Aerogel insulation was added to this section. A new gaseous nitrogen system was added to the CryoCart to purge the vacuum tube which will serve as the test chamber. Once the CryoCart was completed, construction of parts of the vacuum tube began. A flange was manufactured with welded fittings to hold the line of the vacuum pump as well as some extra fittings which will serve as extra inlets used to introduce fluid lines to the vacuum tube. Stress analysis was ran in this flange to ensure that it would not fail under vacuum conditions. The fluid lines leading from the air side of the vacuum to the test article were also constructed and added to the mount that had already been manufactured. Three different sets of tubing were constructed to accommodate the seven different RCS thruster and the main engine igniter that are going to be tested. Full electrical system checks were completed to ensure that all the wire harnesses and valves were functioning. Upon the completion of the CryoCart and the vacuum tube, hot fire testing for the RCS thrusters and the main engine igniter are going to begin. During this time any issues encountered with the engines or igniter will be addressed to ensure that the components function under vacuum conditions. After successful completion of testing, the vertical test bed, Morpheus, will be rebuilt and prepared to be sent to Plum Brook. In Plum Brook, the vehicle will be tested in the thermal-vacuum chamber to demonstrate that integrated lox-methane propulsion system operation in space-like conditions. This internship has allowed me the opportunity to gain valuable hands on experience and to develop skills that will aid in my education as well as in the workforce, while at the same time helping me determine that I would like to further pursue a career in propulsion engineering.

Chaidez, Mariana↗
NASA NTRS2014

Data Analysis & Statistical Methods for Command File Errors

This paper explains current work on modeling for managing the risk of command file errors. It is focused on analyzing actual data from a JPL spaceflight mission to build models for evaluating and predicting error rates as a function of several key variables. We constructed a rich dataset by considering the number of errors, the number of files radiated, including the number commands and blocks in each file, as well as subjective estimates of workload and operational novelty. We have assessed these data using different curve fitting and distribution fitting techniques, such as multiple regression analysis, and maximum likelihood estimation to see how much of the variability in the error rates can be explained with these. We have also used goodness of fit testing strategies and principal component analysis to further assess our data. Finally, we constructed a model of expected error rates based on the what these statistics bore out as critical drivers to the error rate. This model allows project management to evaluate the error rate against a theoretically expected rate as well as anticipate future error rates.

Correlation Analysis↗
NASA NTRS2016

Alternative Method of On-Orbit Response-Versus-Scan-Angle Characterization for MODIS Reflective Solar Bands

The moderate resolution imaging spectroradiometer (MODIS) has 20 reflective solar bands (RSB), covering a spectral range from 0.41 to 2.2 microns, which are calibrated on-orbit using its onboard calibrators, which include a solar diffuser, a solar diffuser stability monitor, and a spectroradiometric calibration assembly. A space view (SV) port is used to provide a background reference and also facilitates near-monthly lunar observations through a spacecraft roll. In every scan, the Earth's surface, SV, and onboard calibrators are viewed via a two-sided scan mirror, the reflectance of which depends on the angle of incidence (AOI) as well as the wavelength of the incident light. Response-versus-scan-angle (RVS) is defined as a dependence function of the scan mirror's reflectance over AOI. An initial RVS for each RSB was measured prelaunch for both Terra and Aqua MODIS. Algorithms have been developed to track the on-orbit RVS variation using the measurements from the onboard calibrators, supplemented with the earth view (EV) trends from pseudoinvariant desert targets obtained at different AOI. Since the mission beginning, the MODIS characterization support team (MCST) has dedicated efforts in evaluating approaches of characterizing the on-orbit RVS. A majority of the approaches focused on fitting the data at each AOI over time and then deriving the relative change at different AOI. The current version of the on-orbit RVS algorithm, as implemented in the collection 6 (C6) level-1B (L1B), is also based on the above rationale. It utilizes the EV response trends from the pseudoinvariant Libyan desert targets to supplement the gain derived from the onboard calibrators. The primary limitation of this approach is the assumption of the temporal stability of these desert sites. Consequently, MCST developed an approach that derives the on-orbit RVS change using measurements from a single desert site, combined with the on-orbit lunar measurements. In addition, the EV and onboard responses are fit first as a function of AOI before fitting temporally in order to eliminate the dependence on the stability of the desert site. Comprehensive comparisons are performed with current C6 RVS results for both Terra and Aqua MODIS. Results demonstrate that this alternative method provides a supplemental means to monitor the on-orbit RVS for MODIS RSB.

solar bands↗
NASA NTRS2016

Predictive Simulation of Gas Adsorption in Fixed-Beds and Limitations due to the Ill-Posed Danckwerts Boundary Condition

The 1-D axially dispersed plug flow model is a mathematical model widely used for the simulation of adsorption processes. Lumped mass transfer coefficients such as the Glueckauf linear driving force (LDF) term and the axial dispersion coefficient are generally obtained by fitting simulation results to the experimental breakthrough test data. An approach is introduced where these parameters, along with the only free parameter in the energy balance equations, are individually fit to specific test data that isolates the appropriate physics. It is shown that with this approach this model provides excellent simulation results for the C02 on zeolite SA sorbent/sorbate system; however, for the H20 on zeolite SA system, non-physical deviations from constant pattern behavior occur when fitting dispersive experimental results with a large axial dispersion coefficient. A method has also been developed that determines a priori what values of the LDF and axial dispersion terms will result in non-physical simulation results for a specific sorbent/sorbate system when using the one-dimensional axially dispersed plug flow model. A relationship between the steepness of the adsorption equilibrium isotherm as indicated by the distribution factor, the magnitude of the axial dispersion and mass transfer coefficient, and the resulting non-physical behavior is derived. This relationship is intended to provide a guide for avoiding non-physical behavior by limiting the magnitude of the axial dispersion term on the basis of the mass transfer coefficient and distribution factor.

Knox, James Clinton↗
NASA NTRS2016

Spectral and Temporal Properties of the Ultra-Luminous X-Ray Pulsar in M82 from 15 Years of Chandra Observations and Analysis of the Pulsed Emission Using NuSTAR

The recent discovery by Bachetti et al. of a pulsar in M82 that can reach luminosities of up to 10(exp 40) erg s(exp -1), a factor of approximately 100 times the Eddington luminosity for a 1.4 solar mass compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.58 keV X-ray emission from this source from 15 years of Chandra observations. We analyze 19 ACIS observations where the point-spread function (PSF) of the pulsar is not contaminated by nearby sources. We fit the Chandra spectra of the pulsar with a power-law model and a disk blackbody model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where four observations have a pile-up fraction of 10, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (LX greater than 10(exp 39) erg s(exp -1) is equal to gamma 1.33 +/-.0.15. For the disk blackbody model, the average temperature is T(sub in) 3.24 +/- 0.65 keV, the spectral shape being consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cutoff, where gamma is equal to 0.6 +/- 0.3 and E(sub C) is equal to 14(exp +5) (sub -3)) keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10( exp (39) erg s (exp −1), greater than 10 times its Eddington limit.

individual (M82) - stars: neutron X rays; binaries↗
NASA NTRS2016

Analysis of One-Way Laser Ranging Data to LRO, Time Transfer and Clock Characterization

We processed and analyzed one-way laser ranging data from International Laser Ranging Service ground stations to NASA's Lunar Reconnaissance Orbiter (LRO), obtained from June 13, 2009 until September 30, 2014. We pair and analyze the one-way range observables from station laser fire and spacecraft laser arrival times by using nominal LRO orbit models based on the GRAIL gravity field. We apply corrections for instrument range walk, as well as for atmospheric and relativistic effects. In total we derived a tracking data volume of approximately 3000 hours featuring 64 million Full Rate and 1.5 million Normal Point observations. From a statistical analysis of the dataset we evaluate the experiment and the ground station performance. We observe a laser ranging measurement precision of 12.3 centimeters in case of the Full Rate data which surpasses the LOLA (Lunar Orbiting Laser Altimeter) timestamp precision of 15 centimeters. The averaging to Normal Point data further reduces the measurement precision to 5.6 centimeters. We characterized the LRO clock with fits throughout the mission time and estimated the rate to 6.9 times10 (sup -8), the aging to 1.6 times 10 (sup -12) per day and the change of aging to 2.3 times 10 (sup -14) per day squared over all mission phases. The fits also provide referencing of onboard time to the TDB (Barycentric Dynamical Time) time scale at a precision of 166 nanoseconds over two and 256 nanoseconds over all mission phases, representing ground to space time transfer. Furthermore we measure ground station clock differences from the fits as well as from simultaneous passes which we use for ground to ground time transfer from common view observations. We observed relative offsets ranging from 33 to 560 nanoseconds and relative rates ranging from 2 times 10 (sup -13) to 6 times 10 (sup -12) between the ground station clocks during selected mission phases. We study the results from the different methods and discuss their applicability for time transfer.

Bauer, S.↗
NASA NTRS2017

New-Generation NASA Aura Ozone Monitoring Instrument (OMI) Volcanic SO2 Dataset: Algorithm Description, Initial Results, and Continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS)

Since the fall of 2004, the Ozone Monitoring Instrument (OMI) has been providing global monitoring of volcanic SO2 emissions, helping to understand their climate impacts and to mitigate aviation hazards. Here we introduce a new-generation OMI volcanic SO2 dataset based on a principal component analysis (PCA) retrieval technique. To reduce retrieval noise and artifacts as seen in the current operational linear fit (LF) algorithm, the new algorithm, OMSO2VOLCANO, uses characteristic features extracted directly from OMI radiances in the spectral fitting, thereby helping to minimize interferences from various geophysical processes (e.g., O3 absorption) and measurement details (e.g., wavelength shift). To solve the problem of low bias for large SO2 total columns in the LF product, the OMSO2VOLCANO algorithm employs a table lookup approach to estimate SO2 Jacobians (i.e., the instrument sensitivity to a perturbation in the SO2 column amount) and iteratively adjusts the spectral fitting window to exclude shorter wavelengths where the SO2 absorption signals are saturated. To first order, the effects of clouds and aerosols are accounted for using a simple Lambertian equivalent reflectivity approach. As with the LF algorithm, OMSO2VOLCANO provides total column retrievals based on a set of predefined SO2 profiles from the lower troposphere to the lower stratosphere, including a new profile peaked at 13 km for plumes in the upper troposphere. Examples given in this study indicate that the new dataset shows significant improvement over the LF product, with at least 50% reduction in retrieval noise over the remote Pacific. For large eruptions such as Kasatochi in 2008 (approximately 1700 kt total SO2/ and Sierra Negra in 2005 (greater than 1100DU maximum SO2), OMSO2VOLCANO generally agrees well with other algorithms that also utilize the full spectral content of satellite measurements, while the LF algorithm tends to underestimate SO2. We also demonstrate that, despite the coarser spatial and spectral resolution of the Suomi National Polar-orbiting Partnership (Suomi-NPP) Ozone Mapping and Profiler Suite (OMPS) instrument, application of the new PCA algorithm to OMPS data produces highly consistent retrievals between OMI and OMPS. The new PCA algorithm is therefore capable of continuing the volcanic SO2 data record well into the future using current and future hyperspectral UV satellite instruments.

OMI↗
NASA NTRS2016

Polarimetry and Flux Distribution in the Debris Disk Around HD 32297

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at > or = 5(sigma) levels from approx. 0.45" to approx.1.7" (50-192AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of approx. 0.75" (NE side) and approx. 0.65" (SW side). Global forward-modeling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from approx. 0.25-1.6", although the central region is quite noisy and high S/N are only found in the range approx. 0.75-1.2". The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from approx. 10% at 0.55" to approx. 25% at 1.6". The maximum is found at scattering angles of 90, either from the main components of the disk or from dust grains blown out to larger radii.

high-contrast angular differential imaging (ADI)↗
NASA NTRS2017

HST WFC3 Observations of Uranus' 2014 Storm Clouds and Comparison with VLT/SINFONI and IRTF/Spex Observations

In November 2014 Uranus was observed with the Wide Field Camera 3 (WFC3) instrument of the Hubble Space Telescope as part of the Hubble 2020: Outer Planet Atmospheres Legacy program, OPAL. OPAL annually maps Jupiter, Uranus and Neptune (and will also map Saturn from 2018) in several visible near- infrared wavelength filters. The Uranus 2014 OPAL observations were made on the 89th November at a time when a huge cloud complex, first observed by de Pater et al. (2015) and subsequently tracked by professional and amateur astronomers (Sayanagi et al., 2016), was present at 30-40deg N. We imaged the entire visible atmosphere, including the storm system, in seven filters spanning 467924 nm, capturing variations in the coloration of Uranus clouds and also vertical distribution due to wavelength dependent changes in Rayleigh scattering and methane absorption optical depth. Here we analyse these new HST observations with the NEMESIS radiative-transfer and retrieval code in multiple-scattering mode to determine the vertical cloud structure in and around the storm cloud system. The same storm system was also observed in the H-band (1.4-1.8 micrometers) with the SINFONI Integral Field Unit Spectrometer on the Very Large Telescope (VLT) on 31st October and 11th November, reported by Irwin et al. (2016, 10.1016j.icarus.2015.09.010). To constrain better the cloud particle sizes and scattering properties over a wide wavelength range we also conducted a limb-darkening analysis of the background cloud structure in the 30-40deg N latitude band by simultaneously fitting: a) these HSTOPAL observations at a range of zenith angles; b) the VLTSINFONI observations at a range of zenith angles; and c) IRTFSpeX observations of this latitude band made in 2009 at a single zenith angle of 23deg, spanning the wavelength range 0.8-1.8 micrometers (Irwin et al., 2015, 10.1016j.icarus.2014.12.020). We find that the HST observations, and the combined HSTVLTIRTF observations at all locations are well modelled with a three-component cloud comprised of: 1) a vertically thin, but optically thick deep tropospheric cloud at a pressure of approximately 2 bars; 2) a methane-ice cloud based at the methane-condensation level of approximately 1.23 bar, with variable vertical extent; and 3) a vertically extended tropospheric haze, also based at the methane-condensation level of 1.23 bar. We find that modelling both haze and tropospheric cloud with particles having an effective radius of approximately 0.1 micron provides a good fit the observations, although for the tropospheric cloud, particles with an effective radius as large as 1.0 micron provide a similarly good fit. We find that the particles in both the tropospheric cloud and haze are more scattering at short wave- lengths, giving them a blue color, but are more absorbing at longer wavelengths, especially for the tropospheric haze. We find that the spectra of the storm clouds are well modelled by localized thickening and vertical extension of the methane-ice cloud. For the particles in the storm clouds, which we assume to be composed of methane ice particles, we find that their mean radii must lie somewhere in the range 0. 1 1. 0 m. We find that the high clouds have low integrated opacity, and that streamers reminiscent of convective thunderstorm anvils are confined to levels deeper than 1 bar. These results argue against vigorous moist convective origins for the cloud features.

atmospheres↗
NASA NTRS2016

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗