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.
SUBAERF2 program developed to provide for analysis of subsonic aerodynamics and design by iteration of low-speed wing/flap systems. Based on linearized-theory lifting-surface solution but also accounts for some nonlinear characteristics. Effects of leading- and trailing-edge flaps included. New and improved version of programs described in LAR-13116 and LAR-12987, and replaces both of them. Written in FORTRAN V.
Detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity (EPRV) measurements. The required precision of these measurements necessitates developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model was used to determine the statistical properties of the turbulent plasma and to obtain synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations. Presented are the initial results of our effort to characterize stellar jitter properties as well as future objectives.
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code, are used with a blowing boundary condition in the DPLR hypersonic CFD code. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver, is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
A number of modifications have been proposed in order to improve the jet noise prediction capabilities of the MGB code. This code which was developed at General Electric, employees the concept of acoustic analogy for the prediction of turbulent mixing noise. The source convection and also refraction of sound due to the shrouding effect of the mean flow are accounted for by incorporating the high frequency solution to Lilley's equation for cylindrical jets (Balsa and Mani). The broadband shock-associated noise is estimated using Harper-Bourne and Fisher's shock noise theory. The proposed modifications are aimed at improving the aerodynamic predictions (source/spectrum computations) and allowing for the non- axisymmetric effects in the jet plume and nozzle geometry (sound/flow interaction). In addition, recent advances in shock noise prediction as proposed by Tam can be employed to predict the shock-associated noise as an addition to the jet mixing noise when the flow is not perfectly expanded. Here we concentrate on the aerodynamic predictions using the PARC code with a k-E turbulence model and the ensuing turbulent mixing noise. The geometry under consideration is an axisymmetric convergent-divergent nozzle at its design operating conditions. Aerodynamic and acoustic computations are compared with data as well as predictions due to the original MGB model using Reichardt's aerodynamic theory.
A description of the CELCAP cogeneration analysis program is presented. A detailed description of the methodology used by the Naval Civil Engineering Laboratory in developing the CELCAP code and the procedures for analyzing cogeneration systems for a given user are given. The four engines modeled in CELCAP are: gas turbine with exhaust heat boiler, diesel engine with waste heat boiler, single automatic-extraction steam turbine, and back-pressure steam turbine. Both the design point and part-load performances are taken into account in the engine models. The load model describes how the hourly electric and steam demand of the user is represented by 24 hourly profiles. The economic model describes how the annual and life-cycle operating costs that include the costs of fuel, purchased electricity, and operation and maintenance of engines and boilers are calculated. The CELCAP code structure and principal functions of the code are described to how the various components of the code are related to each other. Three examples of the application of the CELCAP code are given to illustrate the versatility of the code. The examples shown represent cases of system selection, system modification, and system optimization.
Satellite-based solar resource data are often developed and validated by using binary cloudiness categories: clear sky or overcast cloudy sky. To investigate the reliability of solar resource data in partially cloudy conditions, we estimate cloud fraction using two distinct algorithms: a physical retrieval model using surface observed global horizontal irradiance (GHI) and direct normal irradiance (DNI) and a temporal average of cloud mask data estimated by the observed DNI. Our analysis reveals a significant presence of scattered clouds, broken clouds, and mismatches between satellite- and surface-based cloud data at 17 surface sites across the contiguous United States, though confidently clear and cloudy conditions collectively account for more than 70 % of the data. Solar radiation is computed using the National Solar Radiation Database (NSRDB) algorithm and validated using surface observations. Here, our findings suggest that, in the presence of scattered clouds, NSRDB data for clear-sky conditions can be subject to significant overestimation. In cloudy-sky conditions classified by satellite data, DNI computed by the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) can be underestimated when limited clouds are detected by surface observations. The bias observed in several cloudiness categories indicates that the NSRDB is exceptionally accurate in confidently clear conditions. However, clear-sky conditions with scattered clouds and mismatched cloud data contribute significantly to the overall uncertainties in the NSRDB. Therefore, future improvements in solar resource data should involve development and implementation of satellite-derived cloud fraction and should consider a novel radiative transfer model accounting for amplified cloud reflection. The evaluation within cloudiness categories also provides a physical rationale for the superior performance of FARMS-DNI compared to the Direct Insolation Simulation Code (DISC) in both cloudy-sky and all-sky conditions.
Here, this paper presents the development of a lattice physics–to–core simulator two-step procedure for the rapid analysis of the Advanced High Temperature Reactor (AHTR). Lattice physics, reflector, and control blade models were developed from which cross-section libraries could be generated for a nodal core simulator. Few-group structures for the core simulator were also generated to account for the neutronic characteristics of AHTR. After developing the AHTR two-step procedure, cross-section libraries were generated using the SERPENT continuous-energy Monte Carlo code. These libraries were then used in the core simulator NESTLE to perform full-core calculations, which were in turn benchmarked against reference SERPENT full-core models. Benchmarking results showed reasonable accuracy of the developed two-step procedure but revealed an inherent inadequacy in the one-dimensional radial reflector model and showed a likely need for a greater number of energy groups than were used in this study.
The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.
Upwind algorithms are developed for the numerical solution of the multidimensional Euler equations for real gases. Flux-splitting methods are derived which account for a general equation of state. Approximations to the state equation based on physical arguments result in simplified algorithms which may be implemented into existing perfect-gas codes. Applications of the method to several high-Mach-number high-temperature flows are presented for two and three space dimensions.
A new rotating stellar evolution code is developed and applied to the sun. A hydrostatic fully convective premain-sequence model is evolved to the age of the sun. As the model evolves, it accounts for angular-momentum loss via a magnetic wind and angular-momentum redistribution by rotationally induced instabilities. The resulting models have an oblateness in agreement with observed upper limits. The rotation curves show two main features: the outer layers exhibit minimal radial differential rotation, and a rapidly rotating central core is preserved. These basic features persist through a wide range of model parameters.
Further improvements made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods by accounting for condensed phase, thick boundary layer and free stream turbulence effects are discussed. The condensed phase effects were included through species composition effect considered via input to the code and through particle damping effect considered via a turbulence model. The thrust loss calculation procedure for thick boundary layer effects was improved and the optimization of net thrust with respect to nozzle length was performed. The effects of free stream turbulence were approximately modeled in the turbulence model.
Existing and proposed methods for the prediction of plume radiation are discussed in terms of their application to the NASA Advanced Solid Rocket Motor (ASRM) and Space Shuttle Main Engine (SSME) projects. Extrapolations of the Solid Rocket Motor (SRM) are discussed with respect to preliminary predictions of the primary and secondary radiation environments. The methodology for radiation and initial plume property predictions are set forth, including a new code for scattering media and independent secondary source models based on flight data. The Monte Carlo code employs a reverse-evaluation approach which traces rays back to their point of absorption in the plume. The SRM sea-level plume model is modified to account for the increased radiation in the ASRM plume due to the ASRM's propellant chemistry. The ASRM cycle-1 environment predictions are shown to identify a potential reason for the shutdown spike identified with pre-SRM staging.
Using new XMM and Chandra observations, we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow (M = 4.5 +/- 0.2 solar mass/yr) models within the central approx.30 kpc. Alternatively, the data can be fitted equally well if the temperature within each spherical shell varies continuously from approx.T(sub h) to T(sub c) approx. T(sub h)/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion T(sub h) approaches T(sub c) required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multiphase gas having a limited temperature range. We do not find any evidence that azimuthal temperature variations within each annulus on the sky can account for the range in temperatures within each shell. We provide a detailed investigation of the systematic errors on the derived spectral models considering the effects of calibration, plasma codes, bandwidth, variable NH, and background rate. We find that the RGS gratings and the EPIC and ACIS CCDs give fully consistent results when the same models are fitted over the same energy ranges for each instrument. The cooler component of the 2T model has a temperature (T(sub c) approx. 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (T(sub h) approx. 1.4 keV) characteristic of the virial temperature of the solar mass halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R approx. equals 10 kpc, bubbles of gas heated to approx.T(sub h) in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R approx. 50 kpc, as suggested by the sharp edge in the EPIC images. Subject headings: cooling flows - galaxies: formation - galaxies: halos - intergalactic medium - X-rays: galaxies: clusters
A new model for the radiation environment to be found on the planet Mars due to Galactic Cosmic Rays (OCR) has been developed at the NASA Langley Research Center. Solar modulated primary particles rescaled for Mars conditions are transported through the Martian atmosphere, with temporal properties modeled with variable timescales, down to the surface, with altitude and backscattering patterns taken into account. The Martian atmosphere has been modeled by using the Mars Global Reference Atmospheric Model--version 2001 (Mars-GRAM 2001). The altitude to compute the atmospheric thickness profile has been determined by using a model for the topography based on the data provided by the Mars Orbiter Laser Altimeter (MOLA) instrument on board the Mars Global Surveyor (MGS) spacecraft. The Mars surface composition has been modeled based on averages over the measurements obtained from orbiting spacecraft and at various landing sites, taking into account the possible volatile inventory (e.g., CO2 ice, H2O ice) along with its time variation throughout the Martian year. Particle transport has been performed with the HZETRN heavy ion code. The Mars Radiation Environment Model has been made available worldwide through the Space Ionizing Radiation Effects and Shielding Tools (SIREST) website, a project of NASA Langley Research Center. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.
An assessment was made of the capability of jet noise prediction codes over a broad range of jet flows, with the objective of quantifying current capabilities and identifying areas requiring future research investment. Three separate codes in NASA s possession, representative of two classes of jet noise prediction codes, were evaluated, one empirical and two statistical. The empirical code is the Stone Jet Noise Module (ST2JET) contained within the ANOPP aircraft noise prediction code. It is well documented, and represents the state of the art in semi-empirical acoustic prediction codes where virtual sources are attributed to various aspects of noise generation in each jet. These sources, in combination, predict the spectral directivity of a jet plume. A total of 258 jet noise cases were examined on the ST2JET code, each run requiring only fractions of a second to complete. Two statistical jet noise prediction codes were also evaluated, JeNo v1, and Jet3D. Fewer cases were run for the statistical prediction methods because they require substantially more resources, typically a Reynolds-Averaged Navier-Stokes solution of the jet, volume integration of the source statistical models over the entire plume, and a numerical solution of the governing propagation equation within the jet. In the evaluation process, substantial justification of experimental datasets used in the evaluations was made. In the end, none of the current codes can predict jet noise within experimental uncertainty. The empirical code came within 2dB on a 1/3 octave spectral basis for a wide range of flows. The statistical code Jet3D was within experimental uncertainty at broadside angles for hot supersonic jets, but errors in peak frequency and amplitude put it out of experimental uncertainty at cooler, lower speed conditions. Jet3D did not predict changes in directivity in the downstream angles. The statistical code JeNo,v1 was within experimental uncertainty predicting noise from cold subsonic jets at all angles, but did not predict changes with heating of the jet and did not account for directivity changes at supersonic conditions. Shortcomings addressed here give direction for future work relevant to the statistical-based prediction methods. A full report will be released as a chapter in a NASA publication assessing the state of the art in aircraft noise prediction.
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases [1]. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code [2], are used with a blowing boundary condition in the DPLR hypersonic CFD code [3]. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver [4], is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
Companion acoustic liner attenuation predictions and flight data representative of certification conditions are analyzed to determine the validation quality and gain an understanding of current prediction shortcomings. These comparisons are limited to the inlet broadband component of fan noise. Two prediction methods are assessed: a traditional, semiempirical model, and a computational approach with a parabolic duct propagation code coupled with a Ffowcs Williams-Hawkings solver for free-field propagation. The semiempirical method is widely used for aircraft system noise predictions, but does not explicitly account for several physical parameters important for sound attenuation. The computational approach includes these effects, but requires a significant increase in cycle time. For the three-degree-of-freedom (3DOF) liner tested, both methods have a bias error less than 0.5 dB for cutback and takeoff power settings, and a much larger 3–4 dB bias for approach power setting. The standard deviation of each method varies between 1–2 dB for different conditions. The similar performance of both approaches illustrates that, in this instance, the additional complexity of the computational approach does not provide a clear improvement over the semiempirical method. Planned improvements to the numerical grid and source assumptions may result in reevaluation of this conclusion in later work.
This paper summarizes the underlying theory in the recently developed KiteAeroDyn (KiteAD) module, which is part of KiteFAST (KiteFAST), a numerical simulation tool that is currently being developed at the NREL. KiteFAST is intended for the simulation of airborne wind energy systems (e.g., energy kites). KiteAD implements an innovative vortex step method for the calculation of aerodynamic loads on all of the lifting surfaces of the kite, which is capable of accounting for two-dimensional viscous effects through the use of nonlinear airfoil polars. Considering nonlinear airfoil polar data allows for the modeling of flap or other moving-surface deflections, as well as for the modelling of inflow conditions with larger angles angles of attack. From preliminary verification against other existing codes that can solve for lift distribution on lifting surfaces, the results of the newly implemented method proved to be numerically robust and computationally inexpensive, thus ideal for aeroelastic design and analysis applications.