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.
After performing steady-state Computational Fluid Dynamics (CFD) calculations using OVERFLOW to validate the CFD method against static wind-tunnel data of a box-shaped cargo container, the same setup was used to investigate unsteady flow with a moving body. Results were compared to flight test data previously collected in which the container is spinning.
Carbon nanotubes form a new class of nanomaterials that are presumed to have extraordinary mechanical, electrical and thermal properties. The single wall nanotubes (SWNTs) are estimated to be 100 times stronger than steel with 1/6th the weight; electrical carrying capacity better than copper and thermal conductivity better than diamond. Applications of these SWNTs include possible weight reduction of aerospace structures, multifunctional materials, nanosensors and nanoelectronics. Double pulsed laser vaporization process produces SWNTs with the highest percentage of nanotubes in the output material. The normal operating conditions include a green laser pulse closely followed by an infrared laser pulse. Lasers ab late a metal-containing graphite target located in a flow tube maintained in an oven at 1473K with argon flow of 100 sccm at a 500 Torr pressure. In the present work a number of production runs were carried out, changing one operating condition at a time. We have studied the effects of nine parameters, including the sequencing of the laser pulses, pulse separation times, laser energy densities, the type of buffer gas used, oven temperature, operating pressure, flow rate and inner flow tube diameters. All runs were done using the same graphite target. The collected nanotube material was characterized by a variety of analytical techniques including scanning electron microscopy (SEM), transmission electron microscopy (TEM), Raman and thermo gravimetric analysis (TGA). Results indicate trends that could be used to optimize the process and increase the efficiency of the production process.
Calculations for ice accretion and shedding are presented for a model scale rotor in hover. The aerodynamic characteristics of the rotor are first computed using a combined blade element-momentum theory. The effective angles of attack, and the local flow velocity are used within the NASA Glenn solver LEWICE to estimate the collection efficiency. The computed convection efficiency and the surface pressure distribution from a panel method within LEWICE are used to estimate the ice accretion over the rotor blades over an elapsed time interval. Finally, a force balance approach is used to establish shedding events where the centrifugal force over the ice mass exceeds the adhesive forces at the rotor surface and the cohesive forces between adjacent masses of ice. Preliminary comparisons with test data acquired at the Pennsylvania State Icing Research Tunnel are presented. Sensitivity of the ice shedding events to surface roughness, adhesive strength, cohesive strength, and ambient conditions is discussed.
Fabrication of functional polymer-based particles by crosslinking UV-curable polymer drops in mid-air and collecting crosslinked particles in a solid container, a liquid suspension, or an air flow. The particles can contain different phases in the form or layered structures that contain one to multiple cores, or structures that are blended with dissolved or emulsified smaller domains. A curing system produces ultraviolet rays that are directed onto the particles in the jet stream from one side. A reflector positioned on other side of the jet stream reflects the ultraviolet rays back onto the particles in the jet stream.
Fabrication of functional polymer-based particles by crosslinking UV-curable polymer drops in mid-air and collecting crosslinked particles in a solid container, a liquid suspension, or an air flow. The particles can contain different phases in the form or layered structures that contain one to multiple cores, or structures that are blended with dissolved or emulsified smaller domains. A curing system produces ultraviolet rays that are directed onto the particles in the jet stream from one side. A reflector positioned on other side of the jet stream reflects the ultraviolet rays back onto the particles in the jet stream.
Measurements and simulations show that plasma relaxation processes in the reversed field pinch drive and redistribute both magnetic flux and momentum. To examine this relaxation process, a new 3D Mach B-dot probe has been constructed. This probe collects ion saturation currents through six molybdenum electrodes arranged on the flattened vertices of an octahedron made of boron nitride (BN). The ion saturation current flows through configurable voltage dividers for measurement and returns through one of six selectable return electrodes equally spaced along the 12 cm BN probe arm. In addition, the probe arm houses three B-dot magnetic pickup coils in the BN stalk immediately below to the octahedron, to measure the local magnetic field. Inserted in the Madison Symmetric Torus (MST) during deuterium discharges with 220 kA plasma current, density of 0.8 × 10 13 cm –3 , the probe collects ion saturation currents with sawtooth-like peaks correlated with relaxation events. This compact octahedral design fitting six Mach electrode surfaces within a 1 cm3 cube will enable future multi-point, multi-field probes compatible with the 1.5 in. ports of MST. Such probes will allow for flow circulation, current, and canonical vorticity to be calculated in the center of the finite difference stencil formed by the measurement locations.
Temporally Resolved Particle Image Velocimetry (TR-PIV) is being used to characterize the decay of turbulence in jet flows a critical element for understanding the acoustic properties of the flow. A TR-PIV system, developed in-house at the NASA Glenn Research Center, is capable of acquiring planar PIV image frame pairs at up to 10 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered 6 nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centered along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space-time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number.
Acceleration data taken from the Orbital Acceleration Research Experiment (OARE) during reentry on STS-62 have been analyzed using calibration factors taken on orbit. This is the first Orbiter mission which collected OARE data during the Orbiter reentry phase. The data examined include the flight regime from orbital altitudes down to about 90 km which covers the free-molecule-flow regime and the upper altitude fringes of the rarefied-flow transition into the hypersonic continuum. Ancillary flight data on Orbiter position, orientation, velocity, and rotation rates have been used in models to transform the measured accelerations to the Orbiter center-of-gravity, from which aerodynamic accelerations along the Orbiter body axes have been calculated. Residual offsets introduced in the measurements by unmodeled Orbiter forces are identified and discussed. Direct comparisons are made between the OARE flight data and an independent micro-gravity accelerometer experiment, the High Resolution Accelerometer Package (HiRAP), which also obtained flight data on reentry during the mission down to about 95 km. The resulting OARE aerodynamic acceleration measurements along the Orbiter's body axis, aid the normal to axial acceleration ratio in the free-molecule-flow and transition-flow regimes are presented and compared with numerical simulations from three direct simulation Monte Carlo codes.
This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }
Temporally Resolved Particle Image Velocimetry (TR-PIV) is the newest and most exciting tool recently developed to support our continuing efforts to characterize and improve our understanding of the decay of turbulence in jet flows -- a critical element for understanding the acoustic properties of the flow. A new TR-PIV system has been developed at the NASA Glenn Research Center which is capable of acquiring planar PIV image frame pairs at up to 25 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered 6 nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centered along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space-time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number
Temporally resolved particle image velocimetry (TR-PIV) is the newest and most exciting tool recently developed to support our continuing efforts to characterize and improve our understanding of the decay of turbulence in jet flows—a critical element for understanding the acoustic properties of the flow. A new TR-PIV system has been developed at the NASA Glenn Research Center which is capable of acquiring planar PIV image frame pairs at up to 25 kHz. The data reported here were collected at Mach numbers of 0.5 and 0.9 and at temperature ratios of 0.89 and 1.76. The field of view of the TR-PIV system covered six nozzle diameters along the lip line of the 50.8 mm diameter jet. The cold flow data at Mach 0.5 were compared with hotwire anemometry measurements in order to validate the new TR-PIV technique. The axial turbulence profiles measured across the shear layer using TR-PIV were thinner than those measured using hotwire anemometry and remained centred along the nozzle lip line. The collected TR-PIV data illustrate the differences in the single point statistical flow properties of cold and hot jet flows. The planar, time-resolved velocity records were then used to compute two-point space–time correlations of the flow at the Mach 0.9 flow condition. The TR-PIV results show that there are differences in the convective velocity and growth rate of the turbulent structures between cold and hot flows at the same Mach number.
A system for detecting MITM for SCADA communication networks includes secure substation-substation communication links for providing secure and reliable paths to exchange OT data between substations for OT data consistency check; a SIB in each substation for sampling CT and PT measurements to calculate voltage magnitude and phase angle thereof; a S&C server in each substation coupled to the SIB for receiving the voltage magnitude and phase angle from the SIB and obtaining a packet carrying active power flow in transmission lines between two substations and a time stamp; an IDS server placed in a SCADA center for collecting the packet of each substation sent by the S&C server; analyzing the received packet from every adjacent substation; inspecting the payload of the received packet; and triggering an intrusion alarm to a SCADA operator when the power flow is not the same as the payload of the packets.
An investigation of leading-edge contamination by insects was conducted at the NASA Dryden Flight Research Center with a JetStar airplane instrumented to detect transition on the outboard leading-edge flap and equipped with a system to wash the leading edge in flight. The results of airline-type flights with the JetStar indicated that insects can contaminate the leading edge during take-off and climbout at large jet airports in the United States. The results also showed that the insects collected on the leading edges at 180 knots did not erode at cruise conditions for a laminar flow control airplane and caused premature transition of the laminar boundary layer. None of the superslick and hydrophobic surfaces tested showed any significant advantages in alleviating the insect contamination problem. While there may be other solutions to the insect contamination problem, the results of these tests with a washer system showed that a continuous water spray while encountering the insects is effective in preventing insect contamination of the leading edges.
As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.
This investigation of vapor-liquid phase separation (VLPS) of He 2 is related to long-term storage of cryogenic liquid. The VLPS system utilizes porous plugs in order to generate thermomechanical (thermo-osmotic) force which in turn prevents liquid from flowing out of the cryo-vessel (e.g., Infrared Astronomical Satellite). An apparatus was built and VLPS data were collected for a 2 and a 10 micrometer sintered stainless steel plug and a 5 to 15 micrometer sintered bronze plug. The VLPS data obtained at high temperature were in the nonlinear turbulent regime. At low temperature, the Stokes regime was approached. A turbulent flow model was developed, which provides a phenomenological description of the VLPS data. According to the model, most of the phase separation data are in the turbulent regime. The model is based on concepts of the Gorter-Mellink transport involving the mutual friction known from the zero net mass flow (ZNMF) studies. The latter had to be modified to obtain agreement with the present experimental VLPS evidence. In contrast to the well-known ZNMF mode, the VLPS results require a geometry dependent constant (Gorter-Mellink constant). A theoretical interpretation of the phenomenological equation for the VLPS data obtained, is based on modelling of the dynamics of quantized vortices proposed by Vinen. In extending Vinen's model to the VLPS transport of He 2 in porous media, a correlation between the K*(GM) and K(p) was obtained which permits an interpretation of the present findings. As K(p) is crucial, various methods were introduced to measure the permeability of the porous media at low temperatures. Good agreement was found between the room temperature and the low temperature K(p)-value of the plugs.
The setup and initial operation of a set of specialized meteorological data collection hardware are described. To study the life cycle of advection fogs at a lake test site, turbulence levels in the fog are identified, and correlated with the temperature gradients and mean wind profiles. A meteorological tower was instrumented to allow multiple-level measurements of wind and temperature on a continuous basis. Additional instrumentation was: (1)hydrothermograph, (2)microbarograph, (3)transmissometers, and (4)a boundary layer profiler. Two types of fogs were identified, and important differences in the turbulence scales were noted.
In the context of managing distributed energy resources (DERs) within distribution networks (DNs), this work focuses on the task of developing local controllers. We propose an unsupervised learning framework to train functions that can closely approximate optimal power flow (OPF) solutions. The primary aim is to establish specific conditions under which these learned functions can collectively guide the network towards desired configurations asymptotically, leveraging an incremental control approach. The flexibility of the proposed methodology allows to integrate fairness-driven components into the cost function associated with the OPF problem. This addition seeks to mitigate power curtailment disparities among DERs, thereby promoting equitable power injections across the network. To demonstrate the effectiveness of the proposed approach, power flow simulations are conducted using the IEEE 37-bus feeder. The findings not only showcase the guaranteed system stability but also underscore its improved overall performance.