Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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,441 records · Page 80

Centaur zero gravity coast and engine restart demonstration on the Titan/Centaur (TC-2) extended mission

The Centaur propellant management and thermal control techniques required for zero gravity coasting were successfully demonstrated during an extended mission following spacecraft separation on the TC-2 flight. As part of the demonstration, two successful engine starts were accomplished. The first engine start followed a one-hour zero gravity coast, and the second engine start followed a three-hour zero gravity coast. All of the Centaur systems performed satisfactorily, the design parameters for zero gravity coasting were verified, and no significant problems were encountered. The flight results showed that the propellant location and behavior, propellant heating, and tank pressure rise rates observed during the zero-gravity coasts were less severe than expected. Consequently, the majority of the propellants remained at the tank bottom, the propellant collection times were very short, and more than 7 hours of coast could have been achieved before a tank venting was required. The tank pressurization prior to the engine starts provided boost pump Net Positive Suction Head values well in excess of the values required. The LO2 tank pressurization was accomplished by a new bubbler method that greatly reduced the helium usage.

Lacovic, R. F.↗

Far ultraviolet observations of the expanding shell in Eridanus

After correcting for hot stars in the instrument field of view, 1350-1750 A wavelength band observations of a thin strip of the sky near the south galactic pole have yielded count rates which vary smoothly from a minimum of about 55/cs near l = 30 deg, b = -60 deg to a maximum of about 740/cs in the direction of l = 199 deg, b = -49 deg. These values correspond to specific intensities of 600-8000 photons/sq cm per sec per sr per A, assuming a flat incident spectrum in this wavelength range. The magnitude and spatial variation of the observed signal are incompatible with any known source of terrestrial far-UV airglow, zodiacal light, or with the expected emission of high latitude reflection nebulosities. Attention is given to alternative sources that, together with the reflection nebulae, may explain the high intensities in this region.

Paresce, F.↗

11-year cycle solar modulation of cosmic ray intensity inferred from C-14 content variation in dated tree rings

A liquid scintillation-photomultiplier tube counter system was used to measure the Delta-C-14 values of 60 tree rings, dating from 1866 to 1925, that were taken from a white spruce grown in Canada at 68 deg N, 130 deg W. A 10-percent variation is found which is anticorrelated with sunspot numbers, although the amplitude of the variation is 2-3 times higher than expected in trees grown at lower latitudes. A large dip in the data at about 1875 suggests an anomalously large modulation of cosmic ray intensity during the 1867-1878 AD solar cycle, which was the most active of the 19th century.

Fan, C. Y.↗

Timing the Geminga Pulsar with Gamma-Ray Observations

We present the COS-B/EGRET 1997 ephemeris for the rotation of the Geminga pulsar. This ephemeris is derived from high-energy gamma-ray observations that span 24 yr. The recently obtained accurate position and proper motion are assumed. A cubic ephemeris predicts the rotational phase of Geminga with errors smaller than 50 milliperiods for all existing high-energy gamma-ray observations that span a 24.2 yr timing baseline. The braking index obtained is 17 +/- 1. Further observation is required to ascertain whether this high value truly reflects the rotational energy loss mechanism, or whether it is a manifestation of timing noise. The ephemeris parameters are sufficiently constrained so that timing noise will be the limitation on forward extrapolation. If Geminga continues to rotate without a glitch, as it has for at least 23 yr, we expect this ephemeris to continue to describe the phase, with an error of less than 100 milliperiods, until 2008. Statistically significant timing residuals are detected in the EGRET data that depart from the cubic ephemeris at a level of 30 milliperiods. Although this could simply be an additional manifestation of timing noise, the EGRET timing residuals appear to have a sinusoidal modulation that is consistent with a planet of mass 1.7/sin i solar mass, orbiting Geminga at a radius of 3.3 AU.

Mattox, J. R.↗

Real-Time Simulation for Verification and Validation of Diagnostic and Prognostic Algorithms

To verify that a health management system (HMS) performs as expected, a virtual system simulation capability, including interaction with the associated platform or vehicle, very likely will need to be developed. The rationale for developing this capability is discussed and includes the limited capability to seed faults into the actual target system due to the risk of potential damage to high value hardware. The capability envisioned would accurately reproduce the propagation of a fault or failure as observed by sensors located at strategic locations on and around the target system and would also accurately reproduce the control system and vehicle response. In this way, HMS operation can be exercised over a broad range of conditions to verify that it meets requirements for accurate, timely response to actual faults with adequate margin against false and missed detections. An overview is also presented of a real-time rocket propulsion health management system laboratory which is available for future rocket engine programs. The health management elements and approaches of this lab are directly applicable for future space systems. In this paper the various components are discussed and the general fault detection, diagnosis, isolation and the response (FDIR) concept is presented. Additionally, the complexities of V&V (Verification and Validation) for advanced algorithms and the simulation capabilities required to meet the changing state-of-the-art in HMS are discussed.

