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.
We discuss new functionality of the spectral simulation code CLOUDY which allows the user to calculate grids with one or more initial parameters varied and formats the predicted spectra in the standard FITS format. These files can then be imported into the x-ray spectral analysis software XSPEC and used as theoretical models for observations. We present and verify a test case. Finally, we consider a few observations and discuss our results.
This slide presentation reviews the use of computational fluid dynamics in performing analysis of the space shuttle with particular reference to the return to flight analysis and other shuttle problems. Slides show a comparison of pressure coefficient with the shuttle ascent configuration between the wind tunnel test and the computed values. the evolution of the grid system for the space shuttle launch vehicle (SSLv) from the early 80's to one in 2004, the grid configuration of the bipod ramp redesign from the original design to the current configuration, charts with the computations showing solid rocket booster surface pressures from wind tunnel data, calculated over two grid systems (i.e., the original 14 grid system, and the enhanced 113 grid system), and the computed flight orbiter wing loads are compared with strain gage data on STS-50 during flight. The loss of STS-107 initiated an unprecedented review of all external environments. The current SSLV grid system of 600+ grids, 1.8 Million surface points and 95+ million volume points is shown. The inflight entry analyses is shown, and the use of Overset CFD as a key part to many external tank redesign and debris assessments is discussed. The work that still remains to be accomplished for future shuttle flights is discussed.
Projects: (1) Boeing Launch Abort Analysis My first project for the summer was analyzing the Boeing CCDev Vehicle's abort aerodynamics using an inviscid solver (CART3D). The goal of the project was to develop the grid and CFD inputs necessary to use CART3D as a quick tool for investigating loading trends at various points along abort trajectories. As a supplementary task, I analyzed a few cases and compared them to the aerodatabase from the last generation geometry. (2) My second project for the summer dealt with investigating how heating changes as the height of a protuberance on top of a flat plate changes. The goal of this investigation is to better understand how to properly model heating on and around a protuberance. This is one of the biggest challenges when designing a re ]entry vehicle because very small changes in the shape and conditions leading up to a protuberance, not to mention the protuberance geometry, will greatly impact the local heating.
Deep convection that penetrates the tropopause, referred to here as overshooting convection, is capable of lifting tropospheric air well into the stratosphere. In addition to water, these overshoots also transport various chemical species, affecting chemistry and radiation in the stratosphere. It is not currently known, however, how much transport is a result of this mechanism. To better understand overshooting convection, this study aims to characterize the durations of overshooting events. To achieve this, radar data from the Next Generation Weather Radar (NEXRAD) network is composited onto a three-dimensional grid at 5-min intervals. Overshoots are identified by comparing echo-top heights with tropopause estimates derived from ERA5 reanalysis data. These overshoots are linked in space from one analysis time to the next to form tracks. This process is performed for 12 four-day sample windows in the months May–August of 2017–19. Track characteristics such as duration, overshoot area, tropopause-relative altitude, and column-maximum reflectivity are investigated. Positive correlations are found between track duration and other track characteristics. Integrated track volume is found as a product of the overshoot area, depth, and duration, and provides a measure of the potential stratospheric impact of each track. Short-lived tracks are observed to contribute the most total integrated volume when considering track duration, while tracks that overshoot by 2–3 km show the largest contribution when considering overshoot depth. A diurnal cycle is observed, with peak track initiation around 1600–1700 local time. Track-mean duration peaks a few hours earlier, while track-mean area and tropopause-relative height peak a few hours later.
Satellite based altimetric data taken by GEOS-3, SEASAT, and GEOSAT over the Aral Sea are analyzed and a least squares collocation technique is used to predict the geoid undulations on a 0.25 by 0.25 degree grid and to transform these geoid undulations to free air gravity anomalies. Rapp's 180 by 180 geopotential model is taken as the empirical model. The collocation procedure is performed with a set of local residual empirical covariance functions. For comparison, Rapp's global covariance functions and Jordan's self-consistent theoretical covariance functions based on Jordan's formulation and on locally derived parameters are also used to grid geoid undulations and to predict gravity anomalies. The sensitivity of the collocation results to the choice of covariance functions is discussed.
A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).
This paper presents a data analysis and modeling framework to tailor and develop linear parameter-varying (LPV) aeroservoelastic (ASE) model database for flexible aircrafts in broad 2D flight parameter space. The Kriging surrogate model is constructed using ASE models at a fraction of grid points within the original model database, and then the ASE model at any flight condition can be obtained simply through surrogate model interpolation. The greedy sampling algorithm is developed to select the next sample point that carries the worst relative error between the surrogate model prediction and the benchmark model in the frequency domain among all input-output channels. The process is iterated to incrementally improve surrogate model accuracy till a pre-determined tolerance or iteration budget is met. The methodology is applied to the ASE model database of a flexible aircraft currently being tested at NASA/AFRC for flutter suppression and gust load alleviation. Our studies indicate that the proposed method can reduce the number of models in the original database by 67%. Even so the ASE models obtained through Kriging interpolation match the model in the original database constructed directly from the physics-based tool with the worst relative error far below 1%. The interpolated ASE model exhibits continuously-varying gains along a set of prescribed flight conditions. More importantly, the selected grid points are distributed non-uniformly in the parameter space, a) capturing the distinctly different dynamic behavior and its dependence on flight parameters, and b) reiterating the need and utility for adaptive space sampling techniques for ASE model database compaction. The present framework is directly extendible to high-dimensional flight parameter space, and can be used to guide the ASE model development, model order reduction, robust control synthesis and novel vehicle design of flexible aircraft.
A highly efficient computer analysis was developed for predicting transonic nacelle/inlet flowfields. This algorithm can compute the three dimensional transonic flowfield about axisymmetric (or asymmetric) nacelle/inlet configurations at zero or nonzero incidence. The flowfield is determined by solving the full-potential equation in conservative form on a body-fitted curvilinear computational mesh. The difference equations are solved using the AF2 approximate factorization scheme. This report presents a discussion of the computational methods used to both generate the body-fitted curvilinear mesh and to obtain the inviscid flow solution. Computed results and correlations with existing methods and experiment are presented. Also presented are discussions on the organization of the grid generation (NGRIDA) computer program and the flow solution (NACELLE) computer program, descriptions of the respective subroutines, definitions of the required input parameters for both algorithms, a brief discussion on interpretation of the output, and sample cases to illustrate application of the analysis.
The In-Space Propulsion Technology (ISPT) Project, funded by NASA s Science Mission Directorate (SMD), is continuing to invest in propulsion technologies that will enable or enhance NASA robotic science missions. This overview provides development status, near-term mission benefits, applicability, and availability of in-space propulsion technologies in the areas of aerocapture, electric propulsion, advanced chemical thrusters, and systems analysis tools. Aerocapture investments improved: guidance, navigation, and control models of blunt-body rigid aeroshells; atmospheric models for Earth, Titan, Mars, and Venus; and models for aerothermal effects. Investments in electric propulsion technologies focused on completing NASA s Evolutionary Xenon Thruster (NEXT) ion propulsion system, a 0.6 to 7 kW throttle-able gridded ion system. The project is also concluding its High Voltage Hall Accelerator (HiVHAC) mid-term product specifically designed for a low-cost electric propulsion option. The primary chemical propulsion investment is on the high-temperature Advanced Material Bipropellant Rocket (AMBR) engine providing higher performance for lower cost. The project is also delivering products to assist technology infusion and quantify mission applicability and benefits through mission analysis and tools. In-space propulsion technologies are applicable, and potentially enabling for flagship destinations currently under evaluation, as well as having broad applicability to future Discovery and New Frontiers mission solicitations.
During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.
Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.
Data from SCR-B on Nimbus 5 have been processed to yield global, orbital temperatures at 10, 5, 2, 1, and 0.4 mb for the period January 1977 through April 1978 under the current task. In addition gridded values at 10 deg latitude by 20 deg longitude were prepared by space-time interpolation for the period January 1975 through April 1978. Temperature retrieval was based on regression of radiances against Meteorological Rocket Network data, with regressions recomputed at approximately six-month intervals. This data now completes a consistent time series from April 1970 to April 1978 for all available radiance data from SCR A and SCR B on Nimbus 4 and 5. The processing details for the current period are discussed, but is also applicable to the previous data periods. The accuracy of the temperature retrievals for each 6-month period for the entire eight years is given in the Appendices.
The In-Space Propulsion Technology (ISPT) Project, funded by NASA's Science Mission Directorate (SMD), is continuing to invest in propulsion technologies that will enable or enhance NASA robotic science missions. This overview provides development status, near-term mission benefits, applicability, and availability of in-space propulsion technologies in the areas of aerocapture, electric propulsion, advanced chemical thrusters, and systems analysis tools. Aerocapture investments improved (1) guidance, navigation, and control models of blunt-body rigid aeroshells, 2) atmospheric models for Earth, Titan, Mars and Venus, and 3) models for aerothermal effects. Investments in electric propulsion technologies focused on completing NASA s Evolutionary Xenon Thruster (NEXT) ion propulsion system, a 0.6-7 kW throttle-able gridded ion system. The project is also concluding its High Voltage Hall Accelerator (HiVHAC) mid-term product specifically designed for a low-cost electric propulsion option. The primary chemical propulsion investment is on the high-temperature Advanced Material Bipropellant Rocket (AMBR) engine providing higher performance for lower cost. The project is also delivering products to assist technology infusion and quantify mission applicability and benefits through mission analysis and tools. In-space propulsion technologies are applicable, and potentially enabling for flagship destinations currently under evaluation, as well as having broad applicability to future Discovery and New Frontiers mission solicitations.
The In-Space Propulsion Technology (ISPT) program is developing spacecraft bus and platform technologies that will enable or enhance NASA robotic science missions. The ISPT program is currently developing technology in four areas that include Propulsion System Technologies (electric and chemical), Entry Vehicle Technologies (aerocapture and Earth entry vehicles), Spacecraft Bus and Sample Return Propulsion Technologies (components and ascent vehicles), and Systems/Mission Analysis. Three technologies are ready for near-term flight infusion: 1) the high-temperature Advanced Material Bipropellant Rocket (AMBR) engine providing higher performance 2) NASAs Evolutionary Xenon Thruster (NEXT) ion propulsion system, a 0.6-7 kW throttle-able gridded ion system and 3) Aerocapture technology development with investments in a family of thermal protection system (TPS) materials and structures guidance, navigation, and control (GN&C) models of blunt-body rigid aeroshells and aerothermal effect models. Two component technologies being developed with flight infusion in mind are the Advanced Xenon Flow Control System, and ultra-lightweight propellant tank technologies. Future direction for ISPT are technologies that relate to sample return missions and other spacecraft bus technology needs like: 1) Mars Ascent Vehicles (MAV) 2) multi-mission technologies for Earth Entry Vehicles (MMEEV) and 3) electric propulsion. These technologies are more vehicle and mission-focused, and present a different set of technology development and infusion steps beyond those previously implemented. The Systems/Mission Analysis area is focused on developing tools and assessing the application of propulsion and spacecraft bus technologies to a wide variety of mission concepts. These in-space propulsion technologies are applicable, and potentially enabling for future NASA Discovery, New Frontiers, and sample return missions currently under consideration, as well as having broad applicability to potential Flagship missions. This paper provides a brief overview of the ISPT program, describing the development status and technology infusion readiness of in-space propulsion technologies in the areas of electric propulsion, Aerocapture, Earth entry vehicles, propulsion components, Mars ascent vehicle, and mission/systems analysis.
New analysis tools and additional unanalyzed observations justify a reanalysis of MR Cygni. The reanalysis applied successively more restrictive physical models, each with an optimization program. The final model assigned separate first and second order limb darkening coefficients, from model atmospheres, to individual grid points. Proper operation of the optimization procedure was tested on simulated observational data, produced by light synthesis with assigned system parameters, and modulated by simulated observational error. The iterative solution converged to a weakly-determined mass ratio of 0.75. Assuming the B3 primary component is on the main sequence, the HR diagram location of the secondary from the light ratio (ordinate) and adjusted T sub eff (abscissa) was calculated. The derived mass ratio, together with a main-sequence mass for the B3 component, implies a main-sequence secondary spectral type of B4. The photometrically-determined secondary radii agree with this spectral type, in marginal disagreement with the B7 type from the HR diagram analysis. The individual masses, derived from the radial velocity curve of the primary component, the photometrically-determined i, and alternative values of derived mass ratio are seriously discrepant with main sequence objects. The imputed physical status of the system is in disagreement with representations that have appeared in the literature.
A steady-state analysis is conducted to examine the basic flow structure of a non-Newtonian fluid in a domain including an inflow region, a contraction region, and an outflow region. A Cartesian grid system is used throughout the entire flow domain, including the contraction region, thus creating an irregular grid cell structure adjacent to the curved boundary. At node points adjacent to the curved boundary symmetry conditions are derived for the different flow variables in order to solve the governing difference equations. Attention is given to the motion and non-Newtonian constitutive equations, the boundary conditions, the numerical modeling of the non-Newtonian equations, the stream function contour lines for the non-Newtonian fluid, the vorticity contour lines for the non-Newtonian fluid, the velocity profile across the contraction, and the shear stress contour lines for the non-Newtonian fluid.
Work done on evaluating the geometric and radiometric quality of early LANDSAT-4 sensor data is described. Band to band and channel to channel registration evaluations were carried out using a line correlator. Visual blink comparisons were run on an image display to observe band to band registration over 512 x 512 pixel blocks. The results indicate a .5 pixel line misregistration between the 1.55 to 1.75, 2.08 to 2.35 micrometer bands and the first four bands. Also a four 30M line and column misregistration of the thermal IR band was observed. Radiometric evaluation included mean and variance analysis of individual detectors and principal components analysis. Results indicate that detector bias for all bands is very close or within tolerance. Bright spots were observed in the thermal IR band on an 18 line by 128 pixel grid. No explanation for this was pursued. The general overall quality of the TM was judged to be very high.
Simulation experiments comparing the relative importance of an idealized LIDAR wind profiling system with idealized temperature and pressure sounding systems on 12 h forecasts are studied for three "nature' fields representing the true evolving atmospheric states. The three fields are obtained respectively from: (1) a long integration of the GLAS 4th Order Model (4 deg lat x 5 deg lon x 9 levels), (2) a continuous sequence of NMC operational analysis and, (3) a long integration of the ECMWF high resolution (1.875 deg 1 lat x 1.875 deg lon x 15 layers) operational forecast model. These fields are interpolated to the grid of the GLAS model and used for simulating the observed global analysed fields of winds, temperature, moisture and surface pressure. The same interpolated fields are also used for verification of forecast impact. The effects of clouds, aerosol concentrations, and instrument accuracies on the simulated observation will be discussed.