Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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.

At least 1,333 records · Page 74

Computer program for calculating laminar and turbulent boundary layer development in compressible flow

A computer program is described which performs a numerical integration of the equations of motion for a compressible two-dimensional boundary layer. Boundary layer calculations may be carried out for both laminar and turbulent flow for arbitrary Reynolds number and free stream Mach number distribution on planar or axisymmetric bodies with wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and a rough or a smooth wall. A variety of options are available as initial conditions. The program can generate laminar initial conditions such as Falkner-Skan similarity solutions (so that initial wedge flows can be simulated including Blasius or stagnation point flow) or approximate equilibrium turbulent profiles. Alternatively, initial profile input data can be utilized.

Herring, H. J.↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

Assimilation of Polar Orbiting Space-Based Data and Its Impact on Weather Prediction

This portion of the Data Assimilation and Numerical Weather Prediction Benefits presentation will describe: the need for remotely sensed data, the types of observations that have been made and will be available, how satellite data have been and are being assimilated, the impact of the data on weather analysis and prediction in both the U.S and abroad and the potential impact of new observing systems on NWP. Specific observing systems that will be discussed include: satellite temperature and moisture soundings/radiances (TOVS, ATOVS, AIRS/AMSU), satellite surface winds from both passive and active microwave atmospheric motion winds (from MODIS or MISR) precipitation measurements clouds and land surface measurements (e.g. from MODIS) lidar wind profiles.

Atlas, Robert M.↗

Compact, Engineered, 2-Micron Coherent Doppler Wind Lidar Prototype: A New NASA Instrument Incubator Program Project

A new project, selected in 2005 by NASA s Science Mission Directorate (SMD) under the Instrument Incubator Program (IIP), will be described. The 3-year effort is intended to design, fabricate, and demonstrate a packaged, rugged, compact, space-qualifiable coherent Doppler wind lidar (DWL) transceiver capable of future validation in an aircraft and/or Unmanned Aerial Vehicle (UAV). The packaged DWL will utilize the numerous advances in pulsed, solid-state, 2-micron laser technology at NASA s Langley Research Center (LaRC) in such areas as crystal composition, architecture, efficiency, cooling techniques, pulse energy, and beam quality. The extensive experience of Raytheon Space and Airborne Systems (RSAS) in coherent lidar systems, in spacebased sensors, and in packaging rugged lidar systems will be applied to this project. The packaged transceiver will be as close to an envisioned space-based DWL system as the resources and technology readiness allow. We will attempt to facilitate a future upgrade to a coherent lidar system capable of simultaneous wind and CO2 concentration profile measurements. Since aerosol and dust concentration is also available from the lidar signal, the potential for a triple measurement lidar system is attractive for both Earth and Mars remote sensing. A key follow on step after the IIP will be to add a telescope, scanner, and software for aircraft validation. This IIP should also put us in a position to begin a parallel formulation study in the 2006-2007 timeframe for a space-based DWL demonstration mission early next decade.

Kavaya, Michael J.↗

Comparison of Mesospheric Winds From a High-Altitude Meteorological Analysis System and Meteor Radar Observations During the Boreal Winters of 2009-2010 and 2012-2013

We present a study of horizontal winds in the mesosphere and lower thermosphere (MLT) during the boreal winters of 2009-2010 and 2012-2013 produced with a new high-altitude numerical weather prediction (NWP) system. This system is based on a modified version of the Navy Global Environmental Model (NAVGEM) with an extended vertical domain up to approximately 116 km altitude coupled with a hybrid four-dimensional variational (4DVAR) data assimilation system that assimilates both standard operational meteorological observations in the troposphere and satellite-based observations of temperature, ozone and water vapor in the stratosphere and mesosphere. NAVGEM-based MLT analyzed winds are validated using independent meteor radar wind observations from nine different sites ranging from 69 deg N-67 deg S latitude. Time-averaged NAVGEM zonal and meridional wind profiles between 75 and 95 km altitude show good qualitative and quantitative agreement with corresponding meteor radar wind profiles. Wavelet analysis finds that the 3-hourly NAVGEM and 1-hourly radar winds both exhibit semi-diurnal, diurnal, and quasi-diurnal variations whose vertical profiles of amplitude and phase are also in good agreement. Wavelet analysis also reveals common time-frequency behavior in both NAVGEM and radar winds throughout the Northern extra tropics around the times of major stratospheric sudden warmings (SSWs) in January 2010 and January 2013, with a reduction in semi-diurnal amplitudes beginning around the time of a mesospheric wind reversal at 60 deg N that precedes the SSW, followed by an amplification of semi-diurnal amplitudes that peaks 10-14 days following the onset of the mesospheric wind reversal. The initial results presented in this study demonstrate that the wind analyses produced by the high altitude NAVGEM system accurately capture key features in the observed MLT winds during these two boreal winter periods.