Aguilar, Robet↗

Sunspot Activity Near Cycle Minimum and What it Might Suggest for Cycle 24, the Next Sunspot Cycle

In late 2008, 12-month moving averages of sunspot number, number of spotless days, number of groups, area of sunspots, and area per group were reflective of sunspot cycle minimum conditions for cycle 24, these values being of or near record value. The first spotless day occurred in January 2004 and the first new-cycle, high-latitude spot was reported in January 2008, although old-cycle, low-latitude spots have continued to be seen through April 2009, yielding an overlap of old and new cycle spots of at least 16 mo. New-cycle spots first became dominant over old-cycle spots in September 2008. The minimum value of the weighted mean latitude of sunspots occurred in May 2007, measuring 6.6 deg, and the minimum value of the highest-latitude spot followed in June 2007, measuring 11.7 deg. A cycle length of at least 150 mo is inferred for cycle 23, making it the longest cycle of the modern era. Based on both the maximum-minimum and amplitude-period relationships, cycle 24 is expected to be only of average to below-average size, peaking probably in late 2012 to early 2013, unless it proves to be a statistical outlier.

Wilson, Robert M.↗

WMO Assessment 2022: Chapter 3 Summary

Efforts for the 2022 Scientific Assessment on Ozone Depletion published by the World Meteorological Organization are well underway. Chapter 3, titled “Update on Global Ozone: Past, Present, and Future”, presents our current understanding of global ozone outside of the polar regions and this talk summarizes some of the highlights and changes from the previous Assessment in 2018. Since the last Assessment, the longer observational records show a small, significant increase in near-global total column ozone. A small significant increase in total column ozone is also seen in the mid-latitudes of both hemispheres, but not yet the tropics. Different processes operating at different altitudes complicate the interpretation of the overall total column trend. However, a significant increase in upper stratospheric ozone noted in the previous Assessment continues, driven by declines in ozone depleting substances and increases in greenhouse gases. Model simulations support our understanding of these trends. Over the coming century, we expect an increase in global stratospheric ozone as the concentrations of ozone depleting substances decline. The future evolution for different latitudes and vertical levels depends on the future concentrations of greenhouse gases and precursors of tropospheric ozone., These other influences may lead to total ozone column levels that remain below 1980 values even after ozone depleting substance concentrations have recovered to pre-1980 levels.

Birgit Hassler↗

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low earth orbit for scientific, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with 8 Teledyne dosimeters over varying shielding effectiveness. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 –AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat 3 unit (U) models. Using NOVICE model radiation analysis coding, the shielding effectiveness's, based on a generic CubeSat 3U structure with 4 electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cmto4.44 cm (0.550 g/cm2–12.0 g/cm2) thick aluminum. For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID)increases by nearly a magnitude for thin-walled aluminum 0.550g/cm2-0.686 g/cm2(0.204 cm –0.254 cm) typical CubeSat structures. The shielding effectiveness by NOVICE Sigma estimates, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from Earth's magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast, to a thin walled generic 3U CubeSat, Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm2from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the spacecraft’s mission value versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗

Development of a Data Fusion Methodology for Lineload Aerodynamic Databases for a Launch Vehicle during Liftoff and Transition

The need for databases for the distributed loading on launch vehicles during the early portion of flight necessitates the use of expensive computational flows in regimes where wake effects dominate. While also being expensive, this is a regime that computational tools tend to historically have problems simulating accurately. To help tackle this problem, a method of data fusion to combine computational results to wind tunnel derived force and moment data is developed. Using this method, significant reduction in computational costs and increases in confidence of the final product is possible and has been used to generate several databases for the Space Launch System (SLS) at NASA. While the full details of database generation are not part of this work, the crucial method at its core is developed here. Two SLS geometries are used throughout the work to demonstrate the techniques. These are two of the larger geometries and represent both planned crewed missions to the Moon as well as potential cargo missions to deep space. The method uses principal component analysis (PCA) to generate a reduced ordered model (ROM) to help fill in the full parameter space. Other similar techniques are explored, but were not found to have a significant result on the predictions of the ROM. Because the full number of components are kept to generate the model, this lack of difference is expected. This method is then extended to ensure that predicted surfaces match trusted force and moment data derived from wind tunnel testing. This extension is done by setting up a constrained optimization problem in order to minimize the deviation from the surface resolved computational data while still integrating to the desired values. When generating the constrained optimization problem, a weighting factor to balance these competing needs is introduced. The work compares previously introduced weighting terms from similar work to the proposed terms and shows that the previously used terms do not have as desirable behavior in this flow regime. This method is then expanded by developing a technique to incorporate uncertainty quantification into the developed data fusion methodology. This expansion takes a two pronged approach. One examines transferring the uncertainties in the force and moment database and characterizes how those adjustments change the predicted lineloads. The second looks at model form error and looks how rebuilding the model using slightly different data changes the predictions. These two terms are then combined in order to create an uncertainty model that takes both effects into account. The limitations of the proposed methods is then discussed as well as possible techniques to address these shortcomings.

