Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error propagation”

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.

790 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral shortwave instruments produce high-resolution and information-dense spectra, thus offering the opportunity to reduce uncertainties in retrieval products. We are in the process of developing a retrieval that uses the very fast Principal Component Radiative Transfer Model in the solar spectral region (PCRTM-Solar) in the forward model calculations. This retrieval will use measured reflectances from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder) instruments. In this manuscript we present progress towards a reference retrieval employing a widely used, verified, accurate, yet computationally slower radiative transfer modeling technique. The reference retrieval, while too slow for using the complete hyperspectral measurement, will allow us to study the behavior of retrieval products and help us verify results from our in-development fast retrieval. Both retrievals will use the optimal estimation retrieval framework. In this manuscript we present results from an uncertainty analysis considering three uncertainty sources for a cirrus cloud retrieval in the form of error covariance matrices: reflectance uncertainty due to water vapor, the reflectance uncertainty due to ice crystal scattering assumptions, and the instrument measurement uncertainty. Results show that the uncertainty due to habit selection is the largest, while that due to water vapor is at most 0.6% relative to channel reflectance. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water.

Cirrus cloud↗
NASA NTRS1997

Computational Considerations for the Simulation of Discontinuous Flows

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. In a previous work the capabilities and deficiencies of shock-capturing methods for aeroacoustic problems were demonstrated using a high-order essentially nonoscillatory (ENO) numerical method. It was shown that first-order results are obtained when simulating time-dependent flows with discontinuities. The present study reaffirms this conclusion by comparing the ENO results with those obtained using a conventional linear scheme. A sixth-order-accurate compact implicit finite difference scheme is used to investigate various discontinuous flows. The design order of accuracy is achieved in the smooth regions of a steady-state, quasi-one-dimensional Euler test case, as well as in the time-dependent Burgers' equation. However, in the unsteady Euler sound-shock interaction, first-order results are obtained downstream of the shock. A comparison is made between the linear and nonlinear results, noting the advantages of each method. A discontinuous linear model problem is then used to identify the cause of the first-order results. Here, the nature of the solution error is quantified as being predominantly a numerical phase shift, and a post-processing procedure is demonstrated which increases the solution accuracy downstream of the discontinuity to second-order.

Carpenter Mark H.↗
NASA NTRS2006

Three Decades of Precision Orbit Determination Progress, Achievements, Future Challenges and its Vital Contribution to Oceanography and Climate Research

Although satellite altimetry has been around for thirty years, the last fifteen beginning with the launch of TOPEX/Poseidon (TP) have yielded an abundance of significant results including: monitoring of ENS0 events, detection of internal tides, determination of accurate global tides, unambiguous delineation of Rossby waves and their propagation characteristics, accurate determination of geostrophic currents, and a multi-decadal time series of mean sea level trend and dynamic ocean topography variability. While the high level of accuracy being achieved is a result of both instrument maturity and the quality of models and correction algorithms applied to the data, improving the quality of the Climate Data Records produced from altimetry is highly dependent on concurrent progress being made in fields such as orbit determination. The precision orbits form the reference frame from which the radar altimeter observations are made. Therefore, the accuracy of the altimetric mapping is limited to a great extent by the accuracy to which a satellite orbit can be computed. The TP mission represents the first time that the radial component of an altimeter orbit was routinely computed with an accuracy of 2-cm. Recently it has been demonstrated that it is possible to compute the radial component of Jason orbits with an accuracy of better than 1-cm. Additionally, still further improvements in TP orbits are being achieved with new techniques and algorithms largely developed from combined Jason and TP data analysis. While these recent POD achievements are impressive, the new accuracies are now revealing subtle systematic orbit error that manifest as both intra and inter annual ocean topography errors. Additionally the construction of inter-decadal time series of climate data records requires the removal of systematic differences across multiple missions. Current and future efforts must focus on the understanding and reduction of these errors in order to generate a complete and consistent time series of improved orbits across multiple missions and decades required for the most stringent climate-related research. This presentation discusses the POD progress and achievements made over nearly three decades, and presents the future challenges, goals and their impact on altimetric derived ocean sciences.

Luthcke, Scott↗
NASA NTRS2022

Precision Measurement of the Nickel Spectrum in Cosmic Rays from 8.8 GeV/n to 240 GeV/n with CALET on the International Space Station

