Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error bound”

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 703 records · Page 39

On the design of turbo codes

In this article, we design new turbo codes that can achieve near-Shannon-limit performance. The design criterion for random interleavers is based on maximizing the effective free distance of the turbo code, i.e., the minimum output weight of codewords due to weight-2 input sequences. An upper bound on the effective free distance of a turbo code is derived. This upper bound can be achieved if the feedback connection of convolutional codes uses primitive polynomials. We review multiple turbo codes (parallel concatenation of q convolutional codes), which increase the so-called 'interleaving gain' as q and the interleaver size increase, and a suitable decoder structure derived from an approximation to the maximum a posteriori probability decision rule. We develop new rate 1/3, 2/3, 3/4, and 4/5 constituent codes to be used in the turbo encoder structure. These codes, for from 2 to 32 states, are designed by using primitive polynomials. The resulting turbo codes have rates b/n (b = 1, 2, 3, 4 and n = 2, 3, 4, 5, 6), and include random interleavers for better asymptotic performance. These codes are suitable for deep-space communications with low throughput and for near-Earth communications where high throughput is desirable. The performance of these codes is within 1 dB of the Shannon limit at a bit-error rate of 10(exp -6) for throughputs from 1/15 up to 4 bits/s/Hz.

Divsalar, D.↗

Early Information Parameter-Set Analysis for Satellite Close Approaches using Machine Learning

Understanding orbital mechanics is essential in space flight and navigation applications, and leveraging modern force models for flight path projection remains an important aspect in space mission design and operation. However, force models do not capture all the dynamics or perturbations in the space environment and thus are subject to errors in predicting the state vectors. The further out the predicted miss distance between spacecraft is from the time of closest approach (TCA), the larger the propagated errors in the predicted miss distance at TCA is. The dependency on these force models for spacecraft flight state prediction calls for a more reliable method that can quantify, or even reduce, these propagated errors. With recent advances in the field artificial intelligence, specifically in machine and deep learning algorithms, a model that implements these approaches can improve on the modern force model approach. The goal for this work is to provide an early-information decision-making threshold, in order to prioritize risk assessment implementation, given the ongoing increase of space objects. In analyzing the relationship of several parameters from conjunction data messages(CDMs) and solar information, early information becomes viable in miss distance prediction with unsupervised learning techniques, which learn the parameters that are linked together with miss distance and probability of collision (Pc) variables. Another approach implemented for identifying relationships within CDMs is supervised learning, in which a shallow neural network binary classifier learns to distinguish events with Pc values¡108. These parameters detected in the unsupervised process are then applied to a regression neural network, which predicts the miss distance at TCA for a specific event within a given uncertainty bound. For the regression neural network, a Long Short Term Memory (LSTM) neural network is implemented, which yields memory about each time step in an event. Using an LSTM network, the model learns to predict miss distance within 0.2km of the value measured at TCA. Although there is a limited amount of "close miss" data to train a network, the network learns to associate parameters, like large energy dissipation rates with the secondary object, with an elevated Pc

Brianna I. Robertson↗

Deflections of anisotropic sandwich beams with variable face sheets and core thicknesses

