SEARCH · Search NASA
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.
The semiannual report is included. Topics covered include communication, information science, data compression, remote sensing, color mapped images, robust coding scheme for packet video, recursively indexed differential pulse code modulation, image compression technique for use on token ring networks, and joint source/channel coder design.
Automated multidisciplinary design of aircraft and other flight vehicles requires the optimization of complex performance objectives with respect to a number of design parameters and constraints. The effect of these independent design variables on the system performance criteria can be quantified in terms of sensitivity derivatives which must be calculated and propagated by the individual discipline simulation codes. Typical advanced CFD analysis codes do not provide such derivatives as part of a flow solution; these derivatives are very expensive to obtain by divided (finite) differences from perturbed solutions. It is shown that sensitivity derivatives can be obtained accurately and efficiently using the ADIFOR source translator for automatic differentiation. In particular, it is demonstrated that the 3-D, thin-layer Navier-Stokes, multigrid flow solver called TLNS3D is amenable to automatic differentiation in the forward mode even with its implicit iterative solution algorithm and complex turbulence modeling. It is significant that by using computational differentiation, consistent discrete nongeometric sensitivity derivatives have been obtained from an aerodynamic 3-D CFD code in a relatively short time, e.g., O(man-week) not O(man-year).
In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.
A method of determining certain characteristic flexural rigidities and elastic properties of nonlinearly elastic materials is presented. An estimation method utilizing perturbation methods and a least squares fitting technique is used to solve the nonlinear differential equation derived from the moment curvature relation, subject to boundary values representing deflections of the bar at discrete points. Deflection data from numerical simulations of a nonlinearly elastic, prismatic bar are used to demonstrate the estimation method. Numerical experiments relating the accuracy of the identification to the number and accuracy of the boundary values are presented. Conclusions based on the numerical experiments are included.
The generation and use of composite grids for general three-dimensional physical boundary configurations is discussed, and the availability of several codes or procedures is noted. With the composite framework, the physical region is segmented into sub-regions, each bounded by six curved sides, and a grid is generated in each sub-region. These grids may be joined at the interfaces between the sub-regions with various degrees of continuity. This structure allows codes to be constructed to operate on rectangular blocks in computational space, so that existing solution procedures can be readily incorporated in the construction of codes for general configurations. Numerical grid generation is an integral part of the numerical solution of partial differential equations and is one of the pacing items in the development of codes for general configurations. The numerically generated grid frees the computational simulation from restriction to certain boundary shapes and allows general codes to be written in which the boundary shape is specified simply by input. The numerically generated grid allows all computation to be done on a fixed square grid in the computational space, which is always rectangular by construction.
The edge technique is a new and powerful method for measuring small frequency shifts. Range resolved lidar measurements of winds can be made with high accuracy and high vertical resolution using the edge technique to measure the Doppler shift of an atmospheric backscattered signal from a pulsed laser. The edge technique can be used at near-infrared or visible wavelengths using well developed solid state lasers and detectors with various edge filters. In the edge technique, the laser frequency is located on the steep slope of the spectral response function of a high resolution optical filter. Due to the steep slope of the edge, very small frequency shifts cause large changes in measured signal. The frequency of the outgoing laser pulse is determined by measuring its location on the edge of the filter. This is accomplished by sending a small portion of the beam to the edge detection setup where the incoming light is split into two channels - an edge filter and an energy monitor channel. The energy monitor signal is used to normalize the edge filter signal for magnitude. The laser return backscattered from the atmosphere is collected by a telescope and directed through the edge detection setup to determine its frequency (location on the edge) in a similar manner for each range element. The Doppler shift, and thus the wind, is determined from a differential measurement of the frequency of the outgoing laser pulse and the frequency of the laser return backscattered from the atmosphere. We have conducted simulations of the performance of an edge lidar system using an injection seeded pulsed Nd:YAG laser at 1.06 microns. The central fringe of a Fabry-Perot etalon is used as a high resolution edge filter to measure the shift of the aerosol return.
A sensitivity equation method to shape optimization problems is applied. An algorithm is developed and tested on a problem of designing optimal forebody simulators for a 2D, inviscid supersonic flow. The algorithm uses a BFGS/Trust Region optimization scheme with sensitivities computed by numerically approximating the linear partial differential equations that determine the flow sensitivities. Numerical examples are presented to illustrate the method.
The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 2 yr maps produced by the Cosmic Background Explorer Satellite Differential Microwave Radiometer (COBE DMR)). The power spectrum of the real sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power-law form, with P(k) proportional to k(exp n). Within the limited range of spatial scales covered by the COBE DMR, corresponding to spherical harmonic indices 3 less than or = l is less than or approximately = 30, the best-fitting value of the spectral index is n = 1.25(sup +0.39 sub -0.44) with the Harrisson-Zel'dovich value n = 1 approximately 0.5 sigma below the best fit. For 3 less than or = l less than or approximately = 19, the best fit is n = 1.46(sup +0.39 sub -0.44). Comparing the COBE DMR delta-T/T at small l to the delta-T/T at l approximately = 50 from degree scale anisotropy experiments gives a smaller range of acceptable spectral indices which includes n = 1.
This paper highlights the accomplishments on a joint effort between NASA - Marshall Space Flight Center and Texas A and M University to develop accurate seal analysis software for use in rocket turbopump design, design audits and trouble shooting. Results for arbitrary clearance profile, transient simulation, thermal effects solution and flexible seal wall model are presented. A new solution for eccentric seals based on cubic spline interpolation and ordinary differential equation integration is also presented.
Remote sensing of the atmosphere from high earth orbit is very attractive due to the large field of view obtained and a true global perspective. This viewpoint is complicated by earth curvature effects so that slant path enhancement and absorption effects, small from low earth orbit, become dominant even at small nadir view angles. The effect is further complicated by the large range of local times and solar zenith angles in a single image leading to a modulation of the image intensity by a significant portion of the diurnal height variation of the absorbing layer. The latter effect is significant in particular for mesospheric, stratospheric and auroral emissions due to their depth in the atmosphere. As a particular case, the emissions from atomic oxygen (130.4 and 135.6 nm) and molecular nitrogen (two LBH bands, LBHS from 140 to 160 nm and LBHL from 160 to 180 nm) as viewed from the Ultraviolet Imager (UVI) are examined. The LBH emissions are of particular interest since LBHS has significant 02 absorption while LBHL does not, In the case of auroral emissions this differential absorption, well examined in the nadir, gives information about the height of the emission and therefore the energy of the precipitating particles. Using simulations of the viewing geometry and images from the UVI we examine these effects and obtain correction factors to adjust to the nadir case with a significant improvement of the derived characteristic energy. There is a surprisingly large effect on the images from the 02 diurnal layer height changes. An empirical compensation to the nadir case is explored based on the local nadir and local zenith angles for each portion of the image. These compensations are demonstrated as applied to the above emissions in both auroral and dayglow images and compared to models. The extension of these findings to other instruments, emissions and spectral regions is examined.
Electric potentials and magnetic fields generated by ensembles of synchronously active neurons in response to external stimuli provide information essential to understanding the processes underlying cognitive and sensorimotor activity. Interpreting recordings of these potentials and fields is difficult as each detector records signals simultaneously generated by various regions throughout the brain. We introduce the differentially Variable Component Analysis (dVCA) algorithm, which relies on trial-to-trial variability in response amplitude and latency to identify multiple components. Using simulations we evaluate the importance of response variability to component identification, the robustness of dVCA to noise, and its ability to characterize single-trial data. Finally, we evaluate the technique using visually evoked field potentials recorded at incremental depths across the layers of cortical area VI, in an awake, behaving macaque monkey.
This paper presents a receding horizon model predictive control variation of the combined Bernstein polynomial optimal reciprocal collision avoidance (ORCA) differential dynamic programming (COBRA-DDP) algorithm for AAM vehicles. Collision avoidance in combination with effective trajectory replanning are expected to be core components of AAM vehicles operating within a crowded airspace. This environment necessitates the use of real-time trajectory planning algorithms that are capable of planning around large amounts of stationary and moving obstacles. Previous work on COBRA-DDP demonstrated the capability of the algorithm to produce dynamically feasible trajectories for AAM vehicles and general collision avoidance. This paper improves upon the previous work by increasing the number of stationary and moving obstacles, implementing a variation of COBRA-DDP that lends itself to real-time application. These advancements are demonstrated on a vertical takeoff and landing (VTOL) vehicle simulation with highly nonlinear vehicle dynamics.
This paper presents a receding horizon model predictive control variation of the combined Bernstein polynomial optimal reciprocal collision avoidance (ORCA) differential dynamic programming (COBRA-DDP) algorithm for AAM vehicles. Collision avoidance in combination with effective trajectory replanning are expected to be core components of AAM vehicles operating within a crowded airspace. This environment necessitates the use of real-time trajectory planning algorithms that are capable of planning around large amounts of stationary and moving obstacles. Previous work on COBRA-DDP demonstrated the capability of the algorithm to produce dynamically feasible trajectories for AAM vehicles and general collision avoidance. This paper improves upon the previous work by increasing the number of stationary and moving obstacles, implementing a variation of COBRA-DDP that lends itself to real-time application. These advancements are demonstrated on a vertical takeoff and landing (VTOL) vehicle simulation with highly nonlinear vehicle dynamics.
Little is known about the effects of spaceflight conditions on epigenetics. The term epigenetics describes changes to the genome that can affect expression of a gene without changes to the sequence of DNA. Epigenetic processes are thought to underlie cellular differentiation, where transcription of specific genes occurs in response to key stimuli, and may be heritable - passing from one cell to its daughter cell. We hypothesize that the mechanical environment during spaceflight, namely microgravity-induced weightlessness or exercise regulate gene expression in the osteoblast-lineage cells both to control bone formation by osteoblasts and bone resorption by osteoclasts, which continually shapes bone structure throughout life. Similarly we intend to evaluate how radiation regulates these same bone cell activity and differentiation related genes. We further hypothesize that the regulation in bone cell gene expression is at least partially controlled through epigenetic mechanisms of methylation or small non-coding RNA (microRNAs). We have acquired preliminary data suggesting that global genome methylation is modified in response to axial compression of the tibia - a model of exercise. We intend to pursue these hypotheses wherein we will evaluate changes in gene expression and, congruently, changes in epigenetic state in bones from mice subjected to the aforementioned conditions: hindlimb unloading to simulate weightlessness, axial compression of the tibia, or radiation exposure in order to gain insight into the role of epigenetics in spaceflight-induced bone loss.
Abstract Among planetary dynamos, the magnetic field of Saturn stands out in its exceptional level of axisymmetry. One of its peculiar features is that the magnetic dipole mode is tilted with respect to the planetary rotation axis by only ≈0.007° or less. Numerical dynamo simulations performed in this context have had great difficulty in producing such small dipole tilt angles without introducing ad hoc ingredients such as a latitudinally varying heat flux pattern in the outer layers or stably stratified layers (SSLs). Here we present a numerical dynamo simulation that generates a highly axisymmetric dynamo with a dipole tilt of about ≈0.0008° on average. The model consists of a deep dynamo layer and an overlying low-conductivity layer but without any SSLs. We highlight a novel mechanism where strong differential rotation generated in the atmospheric layer penetrates into the dynamo region, helping to maintain a very small magnetic dipole tilt. Plain Language Summary Saturn's dipole-dominant magnetic field exhibits a very peculiar feature: the dipole component of the planetary magnetic field is tilted by less than ≈0.007° with respect to the planetary spin axis. Numerical simulations performed in this context suggest that if a spatial heat-flux variation is imposed, along with a stably stratified region, on top of an active dynamo layer, then small dipole tilt values can be realized. Here we present a model where extremely small dipole tilt values can be achieved without these ad hoc ingredients. Our simulations demonstrate that dynamo theory allows extremely small dipole tilt values in a relatively simple model configuration.
To assess prenatal and postnatal developmental status in the offspring of a group of animals, it is typical to examine fetuses from some of the dams as well as infants born to the remaining dams. Statistical limitations often arise, particularly when the animals are rare or especially precious, because all offspring of the dam represent only a single statistical observation; littermates are not independent observations (biologically or statistically). We describe a study in which pregnant laboratory rats were laparotomized on day 7 of gestation (GD7) to ascertain the number and distribution of uterine implantation sites and were subjected to a simulated experience on a 10-day space shuttle flight. After the simulated landing on GD18, rats were unilaterally hysterectomized, thus providing a sample of fetuses from 10 independent uteruses, followed by successful vaginal delivery on GD22, yielding postnatal samples from 10 uteruses. A broad profile of maternal and offspring morphologic and physiologic measures indicated that these novel sampling procedures did not compromise maternal well-being and maintained normal offspring development and function. Measures included maternal organ weights and hormone concentrations, offspring body size, growth, organ weights, sexual differentiation, and catecholamine concentrations.
A Second Moment Closure computation (SMC) is compared in detail with the Direct Numerical Simulation (DNS) data of Le and Moin for the backstep flow at Re = 5,000 in an attempt to understand why the intensity of the backflow and, consequently, the friction coefficient in the recirculation bubble are severely underestimated. The data show that this recirculation bubble is far from being laminar except in the very near wall layer. A novel 'differential a priori' procedure was used, in which the full transport equation for one isolated component of the Reynolds stress tensor was solved using DNS data as input. Conclusions are then different from what would have been deduced by comparing a full simulation to a DNS. One cause of discrepancy was traced back to insufficient transfer of energy to the normal stress by pressure strain, but was not cured. A significant finding, confirmed by the DNS data in the core region of a channel flow, is that the coefficient that controls destruction of dissipation, C epsilon(sub 2), should be decreased by a factor of 2 when production is vanishing. This is also the case in the recirculation bubble, and a new formulation has cured 25% of the backflow discrepancy.
Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.