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.
The following records were processed for inclusion into the NA&SD database:Basic Lidar Principles and Techniques; Applications I: Radiative Forcing and Climate; Boundary Layer and Pollution; Scattering of laser light:more than just smoke and mirrors; Applications II: Water Vapor and Atmospheric Dynamics; Regional and Global Data Available Now And In The Future.
A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.
This paper provides a brief history of the F-16XL-1 aircraft, its role in the High Speed Research program and how it was morphed into the Cranked Arrow Wing Aerodynamics Project. Various flight, wind-tunnel and Computational Fluid Dynamics data sets were generated as part of the project. These unique and open flight datasets for surface pressures, boundary-layer profiles and skin-friction distributions, along with surface flow data, are described and sample data comparisons given. This is followed by a description of how the project became internationalized to be known as Cranked Arrow Wing Aerodynamics Project International and is concluded by an introduction to the results of a four year computational predictive study of data collected at flight conditions by participating researchers.
Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.
Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are characterized by high melting point, good strength, and reasonable oxidation resistance. These materials are of interest for use as sharp leading edges for hypersonic vehicles among other applications. Progress in computational modeling of UHTCs has been limited in part due to the absence of suitable interatomic potentials. We present a Tersoff style parameterization of such potentials for ZrB2 and HfB2 appropriate for atomistic simulations. Parameters are fit to data generated from ab initio calculations. The accuracy of the potentials is assessed against further ab initio data. As a first non ]trivial application, molecular dynamics simulations are performed to evaluate the thermal conductivity of single crystals and the thermal resistance of high symmetry grain boundaries.
This report is the documentation of the work performed under the Hypersonic Project of the NASA's Fundamental Aeronautics Program. It was funded through Task Number NNC10E444T under GESS-2 Contract NNC06BA07B. The objective of the task is to develop advanced computational tools for the simulation of multi-stage turbomachinery in support of aeropropulsion. This includes work elements in extending the TURBO code and validating the multi-stage URANS (Unsteady Reynolds Averaged Navier Stokes) simulation results with the experimental data. The unsteady CFD (Computation Fluid Dynamics) calculations were performed in full wheel mode with and without screen generated total pressure distortion at the computational inflow boundary, as well as in single passage phase lag mode for uniform inflow. The experimental data were provided by NASA from the single stage RTA (Revolutionary Turbine Accelerator) fan test program.Significant non-uniform flow condition at the fan-face of the aeropropulsion system is frequentlyencountered in many of the advanced aerospace vehicles. These propulsion systems can be eithera podded or an embedded design employed in HWB (Hybrid Wing Body) airframe concept. It isalso a topic of interest in military applications, in which advanced air vehicles have already deployedsome form of embedded propulsion systems in their design because of the requirementsof compact and low observable inlets. Even in the conventional airframe/engine design, the fancould operate under such condition when the air vehicle is undergoing rapid maneuvering action.It is believed that a better understanding of the fan’s aerodynamic and aeromechanical responseto this type of operating condition or off design operation would be beneficial to designing distortiontolerant blades for improved engine operability.The objective for this research is to assess the capability of turbomachinery code as an analysistool in understanding the effects and evaluating the impact of flow distortion on the aerodynamicand aeromechanical performance of the fan in advanced propulsion systems. Results from thetesting of an advanced fan stage released by NASA are available and will be used here for CFDcode validation. The experiment was performed at NASA’s high speed compressor facility aspart of the RTA (Revolutionary Turbine Accelerator) demonstration project, a joint effort ofNASA Glenn Research Center and GE Aircraft Engines in developing an advanced Mach 4TBCC (Turbine Based Combined Cycle) turbofan/ramjet engine for access to space. Part of thetest was to assess the aerodynamic performance and operability of the fan stage under nonuniforminflow condition. Various flow distortion patterns were created at the fan-face by manipulatingsets of screens placed upstream of the wind tunnel. Measurements at the fan-face willprovide the necessary distortion flow information as the inflow boundary condition for the CFDin a full wheel simulation. Therefore the purpose of this work is to demonstrate the NASA supportedmulti-stage turbomachinery code, TURBO [1-5], in the aerodynamic performance analysisof a modern fan design operating under off design condition, and in particular to validate theCFD results with the RTA fan test data.A brief description of the RTA fan rig configuration is given in the next section, explaining onhow flow distortion were measured in the test and constructed for the CFD at the fan-face. It isfollowed by a section summarizing previous CFD work performed at NASA relevant to the currentfan configuration. A short description of the TURBO code is given next, followed by detailsin the computational model of the fan rig, the required computing resources, and the numericalprocedure for the simulations. The CFD results are presented in the discussion section and finallyconcluding remarks are summarized.
Adaptive Mesh Refinement (AMR) promises a much more computationally efficient means to obtain a discrete approximation to a continuous boundary value problem of a specified accuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW (a computational fluid dynamics (CFD) code) is utilized to provide estimates of the exact analytical solutions to problems of interest to turbulence modeling. Predictions of surface pressure and skin friction, essentially the state of stress at the surface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on the other hand, show marked differences in flows with shocks. The AMR method, as implemented in OVERFLOW 2.2k, appears to provide the ability to produce arbitrarily accurate solutions at a predictable cost much smaller than classic uniform mesh refinement.
The lunar environment heliospheric X-ray imager (LEXI) and solar wind–magnetosphere–ionosphere link explorer (SMILE) will observe the magnetopause motion in soft X-rays to understand dayside reconnection modes as a function of solar wind conditions after their respective launches in the near future. To support their successful science mission, we investigate the relationship between the magnetopause position and the dayside reconnection rate by utilizing super dual auroral radar network (SuperDARN) observations and widely used empirical models of magnetopause position (Shue et al. in J Geophys Res 103:17691–17700. https://doi.org/10.1029/98JA01103, 1998 and Lin et al. in J Geophys Res 115:A04207. https://doi.org/10.1029/2009JA014235, 2010). We select three cases when the interplanetary magnetic field rotates during periods of quasi-steady solar wind dynamic pressure. We first estimate the dayside reconnection rate by calculating the electric field along the open–closed magnetic field boundary (OCB) in the OCB moving reference frame. Then, we estimate the magnetopause position near the local noon by inputting NASA OMNI solar wind data into the empirical magnetopause models. The reconnection rate shows anti-correlation with the magnetopause position that it generally increases as the magnetopause located closer to Earth and vice versa. Our result also confirms that the reconnection rate increases as the empirical coupling efficiency between solar wind and the magnetosphere increases.
Assessing energy resources under future scenarios requires high-resolution meteorological information that is physically consistent and suitable for regional-scale analysis. While Earth system model (ESM) projections provide valuable large-scale information, their coarse resolution and systematic biases limit direct applicability for energy system modeling and planning. In this study, we develop a high-resolution dynamical downscaling framework based on the Weather Research and Forecasting (WRF) model to translate global-scale ESM data into energy-relevant regional projections over the contiguous United States (CONUS). The framework identifies an optimized WRF configuration through numerical experiments and evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture (SM) and soil temperature (ST) bias correction implemented as an integral part of the bias-corrected ESM forcing to improve land-atmosphere coupling prior to WRF dynamical downscaling. Using an optimized WRF configuration at 4-km resolution, we show that raw ESM forcing introduces systematic dry and cold soil biases that propagate into pronounced warm biases in near-surface air temperature and positive biases in solar irradiance, particularly during summer. Applying bias-corrected atmospheric forcing together with bias-corrected SM and ST substantially reduces these downstream biases and improves the surface energy balance and near-surface atmospheric fields. These results demonstrate that bias-aware treatment of initial conditions is critical for producing high-resolution downscaled projections suitable for energy system modeling and planning applications.
Dielectric particles in weakly conducting fluids rotate spontaneously when subject to strong electric fields. Such Quincke rotation near a plane electrode leads to particle translation that enables physical models of active matter. In this Letter, we show that Quincke rollers can also exhibit oscillatory dynamics, whereby particles move back and forth about a fixed location. We explain how oscillations arise for micron-scale particles commensurate with the thickness of a field-induced boundary layer in the nonpolar electrolyte. This work enables the design of colloidal oscillators.
Lead halide perovskites are strong candidates for high-performance low-cost photovoltaics, light emission, and detection applications. A hot-phonon bottleneck effect significantly extends the cooling time of hot charge carriers, which thermalize through carrier–optic phonon scattering, followed by optic phonon decay to acoustic phonons and finally thermal conduction. To understand these processes, we adjust the lattice dynamics independently of electronics by changing isotopes. We show that doubling the mass of hydrogen in methylammonium lead iodide by replacing protons with deuterons causes a large 20 to 50% softening of the longitudinal acoustic phonons near zone boundaries, reduces thermal conductivity by ~50%, and slows carrier relaxation kinetics. Phonon softening is attributed to anticrossing with the slowed libration modes of the deuterated molecules and the reduced thermal conductivity to lowered phonon velocities. Our results reveal how tuning the organic molecule dynamics enables control of phonons important to thermal conductivity and the hot-phonon bottleneck.
Boundary element algorithms for the solution of steady-state and transient heat conduction are presented. The algorithms are designed for efficient coupling with computational fluid dynamic discretizations and feature piecewise linear elements with offset nodal points. The steady-state algorithm employs the fundamental solution approach; the integration kernels are computed analytically based on linear shape functions, linear elements, and variably offset nodal points. The analytic expressions for both singular and nonsingular integrands are presented. The transient algorithm employs the transient fundamental solution; the temporal integration is performed analytically and the nonsingular spatial integration is performed numerically using Gaussian quadrature. A series solution to the integration is derived for the instance of a singular integrand. The boundary-only character of the algorithm is maintained by integrating the influence coefficients from initial time. Numerical results are compared to analytical solutions to verify the current boundary element algorithms. The steady-state and transient algorithms are numerically shown to be second-order accurate in space and time, respectively.
The relationship between blocking and multiple equilibria was investigated in a two layer model of the Northern Hemisphere January heating and boundary conditions are imposed. The truncation varies over global wavenumbers 4-8 and the coefficients of horizontal diffusion of heat and potential vorticity are varied about realistic values. It is shown that while severely truncated models can exhibit multiple quasi equilibria under realistic domain, forcing, and boundary conditions, higher resolution is required to allow realistic transitions to be mediated by strong baroclinic wave activity. It is indicated that for extended range forecasting blocking transitions are not likely to be forecast well unless the associated synoptic events can be forecasted statistically.
In this paper we describe an approach for dealing with arbitrary complex, two dimensional geometries, the so-called cartesian boundary method. Conceptually, the cartesian boundary method is quite simple. Solid bodies blank out areas of a background, cartesian mesh, and the resultant cut cells are singled out for special attention. However, there are several obstacles that must be overcome in order to achieve a practical scheme. We present a general strategy that overcomes these obstacles, together with some details of our successful conversion of an adaptive mesh algorithm from a body-fitted code to a cartesian boundary code.
Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.
As a portion of the Benchmark Models Program at NASA Langley, a simple generic model was developed for active controls research and was called BACT for Benchmark Active Controls Technology model. This model was based on the previously-tested Benchmark Models rectangular wing with the NACA 0012 airfoil section that was mounted on the Pitch and Plunge Apparatus (PAPA) for flutter testing. The BACT model had an upper surface spoiler, a lower surface spoiler, and a trailing edge control surface for use in flutter suppression and dynamic response excitation. Previous experience with flutter suppression indicated a need for measured control surface aerodynamics for accurate control law design. Three different types of flutter instability boundaries had also been determined for the NACA 0012/PAPA model, a classical flutter boundary, a transonic stall flutter boundary at angle of attack, and a plunge instability near M = 0.9. Therefore an extensive set of steady and control surface oscillation data was generated spanning the range of the three types of instabilities. This information was subsequently used to design control laws to suppress each flutter instability. There have been three tests of the BACT model. The objective of the first test, TDT Test 485, was to generate a data set of steady and unsteady control surface effectiveness data, and to determine the open loop dynamic characteristics of the control systems including the actuators. Unsteady pressures, loads, and transfer functions were measured. The other two tests, TDT Test 502 and TDT Test 5 18, were primarily oriented towards active controls research, but some data supplementary to the first test were obtained. Dynamic response of the flexible system to control surface excitation and open loop flutter characteristics were determined during Test 502. Loads were not measured during the last two tests. During these tests, a database of over 3000 data sets was obtained. A reasonably extensive subset of the data sets from the first two tests have been chosen for Test Cases for computational comparisons concentrating on static conditions and cases with harmonically oscillating control surfaces. Several flutter Test Cases from both tests have also been included. Some aerodynamic comparisons with the BACT data have been made using computational fluid dynamics codes at the Navier-Stokes level (and in the accompanying chapter SC). Some mechanical and active control studies have been presented. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. Cases for static angles of attack, static trailing-edge and upper-surface spoiler deflections are included for a range of conditions near those for the oscillation cases. Cases for trailing-edge control and upper-surface spoiler oscillations for a range of Mach numbers, angle of attack, and static control deflections are included. Cases for all three types of flutter instability are selected. In addition some cases are included for dynamic response measurements during forced oscillations of the controls on the flexible mount. An overview of the model and tests is given, and the standard formulary for these data is listed. Some sample data and sample results of calculations are presented. Only the static pressures and the first harmonic real and imaginary parts of the pressures are included in the data for the Test Cases, but digitized time histories have been archived. The data for the Test Cases are also available as separate electronic files.
The lack of accurate submodels for in-cylinder radiation and heat transfer has been identified as a key shortcoming in developing truly predictive, physics-based computational fluid dynamics (CFD) models that can be used to develop and design combustion systems for advanced high-efficiency, low-emissions engines. Measurements of wall boundary layers in engines show discrepancies of up to 100% with respect to standard CFD boundary layer models. And analysis of in-cylinder radiation based on the most recent spectral property databases and high-fidelity radiative transfer equation solvers has shown that at operating pressures and exhaust-gas recirculation levels typical of modern compression-ignition truck engines, molecular gas radiation (mainly carbon dioxide and water vapor) can be more important than soot particle radiation, and that a significant fraction of the emitted radiation (50% or more) can be reabsorbed before reaching the walls, thereby acting to redistribute energy in the combustion chamber. In this project, an integrated experimental and simulation approach has been used to generate new physical insight into in-cylinder heat-transfer processes in engines, and to develop computationally practicable CFD-based models that accurately represent the physics. New experimental data have been taken in an optically accessible single-cylinder research engine, and in a multi-cylinder production engine. Those data have provided new insight into the radiative environment in engines. New CFD models have been developed and validated for in-cylinder heat transfer, focusing mainly on radiative heat transfer and coupling between turbulent boundary layer wall heat transfer and radiative heat transfer. The CFD models have been developed and implemented in an open-source software framework that can be shared with other researchers. The models are compatible both with unsteady Reynolds-averaged Navier Stokes (URANS) formulations and with large-eddy simulations (LES) of in-cylinder processes in engines.
Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.