A sandwich construction consists of a low-density core material with high strength face sheets bounded to the top and bottom surfaces. The construction has been widely used in the aerospace and marine industries due to its outstanding characteristics such as noise absorption, weight minimization, heat insulation, and better bending stiffness. In sandwich structures used in high-performance aircraft, the face sheets are often made of fiber-reinforced composite materials and the core is made of honeycomb. The structures may also have variable thickness so as to satisfy aerodynamic requirements. In the stress analysis, the constant-thickness face sheets are usually considered as membrane and the core is assumed to be inextensible but deformable in the thickness direction. The static behavior of variable-thickness, isotropic and homogeneous sandwich beams was successfully studied by employing a constant-thickness theory but allowing stiffnesses to vary in accordance with local thickness variations. It has been recently found in a refined theory that the analyses based on the constant thickness theory locally can lead to significant errors in structural responses if the sandwich beam is thickness-tapered and the cores are deformable in transverse shear. The errors arise mainly from two factors: (1) the transverse shear components of the membrane forces in the face sheets alter the transverse shears carried by the core; and (2) the face-sheet membrane strains arise from transverse shear deformation of the core. In practice the variable thickness may not only exist in core but also in face sheets. The thickness-variations may even be a type of step function. In this case the transverse shear stress in the face sheets and bending stress in the core should be taken into account in the refined theory mentioned. In the present study, energy principles are employed in deriving governing equations for general bending of anisotropic sandwich beams with variable thickness in both face sheets and cores. Solutions to these equations are based on a finite difference scheme. As an example in application, a simply supported thickness-tapered sandwich beam subject to a concentrated load at its center is considered. Let W' be the maximum deflection of the beam in which face sheets are considered as membrane, while W'' is that based on using the modified refined theory. It is found that W' is always larger than W'', however, the magnitude of (W'- W'') appears to be insensitive to the change of the taper of the beam.

Lu, Chu-Ho↗

Determination of convergence rates across the Ventura Basin, Southern California, using GPS and historical triangulation

Comparison of angles from historical triangulation observations dating as far back as 1932 with Global Positions System (GPS) measurements taken in 1987 indicates that rapid convergence may be taking place on decade timescales in the central and eastern part of the Ventura basin, an east-west trending trough bounded by thrust faults. Changes in angles over this time were analyzed using Prescott's modified Frank's method and in terms of a model which assumes that the regions to the north and south of the basin are rigid blocks undergoing relative motion. For the two block model, inversion of the observed angle changes over the last 28 years for the relative motion vector leads to north-south convergence across the basin of 30 + or - 5 mm/yr, with a left lateral component of 10 + or - 1 mm/yr in the Fillmore-Santa Paula area in the central part of the basin. The modified Frank's method yields strain rates of approximately 2 microrad/yr in both the east and central parts of the basin for measurements spanning the 1971 San Fernando earthquake. Assuming no east-west strain yeilds north-south compression of approximately 3.5 + or - .2 cm/yr. Comparison of triangulation data prior to the earthquake shows no strain outside the margin of error. The convergence rates determined by geodetic techniques are consistent with geologic observations in the area. Such large geodetic deformation rates, with no apparent near-surface creep on the major thrust, can be understood if these faults become subhorizontal at relatively shallow depths and if the subhorizontal portions of the faults are creeping. An alternative explanation of the large displacement rates might be that the pumping of oil in the vicinity of the benchmarks caused large horizontal motions, although it is unlikely that meter scale horizontal motions are due to oil withdrawal. These and other hypotheses are evaluated to better constrain the tectonics of this active region.

Donnellan, Andrea↗

New Experiment to Measure the Electron Electric Dipole Moment

An electron can possess an electric dipole moment (edm) only if time reversal symmetry (T) is violated. No edm of any particle has yet been discovered. CP-violation, equivalent to T-violation by the CPT theorem, does occur in Kaon decays and can be accounted for by the standard model. However, this mechanism leads to an electron edm d(sub e) of the order of 10(exp -38) e cm, whereas the current experimental bound on d(sub e) is about 10(exp -27) e cm. However, well-motivated extensions of the standard model such as supersymmetric theories do predict that de could be as large as the current bound. In addition, CP violation in the early universe is required to explain the preponderance of matter over anti-matter, but the exact mechanism of this CP violation is unclear. For these reasons, we are undertaking a new experimental program to determine de to an improved accuracy of 10(exp -29) e cm. Our experiment will use laser-cooled, trapped Cesium atoms to measure the atomic edm d(sub Cs) that occurs if d(sub e) is not zero. In order to do this, we will measure the energy splitting between the atoms spin states in parallel electric and magnetic fields. The signature of an edm would be a linear dependence of the splitting on the electric field E due to the interaction - d(sub Cs) dot E. Our measurement will be much more sensitive than previous measurements because atoms can be stored in the trap for tens of seconds, allowing for much narrower Zeeman resonance linewidths. Also, our method eliminates the most important systematic errors, proportional to atomic velocity, which have limited previous experiments. In this presentation, we will describe the design of our new apparatus, which is presently under construction. An important feature of our experimental apparatus is that magnetic field noise will be suppressed to a very low value of the order of 1 fT/(Hz)1/2. This requires careful attention to the Johnson noise currents in the chamber, which have not been important in previous experiments. In addition we will present estimates of the limits of the various errors that we expect for our experiment.

Kittle, Melanie↗

On complexity of trellis structure of linear block codes

The trellis structure of linear block codes (LBCs) is discussed. The state and branch complexities of a trellis diagram (TD) for a LBC is investigated. The TD with the minimum number of states is said to be minimal. The branch complexity of a minimal TD for a LBC is expressed in terms of the dimensions of specific subcodes of the given code. Then upper and lower bounds are derived on the number of states of a minimal TD for a LBC, and it is shown that a cyclic (or shortened cyclic) code is the worst in terms of the state complexity among the LBCs of the same length and dimension. Furthermore, it is shown that the structural complexity of a minimal TD for a LBC depends on the order of its bit positions. This fact suggests that an appropriate permutation of the bit positions of a code may result in an equivalent code with a much simpler minimal TD. Boolean polynomial representation of codewords of a LBC is also considered. This representation helps in study of the trellis structure of the code. Boolean polynomial representation of a code is applied to construct its minimal TD. Particularly, the construction of minimal trellises for Reed-Muller codes and the extended and permuted binary primitive BCH codes which contain Reed-Muller as subcodes is emphasized. Finally, the structural complexity of minimal trellises for the extended and permuted, and double-error-correcting BCH codes is analyzed and presented. It is shown that these codes have relatively simple trellis structure and hence can be decoded with the Viterbi decoding algorithm.

Lin, Shu↗

Pressure Analysis for LAVA-OVEN

The Lunar Advanced Volatiles Analysis (LAVA) and the Oxygen Volatiles Extraction Node (OVEN) are subsystems included in the Regolith Environment Science, and Oxygen Lunar Volatiles Extraction (RESOLVE) payload bound for the Moon in 2019. This Resource Prospector Mission (RPM) has the objective of landing on a shadowed region of the Moons South Pole to collect data and determine whether the resources could be effectively used for space exploration systems. The quantification of the resources will help understand if it can adequately minimize materials carried from Earth by: providing life support, propellants, construction materials or energy supply to the payload or crew. This paper outlines the procedures done for the pressure analysis of the LAVA-OVEN (LOVEN) Integration Testing. The pressure analysis quantifies how much gases and water are present in the sample tested during the Engineering Testing Unit (ETU) phase of instrument development. Ultimately the purpose of these tests is to improve the estimate of the amount of water in each Lunar sample and reduce the time necessary for this estimate. The governing principle that was used for the analysis is the Ideal Gas Law, PV=nRT where P stands for pressure, V for volume, n for number of moles, R being the gas constant and T for temperature. We also estimate the errors involved in these measured and derived quantities since a key objective of the mission is to estimate the quantity of volatiles present in the lunar samples introduced into OVEN.

Lunar↗

The Utilization Profiles of the CCSDS Unified Space Link Protocol (USLP)

The purpose of this paper is to identify the utilization profiles for interfacing the Data Protocol Sublayer using the Unified Space Link Protocols (USLP) (reference 1) with the space link coding procedures as specified in the CCSDS Coding & Synchronization Blue Books (references 2 through 5), used in both telecommand and telemetry applications. This paper describes how the USLP Protocol utilizes the coding and synchronization sublayer to support: a. Direct to Earth (DTE) telemetry links for engineering and science data b. Direct to Earth (DTE) telemetry links for very high rate science data c. Direct from Earth (DFE) command, sequencing and flight software loads d. Space to Space Links (Proximity) utilized by orbiters for data exchange to/from surface bound assets. The CCSDS has divided the functions of the Data Link Layer into two sublayers: the Data Link Protocol Sublayer (DLP-SL) and the Coding and Synchronization Sublayer (CS-SL). The Data Link Protocol Sublayer (DLP-SL) interfaces to the users, accepting the data that is to be transported, on the sending side of the link, and delivering that data on the receiving end. The Transfer Frame is the data unit that is transferred across the Data Link Protocol Sublayer and the Coding and Synchronization Sublayer boundary. The Coding and Synchronization Sublayer (CS-SL) provides the encoding, randomization, and frame synchronization functions that prepares the USLP Transfer Frame for transport across the space link. The CS-SL is divided into 2 processes: 1) The Frame Interface Processes (FIP) performs the interface functions required to prepare the data for delivery to the Coding/Decoding Process (CDP). This process includes prepending a Frame Start Marker to the provided frame, when management has designated that the frame is not to be aligned to the codeblock or when there is no block code used. 2) The Coding/Decoding Process (CDP) performs the forward error correction processes that are used to optimize the performance of the link and minimize the error rate. The CDP creates the symbol stream that is delivered to the Physical Layer. The transfer of the USLP transfer frames across different types of space links is the focus of this paper. The Protocol Data Unit (PDU) that is passed in both directions between the Data Link Protocol Sublayer (DLP-SL) and Coding and Synchronization Sublayer (CS-SL) is the transfer frame. The USLP frame structure provides flexibility that can be constrained by the functions utilized within the CS-SL that prepare the transfer frame for transit. For example, the USLP transfer frame contains a length field that enables the frame to be of variable length but CS-SL under certain conditions may constrain the frame to be fixed in length. This paper describes 5 operational modes available for use by the Data Link Layer to provide data exchange across the USLP space link. These modes are different because different operational requirements apply to vastly different types of space links and thus the communications implementation requirements differ. The environmental issues include the power or energy available, the distance between the end points of the link, the complexity of the equipment available at those end points, the atmospheric conditions and radiometric frequency selection. The CS-SL utilizes different forward error correcting codes supported by specific operational modes to configure the data for transit. This paper describes all of the operational modes in a series of data models which decompose the functionality between the Data Link Protocol Sublayer and the Coding and Synchronization sublayer. The operational modes described are: 1. Uncoded Mode: has been used for short links that contain significant available power to provide an acceptable frame error rate. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 2. Convolutional Only Mode: is currently the prime forward error correction coding used for the proximity links. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 3. Variable Length Frame Aligned to Variable Length Codeblock (TC): is used for Direct from Earth links were power levels are high and the simple, least complex code i.e., the BCH code is used. This mode has been in use since the early 1970s. The BCH code is a short code and the decoder is easy to implement. 4. Fixed Length Frame Aligned to Fixed Length Codeblock (AOS/TM): was introduced when the concatenated Convolutional and Reed-Solomon Code was formulated to provide significant reduction in link data error rate and the ability to determine if there was an error in the decoded codeblock. The frame is aligned to the codeblock so that there is a one to one relationship of frame errors to codeblock errors without additional error detection coding being added. This mode requires the protocol frames to be the exact size of the message portion of the codeblock. 5. Frames Unaligned to Fixed Length Codeblocks (Currently used for very high rates and space to space links): This mode is currently used for missions that have a very high data rate that can be controlled adaptively as the environment changes and as the next generation operating mode for the proximity link. This mode from a coded data stream point of view is exactly like that described in 4. above, except that the frame need not be aligned to the codeblock. There is no requirement on frame length when using this mode. Thus when using USLP it can be used to support links that require short or long frames. There is also no mandatory requirement that frames cannot be separated by idle data reducing the tight data rate connection requirements between the data link protocol sublayer and the coding & synchronization sublayer. In conclusion, how these operational modes can be put to use in mission operational scenarios is described for Direct from Earth links (DFE), Direct to Earth links (DTE), and Proximity links.

