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.
This excel file provides direct use geothermal data for various wells located in Oregon. Data includes well name, location, heating application, temperature and flow of fluid, capacity, energy use, and load factor. The data is from the State of Oregon Department of Geology and Mineral Industries (linked below), where additional data downloads and metadata may be found.
Supersonic wind tunnel testing of 0.813 m diameter Disk-Gap-Band parachutes is being conducted in the NASA Glenn Research Center (GRC) 10' x 10' wind-tunnel. The tests are conducted in support of the Mars Science Laboratory Parachute Decelerator System development and qualification. Four percent of full-scale parachutes were constructed similarly to the flight-article in material and construction techniques. The parachutes are attached to a 4% scale MSL entry-vehicle to simulate the free-flight configuration. The parachutes are tested from Mach 2 to 2.5 over a Reynolds number (Re) range of 1 to 3 x 10(exp 6), representative of the MSL deployment envelope. Constrained and unconstrained test configurations are investigated to quantify the effects of parachute trim, suspension line interaction, and alignment with the capsule wake. The parachute is constrained horizontally through the vent region, to measure canopy breathing and wake interaction for fixed trim angles of 0 and 10 degrees from the velocity vector. In the unconstrained configuration the parachute is permitted to trim and cone, similar to the free-flight varying its alignment relative to the entry-vehicle wake. Test diagnostics were chosen to quantify parachute performance and to provide insight into the flow field structure. An in-line load cell provided measurement of unsteady and mean drag as a function of Mach and Re. High-speed shadowgraph video of the upstream parachute flow field was used to capture bow-shock motion and stand of distance. Particle image velocimetry of the upstream parachute flow field provides spatially and temporally resolved measurement velocity and turbulent statistics. Multiple high speed video views of targets placed in the interior of the canopy enable photo-grammetric measurement of the fabric motion in time and space from reflective. High speed video is also used to document the supersonic inflation and measure trim angle, projected area, and frequency of area oscillations.
In this article, we describe off-design models and control strategies for a Pumped Thermal Energy Storage (PTES) system that uses liquid thermal energy storage: specifically molten salt for hot storage and methanol for cold storage. Off-design conditions arise when load-following, or due to variations in storage tank temperatures or ambient temperatures. We propose a control strategy that uses inventory control to manage the mass flow rate in the thermodynamic cycles, which facilitates load following. We also propose a control strategy for the storage fluid mass flow rates, which are varied to ensure the molten salt is maintained at its design temperature. This maximizes efficiency and minimizes problems with salt freezing or degradation. The cold storage fluid mass flow rate is varied so that the cold tanks have the same state-of-charge as the hot tanks. This leads to variations in cold fluid temperature, but these variations are shown to be acceptably small (e.g. 7.5% increase), and this control method is shown to be simpler and more efficient than an alternative strategy where tanks become unbalanced. The ambient temperature and storage tank temperatures are moved ±50 °C from the design values and the impact on power, duration, and tank temperatures is quantified. Results demonstrate that the proposed control strategy is stable and self-correcting - that is, storage temperatures converge on stable values after two-to-three charge-discharge cycles. When inputs return to design values, the system returns to its design point after two charge-discharge cycles. We also demonstrate that inventory control enables delivery of the target power output even when off-design conditions exist that would normally reduce the power output.
Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.
Ground crew veterans at Kennedy Space Center still talk about what they call "the summer of hydrogen"-the long, frustrating months in 1990 when the shuttle fleet was grounded by an elusive hydrogen leak that foiled our efforts to fill the orbiter's external fuel tank. Columbia (STS-35) was on Launch Pad A for a scheduled May 30 launch when we discovered the hydrogen leak during - tanking. The external fuel tank is loaded through the orbiter. Liquid hydrogen flows through a 17-inch umbilical between the orbiter and the tank. During fueling, we purge the aft fuselage with gaseous nitrogen to reduce the risk of fire, and we have a leak-detection system in the mobile launch platform, which samples (via tygon tubing) the atmosphere in and around the vehicle, drawing it down to a mass spectrometer that analyzes its composition. When we progressed to the stage of tanking where liquid hydrogen flows through the vehicle, the concentration of hydrogen approached four percent-the limit above which it would be dangerously flammable. We had a leak. We did everything we could think of to find it, and the contractor who supplied the flight hardware was there every day, working alongside us. We did tanking tests, which involved instrumenting the suspected leak sources, and cryo-loaded the external tank to try to isolate precisely where the leak originated. We switched out umbilicals; we replaced the seals between the umbilical and the orbiter. We inspected the seals microscopically and found no flaws. We replaced the recirculation pumps, and we found and replaced a damaged teflon seal in a main propulsion system detent cover, which holds the prevalve-the main valve supplying hydrogen to Space Shuttle Main Engine 3 -in the open position. The seal passed leak tests at ambient temperature but leaked when cryogenic temperatures were applied. We added new leak sensors-up to twenty at a time and tried to be methodical in our placements to narrow down the possible sources of the problem. We even switched orbiters, sending Columbia back to the Vehicle Assembly Building and bringing out Atlantis, scheduled to fly as STS-38. Two shuttles on their mobile launchers passing in the night was a majestic sight, but not one you want to see if you're trying to get an orbiter launched. None of this told us where the leak was, or if we were dealing with more than one leak source.
An assessment of a Hybrid RANS/LES (HRLES) approach for C L,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both free air and wind tunnel configurations of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES is shown to address some of the known shortcomings with RANS and prevent inboard and outboard flow separation particularly in the region of C L,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly underperformed RANS. For the in-tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.
A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.
The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.
Detailed surface heat transfer and pressure distributions have been obtained in three-dimensional shock-wave boundary-layer interactions flow regions. The data described were obtained on fundamental shapes: planar wings with trailing edge flaps or spoilers and planar or cylindrical center bodies, representative of the aft portion of hypersonic aircraft. An overview of the work is presented; details of the projects are available in many reports in the open literature. Analytic, empiric methods are advanced for predicting the extent of separation and the increased heat transfer and pressure loads in three-dimensional separated flow regions.
A class of power converters is described for supplying direct current at one voltage from a source at another voltage. It includes a simple passive circuit arrangement of solid-state switches, inductors, and capacitors by which the output voltage of the converter tends to remain constant in spite of changes in load. The switches are sensitive to the current flowing in the circuit and are employed to permit the charging of capacitance devices in accordance with the load requirements. Because solid-state switches (such as SCR's) may be used with relatively high voltage and because of the inherent efficiency of the invention that permits relatively high switching frequencies, power supplies built in accordance with the invention, together with their associated cabling, can be substantially lighter in weight for a given output power level and efficiency of operation than systems of the prior art.
Tests were performed to evaluate three different methods of utilizing air to provide thermal management control for compliant journal foil air bearings. The effectiveness of the methods was based on bearing bulk temperature and axial thermal gradient reductions during air delivery. The first method utilized direct impingement of air on the inner surface of a hollow test journal during operation. The second, less indirect method achieved heat removal by blowing air inside the test journal to simulate air flowing axially through a hollow, rotating shaft. The third method emulated the most common approach to removing heat by forcing air axially through the bearing s support structure. Internal bearing temperatures were measured with three, type K thermocouples embedded in the bearing that measured general internal temperatures and axial thermal gradients. Testing was performed in a 1 atm, 260 C ambient environment with the bearing operating at 60 krpm and supporting a load of 222 N. Air volumetric flows of 0.06, 0.11, and 0.17 cubic meters per minute at approximately 150 to 200 C were used. The tests indicate that all three methods provide thermal management but at different levels of effectiveness. Axial cooling of the bearing support structure had a greater effect on bulk temperature for each air flow and demonstrated that the thermal gradients could be influenced by the directionality of the air flow. Direct air impingement on the journal's inside surface provided uniform reductions in both bulk temperature and thermal gradients. Similar to the direct method, indirect journal cooling had a uniform cooling effect on both bulk temperatures and thermal gradients but was the least effective of the three methods.
This paper presents results of wind-tunnel and acoustic tests to investigate buffet loads on Shuttle Thermal-Protection-System (TPS) tiles, It also describes the application of these results to the prediction of tile buffet loads for the first shuttle flight into orbit (STS-I). The wind-tunnel tests of tiles were conducted at transonic and supersonic Math numbers simulating flow regions on the Orbiter where shock waves and boundary-layer separations occur. The acoustic tests were conducted in a progressive wave tube at an overall sound pressure level (OASPL) approximately equal to the maximum OASPL measured during the wind-tunnel tests in a region of flow separation. The STS-I buffet load predictions yielded peak tile stresses due to buffeting, that were as much as 20 percent of the total stress for the design-load case when a shock wave was on a tile.
A user's manual is presented for a computer program in which a vortex-lattice lifting-surface method is used to model the wing and multiple flaps. The engine wake model consists of a series of closely spaced vortex rings with rectangular cross sections. The jet wake is positioned such that the lower boundary of the jet is tangent to the wing and flap upper surfaces. The two potential flow models are used to calculate the wing-flap loading distribution including the influence of the wakes from up to two engines on the semispan. The method is limited to the condition where the flow and geometry of the configurations are symmetric about the vertical plane containing the wing root chord. The results include total configuration forces and moments, individual lifting-surface load distributions, pressure distributions, flap hinge moments, and flow field calculation at arbitrary field points. The use of the program, preparation of input, the output, program listing, and sample cases are described.
Experimental modelling of the interaction between a jet and an aircraft wing or fuselage in VTOL aircraft was undertaken using a cold jet exiting perpendicular to a flat plate in a uniform cross-flow. Effects of jet decay rate and jet-to-cross-flow velocity ratio, R, on the induced load distribution were investigated. Jet decay rate was increased by using cylindrical centerbodies submerged in the jet nozzle, which caused nonuniform initial jet velocity profiles. Quicker jet decay rate, corresponding to the presence of a centerbody, resulted in as much as 50% reduction in the induced pressure loads on the plate. This has implications in interpretation of results from earlier VTOL model studies of jet induced loads, where the jets have often had relatively slow decay rates due to uniform initial velocity profiles
Results are presented from an investigation to measure the flow field generated by a forward-swept propfan operating in flutter at a low forward velocity. For comparison to the flutter condition, flow field data are also presented for a slightly reduced rotational speed just below flutter. The forward-swept propfan was tested as the front rotor of a counterrotating pusher propeller. A laser Doppler velocimeter (LDV) was used to measure the velocity field in planes normal to the model centerline downstream of the rotor and in planes of constant radius within the blade passages at each operating condition. The mean time-averaged flow about the blades did not change drastically as the propfan rotational speed was increased from the stable operating point to the flutter condition. No regions of flow separation could be identified in the data plots of the mean intrablade flow field. The relative flow about the blades remained subsonic during flutter operation. The blades were found to have a higher than expected tip loading at both operating conditions. This is thought to have been caused by the outer blade sections twisting under load to higher than expected effective blade angles. This high tip loading resulted in strong vortices and a very nonuniform flow downstream of the tips of the forward-swept blades. This high tip loading may also have caused the blade flutter.
Results are presented from an investigation to measure the flow field generated by a forward-swept propfan operating in flutter at a low forward velocity. For comparison to the flutter condition, flow field data are also presented for a slightly reduced rotational speed just below flutter. The forward-swept propfan was tested as the front rotor of a counterrotating pusher propeller. A laser Doppler velocimeter (LDV) was used to measure the velocity field in planes normal to the model centerline downstream of the rotor and in planes of constant radius within the blade passages at each operating condition. A comparison of the data taken at the two different operating conditions indicated that the mean, time-averaged flow about the blades did not change drastically as the propfan rotational speed was increased from the stable operating point to the flutter condition. No regions of flow separation could be identified in the data plots of the mean intrablade flow field. The data also indicate that the relative flow about the blades remained subsonic during flutter operation. The blades were found to have a higher than expected tip loading at both operating conditions. This is thought to have been caused by the outer blade sections twisting under load to higher than expected effective blade angles. This high tip loading resulted in strong vortices and a very nonuniform flow downstream of the tips of the forward-swept blades. This high tip loading may also have caused the blade flutter.
Unsteady flow computations are being performed with the P&W (ATD) and the Rocketdyne baseline configurations of the SSME LO2 turbine turnaround duct (TAD) and heat exchanger (HEX). The work is in support of the HEX inner turning vane cracking investigation. Fatigue cracking has occurred during hot firings with the P&W configuration on the HEX inner vane, and it appears the fix will involve changes to the TAD splitter vane position and to the TAD inner wall curvature to reduce the dynamic loading on the inner vane. Unsteady flow computations on the P&W baseline and fix and on the Rocketdyne baseline reference follow steady-flow screening computations done by MSFC/ED32 on several trial configurations arriving at the fix. The P&W TAD inlet velocity profile has a strong radial velocity component that directs the flow toward the inner wall and raises the local velocity a factor of two and the dynamic pressure a factor, of four. The fix is intended to redistribute the flow more evenly across the HEX inner and outer vanes like the Rocketdyne baseline reference. Vane buffeting at frequencies around 4,000 Hz is the leading suspected cause of the problem. Our simulations (work in progress) are being done with the USA 2D axisymmetric code approximating the flow as axisymmetric u+v 2D (axial, u, and radial, v, components only). The HEX coils are included in the model to make sure the fix does not adversely affect the HEX environment. Turbulent kinetic energy, k, levels where k = 1/2 v' rms2 are locally as high as 10,000 ft2/sec2 for the P&W baseline at the engine interface (between the TAD and HEX) at the HEX inner vane location. However, k is less than 8,000 on the HEX outer vane and only about 4,500 on the HEX inner vane for the Rocketdyne baseline. Unsteady turbulence intensity, v'rms/v, and pressure, p', are being computed in the present computations to compare with steady-flow Reynolds-averaged computations where p'rms = const (pk) for overall rms random turbulence from 0.1 to 12,000 Hz frequency. Random overall static, p'rms fluctuations as large as 1.7 psi are estimated from k on the HEX inner vane for the P&W baseline configuration but only about 0.7 psi for the Rocketdyne configuration.
Explore the source record for details and available documents.