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.
Reactor powered Brayton cycle power conversion system, evaluating performance at different turbine inlet temperatures and power levels by computer simulation
There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of freestream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Reex=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.
The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.
Lunar rocks may be classified into three major groups: (1) coarse-grained igneous rocks, (2) fine-grained igneous rocks, and (3) breccias. Group 1 is interpreted as primitive lunar crustal rocks that display various degrees of crushing and/or annealing. Group 2 is interpreted as volcanic rocks. Group 3 is interpreted as resulting from impacts on the lunar surface and is subdivided on the basis of matrix textures into fragmental breccias, crystalline breccias that have been annealed, and crystalline breccias with igneous matrices. A synthesis of the data concerning lunar highlands polymict breccias compels the prediction that the breccias should have homogeneous matrices from rock to rock within regions of the highlands of limited size where impact mixing has been efficient and extensive. But the returned breccias, even from one landing site, display a wide range in composition. This incompatibility between prediction and observation is a paradox that may be resolved by a process that acts after impact mixing to cause a differentiation of the breccia compositions. Partial melting of the local average crustal composition (as modeled by the average soil composition for each site) and separation of melt and residue in ejecta and/or fall-back blankets are compatible with the reviewed data and may resolve the paradox.
There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of free-stream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Re(ex)=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.
Design features, test data and projected performance levels of the telescope/sun sensor assembly for the Halogen Occultation Experiment (HALOE) on board the Upper Atmosphere Research Satellite are described. HALOE will gather data on stratospheric chemical species with an IR telescope pointed through the atmosphere at the sun during occultation periods. The pointing accuracy will need to be 1 arcmin elevation and 1.6 arcmin azimuthal. Adjustments of the gimbaled instrument are guided by a sun sensor and its associated electronics. The components were subjected to thermal-vacuum, optical boresight stability tests and the data generated were compared with a finite element model of the telescope and sun sensor. The tests consisted of exposure to various thermal gradients while the telescope and sun sensor alignments were tracked by optical deflectometry. The thermal behavior was compared with predictions made with a finite difference model. Alignments were within tolerable ranges and the thermal behavior model was concluded valid for predicting the thermal behavior of orbiting instruments.
The sound generated by viscous flow past a circular cylinder is predicted via the Lighthill acoustic analogy approach. The two dimensional flow field is predicted using two unsteady Reynolds-averaged Navier-Stokes solvers. Flow field computations are made for laminar flow at three Reynolds numbers (Re = 1000, Re = 10,000, and Re = 90,000) and two different turbulent models at Re = 90,000. The unsteady surface pressures are utilized by an acoustics code that implements Farassat's formulation 1A to predict the acoustic field. The acoustic code is a 3-D code - 2-D results are found by using a long cylinder length. The 2-D predictions overpredict the acoustic amplitude; however, if correlation lengths in the range of 3 to 10 cylinder diameters are used, the predicted acoustic amplitude agrees well with experiment.
An acoustic analogy is developed based on the Euler equations for broadband shock- associated noise (BBSAN) that directly incorporates the vector Green's function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) as the mean flow. The vector Green's function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation. An adjoint vector Green's function solver is implemented to determine the vector Green's function based on a locally parallel mean flow at streamwise locations of the SRANS solution. However, the developed acoustic analogy could easily be based on any adjoint vector Green's function solver, such as one that makes no assumptions about the mean flow. The newly developed acoustic analogy can be simplified to one that uses the Green's function associated with the Helmholtz equation, which is consistent with the formulation of Morris and Miller (AIAAJ 2010). A large number of predictions are generated using three different nozzles over a wide range of fully expanded Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise labs. In addition, two models for the so-called 'fine-scale' mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include the propagation effects, especially in the upstream direction of the jet.
Predictions of solar activity are an essential part of our Space Weather forecast capability. Users are requiring usable predictions of an upcoming solar cycle to be delivered several years before solar minimum. A set of predictions of the amplitude of Solar Cycle 24 accumulated in 2008 ranged from zero to unprecedented levels of solar activity. The predictions formed an almost normal distribution, centered on the average amplitude of all preceding solar cycles. The average of the current compilation of 105 predictions of the annual-average sunspot number is 106 +/- 31, slightly lower than earlier compilations but still with a wide distribution. Solar Cycle 24 is on track to have a below-average amplitude, peaking at an annual sunspot number of about 80. Our need for solar activity predictions and our desire for those predictions to be made ever earlier in the preceding solar cycle will be discussed. Solar Cycle 24 has been a below-average sunspot cycle. There were peaks in the daily and monthly averaged sunspot number in the Northern Hemisphere in 2011 and in the Southern Hemisphere in 2014. With the rapid increase in solar data and capability of numerical models of the solar convection zone we are developing the ability to forecast the level of the next sunspot cycle. But predictions based only on the statistics of the sunspot number are not adequate for predicting the next solar maximum. I will describe how we did in predicting the amplitude of Solar Cycle 24 and describe how solar polar field predictions could be made more accurate in the future.
A synthesis of the relevant data concerning lunar highland polymict breccias from the fields of petrography, chemistry, photogeology, and impact studies compels the prediction that the breccias should have homogeneous matrices from rock to rock within regions of the highlands of limited size where impact mixing has been efficient and extensive. But the matrices of returned breccias, even from one landing site, display a wide range in composition. This incompatibility between prediction and observation is a paradox that may be resolved by a process that acts after impact mixing to cause a differentiation of the breccia compositions. Partial melting of the local surficial material, and separation of melt and residue in ejecta and/or fall-back blankets, is compatible with the reviewed data and may resolve the paradox.
The ultraviolet spectrum of the magnetic degenerate Grw -70 deg 8247 contains a relatively sharp, but resolved absorption feature near 1345 A. This is interpreted as the sigma(-), 1s0 to 2p-1 component of hydrogen Ly-alpha in a very strong magnetic field, near 350 megagauss. A cyclotron absorption may also exist near 2300 A, corresponding to such a field. Review of theoretical predictions of these hydrogen levels shows that a relatively sharp feature near 1335 A would be produced in fields within the range from 280 to 800 megagauss. The predicted wavelength has considerable uncertainty, perhaps + or - 10 A. Diffuse features in the red region of Palomar CCD spectra may arise from two sigma(-) and a pi component, based on even less certain computations for H-alpha. The large field is consistent with polarization data. The original bands discovered by Minkowski are not yet identifiable. But this important, magnetic star contains hydrogen even on the basis of Ly-alpha alone.
The original research on the use on nonlinear vibration technique to solve for the hypervelocity ballistic limit for double plates is examined. Such structure is commonly found in typical Space Station design where the incoming space or man-made debris would be fragmented upon hitting the outer plate (shield) and the subsequent impact on the main wall would result in a much reduced damge of the space station or spacecraft. The existing few theoretical impact equations do not agree well with each other. The existing computer code bumper used at NASA-Johnson Space Center appears to predict an unconservative ballistic limit when compared with experimental data where the velocity ranges from 3 km/s to 8 km/s. Such unconservative prediction is unacceptable from a practical safe design point of view. The bumper code is based on Wilkinson's (1968) paper and his equations have not been improved nor modified even though they are viewed with suspicion due to lack of agreement with experiments. To make matters worse, there is not another theory which is better than Wilkinson's equation and the designers are forced to use purely empirical Nysmith (1969) or semiempirical equations developed by Cour-Palais in 1969. The Cour-Palais equations were later modified empirically in 1989. The purpose of the present investigation is to examine the many assumptions of the Wilkinson equation. An attempt is made to present design charts based on the modified-Wilkinson equation so that the designer can get a feel of the ranges of the parameters which are of interest and discard a huge range of parameters, thus, significantly reducing the number of test shots required. The analysis is based on a solution of the governing nonlinear differential equations for a plate, assuming axisymmetric behavior using polar coordinates.
This paper reports an analysis and prediction development of rotor broadband noise. The two primary components of this noise are Blade-Wake Interaction (BWI) noise, due to the blades' interaction with the turbulent wakes of the preceding blades, and "Self" noise, due to the development and shedding of turbulence within the blades' boundary layers. Emphasized in this report is the new code development for Self noise. The analysis and validation employs data from the HART program, a model BO-105 rotor wind tunnel test conducted in the German-Dutch Wind Tunnel (DNW). The BWI noise predictions are based on measured pressure response coherence functions using cross-spectral methods. The Self noise predictions are based on previously reported semiempirical modeling of Self noise obtained from isolated airfoil sections and the use of CAMRAD.Modl to define rotor performance and local blade segment flow conditions. Both BWI and Self noise from individual blade segments are Doppler shifted and summed at the observer positions. Prediction comparisons with measurements show good agreement for a range of rotor operating conditions from climb to steep descent. The broadband noise predictions, along with those of harmonic and impulsive Blade-Vortex Interaction (BVI) noise predictions, demonstrate a significant advance in predictive capability for main rotor noise.
Laser Ranging Simulation Program (LRSP) is a computer program that predicts selected aspects of the performances of a laser altimeter or other laser ranging or remote-sensing systems and is especially applicable to a laser-based system used to map terrain from a distance of several kilometers. Designed to run in a more recent version (5 or higher) of the MATLAB programming language, LRSP exploits the numerical and graphical capabilities of MATLAB. LRSP generates a graphical user interface that includes a pop-up menu that prompts the user for the input of data that determine the performance of a laser ranging system. Examples of input data include duration and energy of the laser pulse, the laser wavelength, the width of the laser beam, and several parameters that characterize the transmitting and receiving optics, the receiving electronic circuitry, and the optical properties of the atmosphere and the terrain. When the input data have been entered, LRSP computes the signal-to-noise ratio as a function of range, signal and noise currents, and ranging and pointing errors.
We use Atacama Cosmology Telescope (ACT) observations at 98 GHz (2015–2019), 150 GHz (2013–2019), and 229 GHz (2017–2019) to perform a blind shift-and-stack search for Planet 9. The search explores distances from 300 au to 2000 au and velocities up to 6\farcm3 per year, depending on the distance (r). For a 5 Earth-mass Planet 9 the detection limit varies from 325 au to 625 au, depending on the sky location. For a 10 Earth-mass planet the corresponding range is 425 au to 775 au. The predicted aphelion and most likely location of the planet corresponds to the shallower end of these ranges. The search covers the whole 18,000 square degrees of the ACT survey. No significant detections are found, which is used to place limits on the millimeter-wave flux density of Planet 9 over much of its orbit. Overall we eliminate roughly 17% and 9% of the parameter space for a 5 and 10 Earth-mass Planet 9, respectively. These bounds approach those of a recent INPOP19a ephemeris-based analysis, but do not exceed it. We also provide a list of the 10 strongest candidates from the search for possible follow-up. More generally, we exclude (at 95% confidence) the presence of an unknown solar system object within our survey area brighter than 4–12 mJy (depending on position) at 150 GHz with current distance 300 au < r < 600 au and heliocentric angular velocity 1\farcm5 per yr < v x (500au/r) < 2\farcs3 per yr, corresponding to low-to-moderate eccentricities. These limits worsen gradually beyond 600 au, reaching 5–15 mJy by 1500 au.
As part of an effort to determine the feasibility of helium transfer in space, a centrifugal pump was tested in He II at a variety of flow rates, pump speeds, and fluid temperatures. The pump, which has a straight bladed impeller 6.86 cm in diameter, generated a maximum pressure rise of 15 kPa and a maximum flow rate of 22 g/s for the conditions of the test. Pump performance seems to be independent of fluid temperature and is in good agreement with the values predicted by the manufacturer. Over the range of flow coefficients, the measured maximum efficiency is around 50 percent. Cavitation is observed in the pump and is thought to be highly dependent on the local heating of the helium in the pump. Preliminary measurements of the noise spectra of the pump suggest a possible mechanism to predict the onset of cavitation.
Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA GRC. The TDU consists of three subsystems: the Reactor Simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the Heat Exchanger Manifold (HXM). An Annular Linear Induction Pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now), is referred to as the RxSim subsystem and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the CAD model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras turbulence model), and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is -0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k-epsilon turbulence model) when compared to measured data.
Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA Glenn Research Center. The TDU consists of three subsystems: the reactor simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the heat exchanger manifold (HXM). An annular linear induction pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now) is referred to as the "RxSim subsystem" and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the computer-aided-design (CAD) model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras (S‒A) turbulence model) and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is ‒0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k- turbulence model) when compared to measured data.