Greenberg, E.↗

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↗

Phase-Controlled Magnetic Mirror for Wavefront Correction

Typically, light interacts with matter via the electric field and interaction with weakly bound electrons. In a magnetic mirror, a patterned nanowire is fabricated over a metallic layer with a dielectric layer in between. Oscillation of the electrons in the nanowires in response to the magnetic field of incident photons causes a re-emission of photons and operation as a "magnetic mirror." By controlling the index of refraction in the dielectric layer using a local applied voltage, the phase of the emitted radiation can be controlled. This allows electrical modification of the reflected wavefront, resulting in a deformable mirror that can be used for wavefront control. Certain applications require wavefront quality in the few-nanometer regime, which is a major challenge for optical fabrication and alignment of mirrors or lenses. The use of a deformable magnetic mirror allows for a device with no moving parts that can modify the phase of incident light over many spatial scales, potentially with higher resolution than current approaches. Current deformable mirrors modify the incident wavefront by using nano-actuation of a substrate to physically bend the mirror to a desired shape. The purpose of the innovation is to modify the incident wavefront for the purpose of correction of fabrication and alignment-induced wavefront errors at the system level. The advanced degree of precision required for some applications such as gravity wave detection (LISA - Laser Interferometer Space Antenna) or planet finding (FKSI - Fourier-Kelvin Stellar Interferometer) requires wavefront control at the limits of the current state of the art. All the steps required to fabricate a magnetic mirror have been demonstrated. The modification is to apply a bias voltage to the dielectric layer so as to change the index of refraction and modify the phase of the reflected radiation. Light is reflected off the device and collected by a phase-sensing interferometer. The interferometer determines the initial wavefront of the device and fore optics. A wavefront correction is calculated, and voltage profile for each nanowire strip is determined. The voltage is applied, modifying the local index of refraction of the dielectric under the nanowire strip. This modifies the phase of the reflected light to allow wavefront correction.