Launch Vehicles↗

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗

Simplified model of statistically stationary spacecraft rotation and associated induced gravity environments

A stochastic model of spacecraft motion was developed based on the assumption that the net torque vector due to crew activity and rocket thruster firings is a statistically stationary Gaussian vector process. The process had zero ensemble mean value, and the components of the torque vector were mutually stochastically independent. The linearized rigid-body equations of motion were used to derive the autospectral density functions of the components of the spacecraft rotation vector. The cross-spectral density functions of the components of the rotation vector vanish for all frequencies so that the components of rotation were mutually stochastically independent. The autospectral and cross-spectral density functions of the induced gravity environment imparted to scientific apparatus rigidly attached to the spacecraft were calculated from the rotation rate spectral density functions via linearized inertial frame to body-fixed principal axis frame transformation formulae. The induced gravity process was a Gaussian one with zero mean value. Transformation formulae were used to rotate the principal axis body-fixed frame to which the rotation rate and induced gravity vector were referred to a body-fixed frame in which the components of the induced gravity vector were stochastically independent. Rice's theory of exceedances was used to calculate expected exceedance rates of the components of the rotation and induced gravity vector processes.

Fichtl, G. H.↗

Applications of parameter estimation in the study of spinning airplanes

Spinning airplanes offer challenges to estimating dynamic parameters because of the nonlinear nature of the dynamics. In this paper, parameter estimation techniques are applied to spin flight test data for estimating the error in measuring post-stall angles of attack, deriving Euler angles from angular velocity data, and estimating nonlinear aerodynamic characteristics. The value of the scale factor for post-stall angles of attack agrees closely with that obtained from special wind-tunnel tests. The independently derived Euler angles are seen to be valid in spite of steep pitch angles. Estimates of flight derived nonlinear aerodynamic parameters are evaluated in terms of the expected fit error.

W Taylor, L., Jr.↗

LN2 spray droplet size measurement via ensemble diffraction technique

The size of subcooled liquified nitrogen droplets are measured with a 5 mW He-Ne laser as a function of pressure difference (delta P) across flat spray and full cone pressure atomizing nozzles. For delta P's of 3 to 30 psid, the spray sauter mean diameter (SMD) ranged between 250 to 50 microns. The pressure range tested is representative of those expected during cryogenic fluid transfer operations in space. The droplet sizes from the flat spray nozzles were greater than those from the full cone nozzle. A power function of the form, SMD varies as delta P(exp a), describes the spray SMD as a function of the delta P very well. The values of a were -0.36 for the flat spray and -0.87 for the full cone. The reduced dependence of the flat spray SMD on the delta P was probably because of: (1) the absence of a swirler that generates turbulence within the nozzle to enhance atomization, and (2) a possible increase in shearing stress resulting from the delayed atomization due to the absence of turbulence. The nitrogen quality, up to 1.5 percent is based on isenthalpic expansion, did not have a distinct and measurable effect on the spray SMD. Both bimodal and monomodal droplet size population distributions were measured. In the bimodal distribution, the frequency of the first mode was much greater than the frequency of the second mode. Also, the frequency of the second mode was low enough such that a monomodal approximation probably would give reasonable results.

Saiyed, N. H.↗

LN2 spray droplet size measurement via ensemble diffraction technique

The size of subcooled liquified nitrogen droplets are measured with a 5 mW He-Ne laser as a function of pressure difference (delta P) across flat spray and full cone pressure atomizing nozzles. For delta P's of 3 to 30 psid, the spray sauter mean diameter (SMD) ranged between 250 to 50 microns. The pressure range tested is representative of those expected during cryogenic fluid transfer operations in space. The droplet sizes from the flat spray nozzles were greater than those from the full cone nozzle. A power function of the form, SMD varies as delta P(exp a), describes the spray SMD as a function of the delta P very well. The values of a were -0.36 for the flat spray and -0.87 for the full cone. The reduced dependence of the flat spray SMD on the delta P was probably because of: (1) the absence of a swirler that generates a turbulence within the nozzle to enhance atomization, and (2) a possible increase in shearing stress resulting from the delayed atomization due to the absence of turbulence. The nitrogen quality, up to 1.5 percent is based on isenthalpic expansion, did not have a distinct and measurable effect on the spray SMD. Both bimodal and monomodal droplet size population distributions were measured. In the bimodal distribution, the frequency of the first mode was much greater than the frequency of the second mode. Also, the frequency of the second mode was low enough such that a monomodal approximation probably would give reasonable results.