The relative abundance of nickel with respect to iron is by far larger than all other trans-iron elements, therefore it provides a favourable opportunity for a low background measurement of its spectrum. Since nickel, as well as iron, is one of the most stable nuclei, the nickel energy spectrum and its relative abundance with respect to iron provide important information to estimate the abundances at the source and to model the propagation of heavy nuclei. However, only a few direct measurements of cosmic-ray nickel at energy larger than∼3 GeV/n are available at present in the literature and they are affected by strong limitations in both energy reach and statistics. In this paper we present a measurement of the differential energy spectrum of nickel from 8.8 GeV/nto 240 GeV/n, carried out with unprecedented precision by the Calorimetric Electron1 Telescope(CALET) in operation on the International Space Station since 2015. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number Z = 40). The particle’s energy is measured by a homogeneous calorimeter (1.2 proton interaction lengths, 27 radiation lengths) preceded by a thin imaging section (3 radiation lengths) providing tracking and energy sampling. This paper follows our previous measurement of the iron spectrum [1] and it extends our investigation on the energy dependence of the spectral index of heavy elements. It reports the analysis of nickel data collected from November 2015 to May 2021 and a detailed assessment of the systematic uncertainties. In the region from 20 GeV/nto 240 GeV/n our present data are compatible within the errors with a single power law with spectral index−2.51±0.07.

O Adriani↗
NASA NTRS2010

Lunar Reconnaissance Orbiter (LRO) Thruster Control Mode Design and Flight Experience

National Aeronautics and Space Administration s (NASA) Goddard Space Flight Center (GSFC) in Greenbelt, MD, designed, built, tested, and launched the Lunar Reconnaissance Orbiter (LRO) from Cape Canaveral Air Force Station on June 18, 2009. The LRO spacecraft is the first operational spacecraft designed to support NASA s return to the Moon, as part of the Vision for Space Exploration. LRO was launched aboard an Atlas V 401 launch vehicle into a direct insertion trajectory to the Moon. Twenty-four hours after separation the propulsion system was used to perform a mid-course correction maneuver. Four days after the mid-course correction a series of propulsion maneuvers were executed to insert LRO into its commissioning orbit. The commission period lasted eighty days and this followed by a second set of thruster maneuvers that inserted LRO into its mission orbit. To date, the spacecraft has been gathering invaluable data in support of human s future return to the moon. The LRO Attitude Control Systems (ACS) contains two thruster based control modes: Delta-H and Delta-V. The design of the two controllers are similar in that they are both used for 3-axis control of the spacecraft with the Delta-H controller used for momentum management and the Delta-V controller used for orbit adjust and maintenance maneuvers. In addition to the nominal purpose of the thruster modes, the Delta-H controller also has the added capability of performing a large angle slew maneuver. A suite of ACS components are used by the thruster based control modes, for both initialization and control. For initialization purposes, a star tracker or the Kalman Filter solution is used for providing attitude knowledge and upon entrance into the thruster based control modes attitude knowledge is provided via rate propagation using a inertial reference unit (IRU). Rate information for the controller is also supplied by the IRU. Three-axis control of the spacecraft in the thruster modes is provided by eight 5-lbf class attitude control thrusters configured in two sets of four thrusters for redundancy purposes. Four additional 20-lbf class thrusters configured in two sets of two thrusters are used for Lunar Orbit Insertion maneuvers. The propulsion system is one the few systems on-board the LRO spacecraft that has built in redundancy. The Delta-H controller consists of a Proportional-Derivative (PD) controller with a structural filter on the thrusters and a Proportional controller on the reaction wheels. The PD control that employs the thrusters is used for attitude and rate control. The Proportional controller on the reaction wheels is used for commanding the wheels to a new momentum state. The ground commands used for the Delta-H controller are the system momentum vector, reaction wheel momentum, maximum expected command time, and which set of attitude control thrusters to use. The ability to command both the system momentum vector and reaction wheel momentum in the Delta-H controller provides both a capability and an additional source of operator error. Large angle slews via the Delta-H controller is achievable via this commands because these commands are used for the exit mode criteria. Setting these commands to non-consistent values prevents the mode from exiting nominally.

Hsu, Oscar C.↗
NASA NTRSDate not supplied

Toward Time Synchronization in Delay Tolerant Network based Solar System Internetworking