Hagopian, John↗

A General Closed-Form Solution for the Lunar Reconnaissance Orbiter (LRO) Antenna Pointing System

The National Aeronautics and Space Administration s (NASA) Lunar Reconnaissance Orbiter (LRO) launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle into a direct insertion trajectory to the Moon LRO, designed, built, and operated by the NASA Goddard Space Flight Center in Greenbelt, MD, is gathering crucial data on the lunar environment that will help astronauts prepare for long-duration lunar expeditions. During the mission s nominal life of one year its six instruments and one technology demonstrator will find safe landing site, locate potential resources, characterize the radiation environment and test new technology. To date, LRO has been operating well within the bounds of its requirements and has been collecting excellent science data images taken from the LRO Camera Narrow Angle Camera (LROC NAC) of the Apollo landing sites have appeared on cable news networks. A significant amount of information on LRO s science instruments is provided at the LRO mission webpage. LRO s Attitude Control System (ACS), in addition to controlling the orientation of the spacecraft is also responsible for pointing the High Gain Antenna (HGA). A dual-axis (or double-gimbaled) antenna, deployed on a meter-long boom, is required to point at a selected Earth ground station. Due to signal loss over the distance from the Moon to Earth, pointing precision for the antenna system is very tight. Since the HGA has to be deployed in spaceflight, its exact geometry relative to the spacecraft body is uncertain. In addition, thermal distortions and mechanical errors/tolerances must be characterized and removed to realize the greatest gain from the antenna system. These reasons necessitate the need for an in-flight calibration. Once in orbit around the moon, a series of attitude maneuvers was conducted to provide data needed to determine optimal parameters to load onboard, which would account for the environmental and mechanical errors at any antenna orientation. The nominal geometry for the HGA involves an outer gimbal axis that is exactly perpendicular to the inner gimbal axis, and a target direction that is exactly perpendicular to the outer gimbal axis. For this nominal geometry, closed-form solutions of the desired gimbal angles are simple to get for a desired target direction specified in the spacecraft body fame. If the gimbal axes and the antenna boresight are slightly misaligned, the nominal closed-form solution is not sufficiently accurate for computing the gimbal angles needed to point at a target. In this situation, either a general closed-form solution has to be developed for a mechanism with general geometries, or a correction scheme has to be applied to the nominal closed-form solutions. The latter has been adopted for Solar Dynamics Observatory (SDO) as can be seen in Reference 1, and the former has been used for LRO. The advantage of the general closed-form solution is the use of a small number of parameters for the correction of nominal solutions, especially in the regions near singularities. Singularities here refer to cases when the nominal closed-form solutions have two or more solutions. Algorithm complexity, however, is the disadvantage of the general closed-form solution.