McCormack, J.↗

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Baurle, R. A.↗

BENEFIT AND CRITICAL FACTORS FOR THE PERFORMANCE OF THE BOUNDARY LAYER INGESTING PROPULSION

The benefit of the boundary layer ingestion(BLI)is described in the perspective of the propulsion and engine development. A power saving map of the BLI engines is derived based on the correlation of the wake velocity ratio of the ingested boundary layer profile and the propulsive efficiency. The ratio of the mass flow rate between BLI and non-BLI propulsors is introduced to quantify the power saving of the BLI engine relative to a clean inlet flow engine which generates same amount of thrust. The wake recovery factor from the jet flow out of the BLI engine is employed to find an adequate sizing of the BLI engine for the given design requirement. The effects of the fan pressure ratio on the power saving are also investigated to explore the feasible range of the BLI engine design. The derived correlation is validated with CFD analyses. A numerical experiment is carried out by varying the wake velocity ratio through different BLI engines sized with respect to an influencing body. Consequently, the propulsor efficiency is quantified and presented by the saving in the actual shaft power. The efficiency penalty, pressure ratio of the BLI fan stage are correlated with the power saving and the correlation is validated through BLI2DTF and R4 fan stage CFD analyses based on rig test data.

B J Lee↗

A Deep Learning Approach to Fast Radiative Transfer

Due to the sheer volume of data, leveraging satellite instrument observations effectively in a data assimilation context for numerical weather prediction or for remote sensing requires a radiative transfer model as an observation operator that is both fast and accurate at the same time. Physics-based line-by-line radiative transfer (RT) models fulfil the requirement for accuracy, but are too slow and too costly in computational terms for operational applications. Therefore, fast methods were developed to be able to perform fast RT calculations using techniques such as spectral sampling or pre-computed look-up tables. The operational fast models currently calculate the absorption and scattering coefficients from the pre-computed regression coefficients and atmospheric state and cloud profiles. As a novel solution to this problem, this work investigates a deep learning approach to replace the regression coefficients in the fast RT models. A selection of hidden-layer neural network configurations is trained against atmospheric transmittance profile data computed by an accurate line-by-line model and their performance is evaluated and their advantages and disadvantages are discussed.

Machine learning↗

Senteniel-6 Radio Occultation Product Released by NASA GES DISC to Supplement Satellite Remote Sensing Datasets for PBL Study

The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) curates hyperspectral atmospheric sounder remote sensing and numerical model reanalysis datasets which have been utilized in the Planetary Boundary Layer (PBL) research and applications. The hyperspectral sounder remote-sensing datasets include the Atmospheric Infrared Sounder (AIRS) on the Aqua satellite to the Cross-track Infrared Sounder (CrIS) on Suomi--National Polar- orbiting Partnership (NPP) and National Oceanic and Atmospheric Administration -20 (N NOAA-20)/ Joint Polar-orbiting Satellite System -1 (JPSS-1). The Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2) global reanalysis product provides a data record commencing in 1980. The sounder remote sensing and reanalysis datasets include temperature, water vapor, and trace gas profile down to the PBL, and also have a derived PBL height as well. A nearly 10-year (June 2006 to December 2015) seasonal and annual PBL height climatology dataset from COSMIC Global Navigation Satellite System (GNSS) radio occultation (RO) measurement is also available from the GES DISC. In collaboration with Sentinel-6 Project, the GES DISC is implementing curation activities for GNSS RO products from the Sentinel-6A/Sentinel-6 Michael Freilich satellite launched on November 21, 2020. Sentinel-6A RO products provide refractivity, temperature, and humidity profile with finer vertical resolution, leveraging PBL research and application as a supplement to the hyperspectral sounder remote sensing and reanalysis products. The public release of Sentinel- 6A RO products is scheduled for mid-October of 2021. In this presentation, we will introduce all Senitnel-6A products and services, and demonstrate use cases studying the PBL by combining these products with other GES DISC archived data products.

Feng Ding↗

Pressure- and buoyancy-driven thermal convection in a rectangular enclosure

Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.

Spradley, L. W.↗

Calculation of shock-separated turbulent boundary layers

