Investigation of nuclear light bulb startup and engine dynamics
Analysis of nuclear light bulb startup and engine dynamics including transient responses to perturbations at full power operating level
SEARCH · Engineering Papers
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.
Analysis of nuclear light bulb startup and engine dynamics including transient responses to perturbations at full power operating level
The results of a study which was designed to monitor, characterize, and evaluate the chemical composition of precipitation (rain) which fell at the Kennedy Space Center, Florida (KSC) during the period July 1977 to March 1979 are reported. Results which were obtained from a soil sampling and associated chemical analysis are discussed. The purpose of these studies was to determine the environmental perturbations which might be caused by NASA space activities.
A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.
The differences in the phase of the velocity oscillations between a pair of chromospheric Ca II lines was measured using the Vacuum Tower Telescope at the Sacramento Peak Observatory. The observed phase differences indicate that the acoustic modes are trapped or envanescent, rather than propagating, in the chromosphere. Systematic distinctions are found in the phase delays between quiet network and cell interior regions for both intensity and velocity oscillations in photospheric and chromospheric lines. The theory of linear perturbations in an isothermal atmosphere is invoked to interpret these differences. From this analysis it is found that one or more of the following explanations is possible: (1) the radiative damping is more effective in the network than in the cell interior; (2) the network features exclude oscillations of large horizontal wavenumber; or (3) the scale height of the chromosphere is larger in the network than in the cell interior.
For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.
The physical roles of 'buoyant' and 'dynamic' pressure components, and the distinction between buoyant and hydrostatic pressure perturbations, are aspects of the pressure perturbation field in strongly sheared convective storms studied by means of two- and three-dimensional anelastic numerical modeling experiments with common environmental profiles. The pressure analysis clarifies the differences between two- and three-dimensional storms. In the main updraft, strong midlevel thermal buoyancy is partly opposed by a downward-perturbed vertical pressure gradient force. This, however, occurs to a much greater extent in two dimensions than in three, contributing to smaller net upward accelerations. While the buoyant and hydrostatic fields are intimately related to the total buoyancy distribution, the buoyant pressure perturbation is smoother and of lower amplitude than its hydrostatic counterpart. For the model experiments, this distinction is far greater in three dimensions than in two, in association with the smaller scale of the active convection in three dimensions.
During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.
A filament in active region AR 2717 was observed in two lines formed at different temperatures (H-alpha at 10 exp 4 K and C IV at 10 exp 5 K) with the multichannel-double-pass (MSDP) spectrograph and the ultraviolet spectrometer and polarimeter (UVSP). The partial disparition brusque (DB) of the filament as observed in H-alpha was due to the heating of a filament section in the vicinity of a bright point. A propagating disturbance followed this event. A detailed analysis of C IV rasters shows that this disturbance was not a passive perturbation, but was itself triggering active phenomena at various locations along its path, resulting in energy releases. It is suggested that this propagation of brightness was due to fast successive reconnections between fine looplike structures of the filament.
To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.
A Precision Departure Release Capability concept is being evaluated by both the National Aeronautics and Space Administration and the Federal Aviation Administration as part of a larger goal of improving throughput, efficiency and capacity in integrated departure, arrival and surface operations. The concept is believed to have the potential of increasing flight efficiency and throughput by avoiding missing assigned slots and minimizing speed increase or path stretch to recover the slot. The main thrust of the paper is determining the impact of early and late departures from the departure runway when an aircraft has a slot assigned either at a meter fix or at the arrival airport. Results reported in the paper are for two scenarios. The first scenario considers flights out of Dallas/Fort Worth destined for Hartsfield-Jackson International Airport in Atlanta flying through the Meridian meter-fix in the Memphis Center with miles-in-trail constraints. The second scenario considers flights destined to George Bush Intercontinental/Houston Airport with specified airport arrival rate constraint. Results show that delay reduction can be achieved by allowing reasonable speed changes in scheduling. It was determined that the traffic volume between Dallas/Fort Worth and Atlanta via the Meridian fix is low and the departures times are spread enough that large departure schedule uncertainty can be tolerated. Flights can depart early or late within 90 minutes without accruing much more delay due to miles-in-trail constraint at the Meridian fix. In the Houston scenario, 808 arrivals from 174 airports were considered. Results show that delay experienced by the 16 Dallas/Fort Worth departures is higher if initial schedules of the remaining 792 flights are kept unaltered while they are rescheduled. Analysis shows that the probability of getting the initially assigned slot back after perturbation and rescheduling decreases with increasing standard deviation of the departure delay distributions. Results show that most Houston arrivals can be expected to be on time based on the assumed zero-mean Normal departure delay distributions achievable by Precision Departure Release Capability. In the current system, airport-departure delay, which is the sum of gate-departure delay and taxi-out delay, is observed at the airports. This delay acts as a bias, which can be reduced by Precision Departure Release Capability.
The Miniaturized Laser Heterodyne Radiometer (mini-LHR) is a passive ground instrument that observes the mole fraction of carbon dioxide (CO) and methane (CH)in the atmospheric column by measuring their absorption of sunlight at 1.6 microns. A laser heterodyne radiometer is similar in design to the super heterodyne radio receiver that is well known by ham radio enthusiasts. While not previously a commercial technique, laser heterodyne radiometers have a history of measuring atmospheric trace gases that started in the 1960s. With the commercial availability of inexpensive, low-power, thumbnail-sized lasers developed for the telecommunications industry, it was possible to miniaturize this technique and ultimately commercialize it. The mini-LHR has been under development at NASA GSFC since 2009. During that time, in addition to signicant technical improvements, processing has also evolved and been streamlined. Here we present details of the processing approach for raw mini-LHR data to produce 30- and 60-minute data products of CH and CO column mole fractions. Processing occurs in two general stages: a python-based pre-processing of raw data, followed by ingestion into a Planetary Spectrum Generator (PSG) retrieval algorithm. Raw data processing involves removal of outliers, correcting for changes in air mass throughout the day, averaging scans, and ultimately converting averaged scans into transmittance vs. wavelength. The PSG retrieval simulates a spectra for the time/day/location of the scan with meteorological inputs from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data set and then perturbs the concentrations of CO and CH to obtain a t based on an iterative least-squares curve fitting procedure.
Results are presented which were obtained with the Winge-Coleman model for theoretical predictions of latitudinal dependencies in the solar wind. A first-order expansion is described which allows analysis of first-order latitudinal variations in the coronal boundary conditions and results in a second-order partial differential equation for the perturbation stream function. Latitudinal dependencies are analytically separated out in the form of Legendre polynomials and their derivative, and are reduced to the solution of radial differential equations. This analysis is shown to supply an estimate of how large the coronal variation in latitude must be to produce an 11 km/sec/deg gradient in the radial velocity of the solar wind, assuming steady-state processes.
Atmospheric photochemistry models have been used to predict the sensitivity of the ozone layer to various perturbations. These same models also predict concentrations of chemical species in the present day atmosphere which can be compared to observations. Model results for both present day values and sensitivity to perturbation depend upon input data for reaction rates, photodissociation rates, and boundary conditions. A method of combining the results of a Monte Carlo uncertainty analysis with the existing set of present atmospheric species measurements is developed. The method is used to examine the range of values for the sensitivity of ozone to chlorine perturbations that is possible within the currently accepted ranges for input data. It is found that model runs which predict ozone column losses much greater than 10 percent as a result of present fluorocarbon fluxes produce concentrations and column amounts in the present atmosphere which are inconsistent with the measurements for ClO, HCl, NO, NO2, and HNO3.
This paper presents a new methodology for adjoint sensitivity analysis, previously developed in general terms by Cacuci, into a form directly applicable 10 meteorological problems. This technique is illustrated by examining the sensitivity of a blocking index in a two-layer isentropic model. The index represents a response function for the sensitivity analysis that unlike previous meteorological applications, is an operator and not a functional, and thus extends the scope of adjoint sensitivity to general operator-type responses depending on time and/or space. The sensitivity of the blocking index to perturbations introduced into the model atmosphere, as well as to model parameters, is discussed. The methodology of generalized adjoint sensitivity analysis described in this paper constitutes a prototype for further applications in the atmospheric and/or oceanic sciences.
The navigation and guidance of a high lift-to-drag ratio sample return vehicle during aerocapture at Mars are investigated. Emphasis is placed on integrated systems design, with guidance algorithm synthesis and analysis based on vehicle state and atmospheric density uncertainty estimates provided by the navigation system. The latter utilizes a Kalman filter for state vector estimation, with useful update information obtained through radar altimeter measurements and density altitude measurements based on IMU-measured drag acceleration. A three-phase guidance algorithm, featuring constant bank numeric predictor/corrector atmospheric capture and exit phases and an extended constant altitude cruise phase, is developed to provide controlled capture and depletion of orbital energy, orbital plane control, and exit apoapsis control. Integrated navigation and guidance systems performance are analyzed using a four degree-of-freedom computer simulation. The simulation environment includes an atmospheric density model with spatially correlated perturbations to provide realistic variations over the vehicle trajectory. Navigation filter initial conditions for the analysis are based on planetary approach optical navigation results. Results from a selection of test cases are presented to give insight into systems performance.
The initial boundary value problem describing the evolution of unsteady linearized perturbations of a steady, uniform subsonic flow is analyzed. The eigenmodes and eigenfrequencies of the system are derived and several examples are presented to illustrate the effect of different boundary conditions on the exponential decay rate of the eigenmodes. The resultant implications for the stability and convergence rates of finite difference computations are discussed.
This research was first proposed in May, 1986, to focus on some of the problems encountered in the analysis of postural responses gathered from crewmembers. The ultimate driving force behind this line of research was the desire to treat, predict, or explain 'Space Adaptation Syndrome' (SAS) and hence circumvent any adverse effects of space motion sickness on crewmember performance. The aim of this project was to develop an easily implemented analysis of the transient responses to platform translation that can be elicited with a protocol designed to force sensorimotor reorganization, utilizing statistically reliable criterion measures. This report will present: (1) a summary of the activity that took place in each of the three funded years of the project; (2) discussion of experimental results and their implications for future research; and (3) a list of presentations and publications resulting from this project.
An investigation of the effect of mean local node-crossing time on the evolution of Sun-synchronous orbits was undertaken during Phase-A orbit analysis for the National Oceanic and Atmospheric Administration (NOAA) O,P,Q environmental spacecraft. That analysis added to the growing body of evidence that individual Sun-synchronous missions, at differing node-crossing times, experience nodal drift rates that can differ in both magnitude and direction. A Sun-synchronous orbit is obtained by means of a nodal drift rate approximating the 0.9856-degree-per-day apparent precession of the position of the mean Sun. This drift rate is achieved through the interaction of the orbital semimajor axis and inclination in Earth's geopotential field. Influencing perturbations include atmospheric drag and, most important, the effects of solar gravitation on inclination. The present analysis examines a series of Sun-synchronous orbits with mean local node-crossing times at 1-hour intervals from 6 a.m. to 6 p.m. It considers the fixed geometry of each orbital plane with respect to both the Sun and the diurnal atmospheric bulge, then analyzes the influence of these features upon the evolution of the semimajor axix and inclination and thus upon the rate of the nodal drift in the course of 1 year.