The expanding presence in space will place an increased dependency on networked communications– a scalable communications infrastructure; that is, the Solar System Internet (SSI). Upcoming developments towards a SSI include NASA’s upcoming LunaNet, or lunar Internet, which provides multi-hop multi-path communications using Delay Tolerant Networking (DTN). DTN has been an active area of research and development, particularly in routing, security, and optimization. DTNs are marked by mobility, disconnection, and a wide variance of latencies (propagation and processing delays). In this paper, we outline progress towards a theory of time synchronization across such a network. An underlying assumption of DTN is that the network is time synchronized already, rather than synchronization being provided as a service. While this is necessary for schedule-based routing, which is necessarily prevalent in DTNs, it is so deeply ingrained as to be built into the primary unit of data in DTNs– the bundle. Indeed, a bundle’s creation timestamp and its time to live (called the lifetime) are based on time, and there are special recommendations for systems that lack accurate clocks. The assumption of time synchronization makes sense when limiting considerations to smaller-scale and more traditional space communication. However, just as end-to-end connectivity cannot be guaranteed in DTNs, neither can access to a reference or authoritative clock. In this more general case, it might be necessary to synchronize over time-varying meshes, and perhaps even to consider relativistic effects. Moreover, by imposing synchronization restrictions in order to sustain a network, the effectiveness of the network to achieve scalability will be necessarily muted. To work towards a time synchronization theory for DTNs, we build upon past successes in modeling DTNs using time-varying graphs and sheaves. This includes error and limitation estimation, which allows one to define domains over which schedule-based routing is possible, up to some threshold sensitivity. Despite the theoretical nature of these results, the approaches taken are also algorithmic, and hence lend themselves to practical implementations. The paper concludes with comparisons of the various methods along with suggestions for future work.

Delay Tolerant Networking↗
NASA NTRS2004

Computation of Gust-Cascade Interaction Using the CE/SE Method

The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.

Wang, X.-Y.↗
NASA NTRS2000

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗
NASA NTRS2007

The Effect of Acoustic Disturbances on the Operation of the Space Shuttle Main Engine Fuel Flowmeter

The Space Shuttle Main Engine (SSME) uses a turbine fuel flowmeter (FFM) in its Low Pressure Fuel Duct (LPFD) to measure liquid hydrogen flowrates during engine operation. The flowmeter is required to provide accurate and robust measurements of flow rates ranging from 10000 to 18000 GPM in an environment contaminated by duct vibration and duct internal acoustic disturbances. Errors exceeding 0.5% can have a significant impact on engine operation and mission completion. The accuracy of each sensor is monitored during hot-fire engine tests on the ground. Flow meters which do not meet requirements are not flown. Among other parameters, the device is screened for a specific behavior in which a small shift in the flow rate reading is registered during a period in which the actual fuel flow as measured by a facility meter does not change. Such behavior has been observed over the years for specific builds of the FFM and must be avoided or limited in magnitude in flight. Various analyses of the recorded data have been made prior to this report in an effort to understand the cause of the phenomenon; however, no conclusive cause for the shift in the instrument behavior has been found. The present report proposes an explanation of the phenomenon based on interactions between acoustic pressure disturbances in the duct and the wakes produced by the FFM flow straightener. Physical insight into the effects of acoustic plane wave disturbances was obtained using a simple analytical model. Based on that model, a series of three-dimensional unsteady viscous flow computational fluid dynamics (CFD) simulations were performed using the MSFC PHANTOM turbomachinery code. The code was customized to allow the FFM rotor speed to change at every time step according to the instantaneous fluid forces on the rotor, that, in turn, are affected by acoustic plane pressure waves propagating through the device. The results of the simulations show the variation in the rotation rate of the flowmeter due to the interaction of the flow straightener wakes and the upstream propagating acoustic waves. A detailed analysis of the acoustic disturbance effects is presented along with an assessment of the impact on measurement accuracy.

Marcu, Bogdan↗
NASA NTRS2011

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described.

Kurum, Mehmet↗
NASA NTRS2016

Impact of Scattering Model on Disdrometer Derived Attenuation Scaling

NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 GHz attenuation from the disdrometer and the 20 GHz timeseries as well as to directly measure the 40 GHz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data.In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.

statistics↗
NASA NTRS2016

Impact of Scattering Model on Disdrometer Derived Attenuation Scaling

NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP#5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 gigahertz attenuation from the disdrometer and the 20 gigahertz time-series as well as to directly measure the 40 gigahertz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data. In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer-derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.

Radio Transmission↗
NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'prop-burn-prop' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons are made with three methods using modified linearized covariance propagation. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. The second method does not model the maneuver but includes a process noise matrix to account for the uncertainty in its magnitude. The third method, which is essentially a hybrid of the first two, includes the nominal portion of the maneuver via the state transition matrix and uses a process noise matrix to account for the magnitude uncertainty. The first method is unable to produce the final orbit covariance except in the case of zero maneuver uncertainty. The second method yields good accuracy for the final covariance matrix but fails to model the final orbital state accurately. Agreement between the simulated covariance data produced by this method and the Monte Carlo truth data fell within 0.5-2.5 percent over a range of maneuver sizes that span two orders of magnitude (0.1-20 m/s). The third method, which yields a combination of good accuracy in the computation of the final covariance matrix and correct accounting for the presence of the maneuver in the nominal orbit, is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, PC. However, applications for the two other methods exist and are briefly discussed. Although the process model ("prop-burn-prop") that was studied is very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
NASA NTRS2021