Saiyed, N. H.↗

The atmospheric load response of the ocean determined using Geosat altimeter data

Sea-surface heights derived from the Geosat altimeter data for the November 1986 - November 1987 period are compared with the Fleet Numerical Oceanographic Center and the NMC sea-level atmospheric pressure fields, to verify the open ocean inverted barometer (IB) response. It is shown that fits along the satellite track consistently yield results that are only 60-70 percent of the IB results. Fits at fixed geographic locations show a clear dependence of the absolute value of the regression coefficient with latitude: a steady decrease from 70 deg to about 20 deg, and then a sudden increase again nearer to the equator, with similar results obtained for the Northern and the Southern Hemispheres. The sign of the departure from IB is found to be consistent with the expected effects of errors in the pressure fields and with wind-induced changes in sea level that are correlated with pressure fluctuations.

Vandam, Tonie M.↗

Solar measurements from the Airglow-Solar Spectrometer Instrument (ASSI) on the San Marco 5 satellite

The analysis of the solar spectral irradiance from the Airglow-Solar Spectrometer Instrument (ASSI) on the San Marco 5 satellite is the focus for this research grant. A pre-print copy of the paper describing the calibrations of and results from the San Marco ASSI is attached to this report. The calibration of the ASSI included (1) transfer of photometric calibration from a rocket experiment and the Solar Mesosphere Explorer (SME), (2) use of the on-board radioactive calibration sources, (3) validation of the ASSI sensitivity over its field of view, and (4) determining the degradation of the spectrometers. We have determined that the absolute values for the solar irradiance needs adjustment in the current proxy models of the solar UV irradiance, and the amount of solar variability from the proxy models are in reasonable agreement with the ASSI measurements. This research grant also has supported the development of a new solar EUV irradiance proxy model. We expected that the magnetic flux is responsible for most of the heating, via Alfen waves, in the chromosphere, transition region, and corona. From examining time series of solar irradiance data and magnetic fields at different levels, we did indeed find that the chromospheric emissions correlate best with the large magnetic field levels.

Woods, Thomas N.↗

Asteroid approach covariance analysis for the Clementine mission

The Clementine mission is designed to test Strategic Defense Initiative Organization (SDIO) technology, the Brilliant Pebbles and Brilliant Eyes sensors, by mapping the moon surface and flying by the asteroid Geographos. The capability of two of the instruments available on board the spacecraft, the lidar (laser radar) and the UV/Visible camera is used in the covariance analysis to obtain the spacecraft delivery uncertainties at the asteroid. These uncertainties are due primarily to asteroid ephemeris uncertainties. On board optical navigation reduces the uncertainty in the knowledge of the spacecraft position in the direction perpendicular to the incoming asymptote to a one-sigma value of under 1 km, at the closest approach distance of 100 km. The uncertainty in the knowledge of the encounter time is about 0.1 seconds for a flyby velocity of 10.85 km/s. The magnitude of these uncertainties is due largely to Center Finding Errors (CFE). These systematic errors represent the accuracy expected in locating the center of the asteroid in the optical navigation images, in the absence of a topographic model for the asteroid. The direction of the incoming asymptote cannot be estimated accurately until minutes before the asteroid flyby, and correcting for it would require autonomous navigation. Orbit determination errors dominate over maneuver execution errors, and the final delivery accuracy attained is basically the orbit determination uncertainty before the final maneuver.

Ionasescu, Rodica↗

Performance of high-recovery recycling reverse osmosis with wash water

Inclusion of a recycling loop for partially-desalted water from second-stage reverse-osmosis permeate has been shown useful for achieving high-recovery at moderate applied pressures. This approach has now been applied to simulated wash waters, to obtain data on retention by the membranes of solutes in a mixture comparable to anticipated spacecraft hygiene wastewaters, and to generate an estimate of the maximum concentration that can be expected without causing membrane fouling. A first experiment set provides selectivity information from a single membrane and an Igepon detergent, as a function of final concentration. A reject concentration of 3.1% Total Organic Carbon has been reached, at a pressure of 1.4 Mega Pascals, without membrane fouling. Further experiments have generated selectivity values for the recycle configuration from two washwater simulations, as a function of applied pump pressure. Reverse osmosis removal has also been tested for washwater containing detergent formulated for plant growth compatibility (containing nitrogen, phosphorous and potassium functional groups.)

Herrmann, Cal C.↗