Numerical solutions of the complete, time-averaged conservation equations using several eddy-viscosity models for the Reynolds shear stress to close the equations are compared with experimental measurements in a compressible, turbulent separated flow. An efficient time-splitting, explicit difference scheme was used to solve the two-dimensional conservation equations. The experiment used for comparison was a turbulent boundary layer that was separated by an incident shock wave in a Mach 2.93 flow with a unit Reynolds number of 5.7 x 10 to the seventh power m. Comparisons of predicted and experimental values of surface pressure, shear stress along the wall, and velocity profiles are shown. One of the tested eddy-viscosity models which allows the shear stress to be out of equilibrium with the mean flow produces substantially better agreement with the experimental measurements than the simpler models. A tool is thereby provided for inferring additional information about the flow, such as static pressures in the stream, which might not be directly obtainable from experiments.

Baldwin, B. S.↗

Spectral analysis of GEOS-3 altimeter data and frequency domain collocation

The mathematical background in spectral analysis as applied to geodetic applications is summarized. The resolution (cut-off frequency) of the GEOS 3 altimeter data is examined by determining the shortest wavelength (corresponding to the cut-off frequency) recoverable. The data from some 18 profiles are used. The total power (variance) in the sea surface topography with respect to the reference ellipsoid as well as with respect to the GEM-9 surface is computed. A fast inversion algorithm for matrices of simple and block Toeplitz matrices and its application to least squares collocation is explained. This algorithm yields a considerable gain in computer time and storage in comparison with conventional least squares collocation. Frequency domain least squares collocation techniques are also introduced and applied to estimating gravity anomalies from GEOS 3 altimeter data. These techniques substantially reduce the computer time and requirements in storage associated with the conventional least squares collocation. Numerical examples given demonstrate the efficiency and speed of these techniques.

Eren, K.↗

A determination of the external forces required to move the benchmark active controls testing model in pure plunge and pure pitch

In view of the strong need for a well-documented set of experimental data which is suitable for the validation and/or calibration of modern Computational Fluid Dynamics codes, the Benchmark Models Program was initiated by the Structural Dynamics Division of the NASA Langley Research Center. One of the models in the program, the Benchmark Active Controls Testing Model, consists of a rigid wing of rectangular planform with a NACA 0012 profile and three control surfaces (a trailing-edge control surface, a lower-surface spoiler, and an upper-surface spoiler). The model is affixed to a flexible mount system which allows only plunging and/or pitching motion. An approximate analytical determination of the forces required to move this model, with its control surfaces fixed, in pure plunge and pure pitch at a number of test conditions is included. This provides a good indication of the type of actuator system required to generate the aerodynamic data resulting from pure plunging and pure pitching motion, in which much interest was expressed. The analysis makes use of previously obtained numerical results.

Dcruz, Jonathan↗

Regional Assimilation of NASA Atmospheric Infrared Sounder (AIRS) Data

The NASA Short-term Prediction Research and Transition (SPORT) Center seeks to accelerate the infusion of NASA Earth Science Enterprise (ESE) observations, data assimilation and modeling research into NW S forecast operations and decision-making. The Atmospheric Infrared Sounder (AIRS), is expected to advance climate research and weather prediction into the 21 st century. It is one of six instruments onboard Aqua, a satellite that is part of NASA s Earth Observing System. AIRS, along with two partner microwave sounding instruments, represents the most advanced atmospheric sounding system ever deployed in space. The system is capable of measuring the atmospheric temperature in the troposphere with radiosonde accuracies of 1 K over 1 km-thick layers under both clear and cloudy conditions, while the accuracy of the derived moisture profiles will exceed that obtained by radiosondes. It is imperative that the scientific community is prepared to take full advantage of next-generation satellite data that will become available within the next decade. The purpose of this paper is to describe a procedure designed to optimally assimilate AIRS data at high spatial resolution over both land and ocean. The assimilation system used in this study is the Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory used extensively around the globe. Results will focus on quality control issues associated with AIRS, optimal assimilation strategies, and the impact of the AIRS data on subsequent numerical forecasts at 12 km produced by the next generation Weather Research and Forecast (WRF) model.

Chou, Shih-Hung↗

Shock waves data for minerals

Shock compression of the materials of planetary interiors yields data which upon comparison with density-pressure and density-sound velocity profiles constrain internal composition and temperature. Other important applications of shock wave data and related properties are found in the impact mechanics of terrestrial planets and solid satellites. Shock wave equation of state, shock-induced dynamic yielding and phase transitions, and shock temperature are discussed. In regions where a substantial phase change in the material does not occur, the relationship between the particle velocity, U(sub p), and the shock velocity, U(sub s), is given by U(sub s) = C(sub 0) + S U(sub p), where C(sub 0) is the shock velocity at infinitesimally small particle velocity, or the ambient pressure bulk sound velocity. Numerical values for the shock wave equation of state for minerals and related materials of the solar system are provided.

Ahrens, Thomas J.↗

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