Shah, Neerav↗

Incorporating Endmember Variability into Spectral Mixture Analysis Through Endmember Bundles

Variation in canopy structure and biochemistry induces a concomitant variation in the top-of-canopy spectral reflectance of a vegetation type. Hence, the use of a single endmember spectrum to track the fractional abundance of a given vegetation cover in a hyperspectral image may result in fractions with considerable error. One solution to the problem of endmember variability is to increase the number of endmembers used in a spectral mixture analysis of the image. For example, there could be several tree endmembers in the analysis because of differences in leaf area index (LAI) and multiple scatterings between leaves and stems. However, it is often difficult in terms of computer or human interaction time to select more than six or seven endmembers and any non-removable noise, as well as the number of uncorrelated bands in the image, limits the number of endmembers that can be discriminated. Moreover, as endmembers proliferate, their interpretation becomes increasingly difficult and often applications simply need the aerial fractions of a few land cover components which comprise most of the scene. In order to incorporate endmember variability into spectral mixture analysis, we propose representing a landscape component type not with one endmember spectrum but with a set or bundle of spectra, each of which is feasible as the spectrum of an instance of the component (e.g., in the case of a tree component, each spectrum could reasonably be the spectral reflectance of a tree canopy). These endmember bundles can be used with nonlinear optimization algorithms to find upper and lower bounds on endmember fractions. This approach to endmember variability naturally evolved from previous work in deriving endmembers from the data itself by fitting a triangle, tetrahedron or, more generally, a simplex to the data cloud reduced in dimension by a principal component analysis. Conceptually, endmember variability could make it difficult to find a simplex that both surrounds the data cloud and has vertices that are realistic endmember spectra with reflectances between 0 and 1. In this paper, we create endmember bundles and bounding fraction images for an AVIRIS subscene simulated with a plant canopy radiative transfer model. The simulated subscene is spatially patterned after a subscene from the AVIRIS image acquired August, 1993 over La Copita, Texas. In addition, for comparison, we performed a traditional unmixing with image endmembers.