Validation of Small Satellite Dynamics Simulation Modules using ASTERIA Flight Data

ASTERIA (Arcsecond Space Telescope Enabling Research in Astrophysics) was a CubeSat space telescope that operated in low-Earth orbit, having been deployed from the International Space Station in 2017. The spacecraft has achieved sub-arcsecond pointing stability and millikelvin thermal stability over 20-minute observations. A key enabling feature of the ASTERIA mission is the low level of pointing error that has been demonstrated under ASTERIA’s fine pointing control mode. Prior to launch, the ASTERIA mission performed analysis and simulations to estimate the in-flight pointing performance. This analysis used reaction wheel models provided by JPL’s Small Satellite Dynamics Testbed (SSDT), which had performed Kistler table testing to characterize the jitter caused by the reaction wheels of the Blue Canyon Technologies XACT attitude control unit within ASTERIA. These models were subsequently incorporated into the SSDT simulation. The main source of jitter, the high-frequency attitude disturbance over the camera exposure time, on ASTERIA is the set of rotating reaction wheels that are used for maintaining fine pointing during observations. The reaction wheels impart disturbance forces and torques continuously that cause unwanted motion during imaging. Increased jitter consequently leads to blurring of the image. Therefore, to benefit all SmallSat missions devoted to photometric or spectroscopic astrophysics applications, an on-orbit flight data acquisition and jitter testing campaign was performed in an attempt to help validate the SSDT’s simulation models. This paper describes the process used to validate several simulation models and attempt to characterize the ASTERIA jitter by analyzing the size of the resulting spot size for the target stars as a function of wheel speed, having commanded four wheel speeds for each of three stars of known brightness. After summarizing the on-orbit observation and jitter level measurement process, this paper compares the results obtained from the flight operational environment to the corresponding set of simulation jitter levels. Though the flight data obtained during this experiment was not sufficient to validate the simulation’s jitter model, other models have been validated with flight data, including the magnetic field, orbit propagation, sun position, and long-duration orbit decay models, which may be used by all other projects that use the SSDT’s simulation to increase system capability knowledge.

Mohan, Swati↗
NASA NTRS2004

Reannealed Fiber Bragg Gratings Demonstrated High Repeatability in Temperature Measurements

Fiber Bragg gratings (FBGs) are formed by periodic variations of the refractive index of an optical fiber. These periodic variations allow an FBG to act as an embedded optical filter, passing the majority of light propagating through a fiber while reflecting back a narrow band of the incident light. The peak reflected wavelength of the FBG is known as the Bragg wavelength. Since the period and width of the refractive index variation in the fiber determines the wavelengths that are transmitted and reflected by the grating, any force acting on the fiber that alters the physical structure of the grating will change the wavelengths that are transmitted and reflected by it. Both thermal and mechanical forces acting on the grating will alter its physical characteristics, allowing the FBG sensor to detect both the temperature variations and the physical stresses and strains placed upon it. This ability to sense multiple physical forces makes the FBG a versatile sensor. To assess the feasibility of using Bragg gratings as temperature sensors for propulsion applications, researchers at the NASA Glenn Research Center evaluated the performance of Bragg gratings at elevated temperatures for up to 300 C. For these purposes, commercially available polyimide-coated high-temperature gratings were used that were annealed by the manufacturer to 300 C. To assure the most thermally stable gratings at the operating temperatures, we reannealed the gratings to 400 C at a very slow rate for 12 to 24 hr until their reflected optical powers were stabilized. The reannealed gratings were then subjected to periodic thermal cycling from room temperature to 300 C, and their peak reflected wavelengths were monitored. The setup shown is used for reannealing and thermal cycling the FBGs. Signals from the photodetectors and the spectrum analyzer were fed into a computer equipped with LabVIEW software. The software synchronously monitored the oven/furnace temperature and the optical spectrum analyzer as well as processed the data. Experimental results presented in the following graph show typical wavelength versus temperature dependence of a reannealed FBG through six thermal cycles (80 hr). The average standard deviation of the temperature-to-wavelength relationship ranged from 1.86 to 2.92 C over the six thermal cycles each grating was subjected to. This is an error of less than 1.0 percent of full scale throughout the entire evaluation temperature range from ambient to 300 C.

Adamovsky, Grigory↗