Bateson, C. Ann↗

Chandra Observations of the X-ray Environs of SN 1998bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the GRB. Eight X-ray point sources were localized, three and five each in the original error boxes - S1 and S2 - assigned for variable X-ray counterparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by CXO on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared, and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S1, Sla, is coincident with the well-determined radio location of SN 1998bw, and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning approximately 1300 days, is obtained by assuming the burst and supernova were coincident at 35.6 Mpc. When this X-ray iight curve is compared with those of the X-ray afterglows of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence emerges for at least two classes of lightcurves, perhaps bounding a continuum. By three to ten years, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998bw and GRB 980425 took place in the same object. One possible explanation for the two classes is a (nearly) standard GRB observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRBs to the ULX source population.

Kouveliotou, C.↗

Chandra Observations of the X-ray Environs of SN 1998 bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the GRB. Eight X-ray point sources were localized, three and five each in the original error boxes--S1 and S2--assigned for variable X-ray counterparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by CXO on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared, and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S1, S1a, is coincident with the well-determined radio location of SN 1998bw, and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning (approx) 1300 days, is obtained by assuming the burst and supernova were coincident at 35. 6 Mpc. When this X-ray light curve is compared with those of the X-ray 'afterglows' of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence emerges for at least two classes of lightcurves, perhaps bounding a continuum.

X-RAY SOURCES↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

Chandra Observations of the X-ray Environs of SN 1998bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the gamma-ray burst (GRB). Eight X-ray point sources were localized, three and five each in the original error boxes, S1 and S2, assigned for variable X-ray counteparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by Chandra on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S 1, S 1 a, is coincident with the well-determined radio location of SN 1998bw and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning approx. 1400 days, is obtained by assuming that the burst and supernova were coincident at 35.6 Mpc. When this X-ray light curve is compared with those of the X-ray af "erglows" of ordinary GRBs, X-ray flashes, and ordinary supernovae, evidence emerges for at least two classes of light curves, perhaps bounding a continuum. By 3-10 yr, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998 bw aid GRB 980425 took place in the same object.One possible explanation for the two classes is that a (nearly) standard GRB was observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRB afterglows to the ultraluminous X-ray source population.

Kouveliotou, C.↗

Application of X-Ray Pulsar Navigation: A Characterization of the Earth Orbit Trade Space

The potential for pulsars as a navigation source has been studied since their discovery in 1967. X-ray pulsar navigation (XNAV) is a celestial navigation system that uses the consistent timing nature of x-ray photons from milli-second pulsars (MSP) to perform space navigation. By comparing the detected arrival of x-ray photons to a reference database of expected pulsar lightcurve timing models, one can infer a range and range rate measurement based on light time delay. Much of the challenge of XNAV comes from the faint signal, availability, and distant nature of pulsars. This is a study of potential pulsar XNAV measurements to measure extended Kalman filter (EKF) tracking performance with a wide trade space of bounded Earth orbits, using a simulation of existing x-ray detector space hardware. An example of an x-ray detector for XNAV is the NASA Station Explorer for X-ray Timing and Navigation (SEXTANT) mission, a technology demonstration of XNAV set to perform on the International Space Station (ISS) in late 2016early 2017. XNAV hardware implementation is driven by trajectory and environmental influences which add noise to the x-ray pulse signal. In a closed Earth orbit, the radiation environment can exponentially increase the signal noise from x-ray pulsar sources, decreasing the quality and frequency of measurements. The SEXTANT mission in particular improves on the signal to noise ratio by focusing an array of 56 x-ray silicon drift detectors at one pulsar target at a time. This reduces timing glitches and other timing noise contributions from ambient x-ray sources to within a 100 nanosecond resolution. This study also considers the SEXTANT scheduling challenges inherent in a single target observation. Finally, as the navigation sources are now relatively inertial targets, XNAV measurements are also subject to periods of occultation from various celestial bodies. This study focuses on the characterization of these drivers in closed Earth orbits and is not a tuning analysis of the EKF. The study shows that the closed Earth orbit for XNAV performance is reliant on the orbit semi-major axis and eccentricity as well as orbit inclination. These parameters are the primary drivers of pulsar measurement availability and significantly influence the natural spacecraft orbit dynamics. Sensitivity to initial orbit determination error growth due to the scarcity of XNAV measurements within an orbital period require appropriate timing of initial XNAV measurements. The orbit angles of argument of perigee and right ascension of the ascending node, alongside the other orbit parameters, complete the initial cadence of XNAV measurements. The performance of initial XNAV measurements then propagates throughout the experimental period. The study provides a basis to missions who wish to consider XNAV as a potential navigation source in a closed Earth orbit design. It provides a family of orbit trajectories as well as other modeling considerations needed to effectively evaluate if XNAV is an effective navigation source for a potential mission. As an EKF is sensitive to a linearized estimated state, this study has a direct benefit of providing effective XNAV measurements to maintain spacecraft tracking, independent of other navigation sources. In the particular use case of the SEXTANT mission, it also provides a novel scheduling algorithm which addresses the need to prioritize and manage pulsar observations for effective navigation.

PULSAR↗

Application of X-Ray Pulsar Navigation: A Characterization of the Earth Orbit Trade Space

The potential for pulsars as a navigation source has been studied since their discovery in 1967. X-ray pulsar navigation (XNAV) is a celestial navigation system that uses the consistent timing nature of x-ray photons from millisecond pulsars (MSP) to perform space navigation. By comparing the detected arrival of x-ray photons to a reference database of expected pulsar light-curve timing models, one can infer a range and range rate measurement based on light time delay. Much of the challenge of XNAV comes from the faint signal, availability, and distant nature of pulsars. This is a study of potential pulsar XNAV measurements to measure extended Kalman filter (EKF) tracking performance with a wide trade space of bounded Earth orbits, using a simulation of existing x-ray detector space hardware. An example of an x-ray detector for XNAV is the NASA Station Explorer for X-ray Timing and Navigation (SEXTANT) mission, a technology demonstration of XNAV set to perform on the International Space Station (ISS) in late 2016early 2017. XNAV hardware implementation is driven by trajectory and environmental influences which add noise to the x-ray pulse signal. In a closed Earth orbit, the radiation environment can exponentially increase the signal noise from x-ray pulsar sources, decreasing the quality and frequency of measurements. The SEXTANT mission in particular improves on the signal to noise ratio by focusing an array of 56 x-ray silicon drift detectors at one pulsar target at a time. This reduces timing glitches and other timing noise contributions from ambient x-ray sources to within a 100 nanosecond resolution. This study also considers the SEXTANT scheduling challenges inherent in a single target observation. Finally, as the navigation sources are now relatively inertial targets, XNAV measurements are also subject to periods of occultation from various celestial bodies. This study focuses on the characterization of these drivers in closed Earth orbits and is not a tuning analysis of the EKF. The study shows that the closed Earth orbit for XNAV performance is reliant on the orbit semi-major axis and eccentricity as well as orbit inclination. These parameters are the primary drivers of pulsar measurement availability and significantly influence the natural spacecraft orbit dynamics. Sensitivity to initial orbit determination error growth due to the scarcity of XNAV measurements within an orbital period require appropriate timing of initial XNAV measurements. The orbit angles of argument of perigee and right ascension of the ascending node, alongside the other orbit parameters, complete the initial cadence of XNAV measurements. The performance of initial XNAV measurements then propagates throughout the experimental period. The study provides a basis to missions who wish to consider XNAV as a potential navigation source in a closed Earth orbit design. It provides a family of orbit trajectories as well as other modeling considerations needed to effectively evaluate if XNAV is an effective navigation source for a potential mission. As an EKF is sensitive to a linearized estimated state, this study has a direct benefit of providing effective XNAV measurements to maintain spacecraft tracking, independent of other navigation sources. In the particular use case of the SEXTANT mission, it also provides a novel scheduling algorithm which addresses the need to prioritize and manage pulsar observations for effective navigation.

x-ray↗