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 National Renewable Energy Laboratory (NREL) is analyzing the rapidly increasing role of energy storage in the electrical grid through 2050 through its Storage Futures Study. In one phase of the study, NREL used the laboratory's Distributed Generation Market (dGen) model to examine the various future distributed storage capacity adoption scenarios, results, and implications.
Abstract Chemomechanical weakening of layered phyllosilicate muscovite mica was studied as a function of chemical environment via in situ liquid-nanoindentation under four conditions (dry, deionized water, and two NaCl solutions of different pH). While traditional indentation analyses of layered materials with extreme mechanical anisotropy have been limited due to displacement bursts (pop-ins), here the bursts were used as proxies for delamination, fracture, and spalling events. Since displacement bursts during an indent represent a stochastic process, 120 indents were conducted for each condition to generate statistically significant amounts of data. In total, over 9000 bursts were assessed using a load–displacement threshold criterion, classifying this as a high-throughput nanoscale fracture technique. For each burst, initiation load, initiation displacement, plastic zone volume at initiation, and energy dissipation were analyzed. A power-law relationship between the burst load and burst energy was noted which separated the bursts into two continuous distributions: (1) bursts due only to the mechanics of the indent and (2) bursts due to both the mechanics of the indent and the environment. By using a cumulative probability distribution, it was found that the NaCl solutions decreased the minimum plastic zone volume necessary to initiate a displacement burst by an order of magnitude relative to the dry condition. Finally, the underlying mechanisms explaining the trends in initiation volume as a function of environment were discussed, with a focus on the chemomechanical degradation processes via chemical attack and cation exchange.
Microgrids and microgrid technologies offer an operations and control framework for realizing the needed aggregation of a growing share of distributed energy resources. As the number of microgrids increases, so will the incentive to interconnect them, leading to a combinatorial growth of potential operating states. At the same time, their growing complexity nudges interconnected microgrids conceptually closer to balancing authorities in the bulk power system. This paper focuses on secondary frequency control and leverages the conceptual similarity between balancing authorities and interconnected microgrids to adopt an automatic generation control design with some adaptations. Frequency control and restoration is demonstrated via several use-cases with switching and reconfiguration events.
H2Integrate(H2I) is an open-source Python package for hybrid systems engineering design and technoeconomic analysis. It models and optimizes hybrid energy plants that produce electricity, hydrogen, ammonia, steel, and other products. H2Integrate is designed to be flexible and extensible, allowing users to create their own components and models for various energy systems. The tool currently includes distributed energy generation (wind, solar, wave, tidal), battery storage, hydrogen, ammonia, methanol, and steel technologies. Other elements such as desalination systems, pipelines, compressors, and storage systems can also be included as developed by users. Some modeling capabilities in H2Integrate are provided by integrating existing tools, such as HOPP, PySAM, ORBIT, and ProFAST. The H2Integrate tool is built on top of NASA's OpenMDAO framework, which provides a powerful and flexible environment for modeling and optimization.
The control and protection functions of inverter-based resources (IBR) have raised concerns with bulk system reliability. Most of the current interest lies with solar photovoltaic generation but increasing amounts of storage would pose the same risks. Newer North American Electric Reliability Corporation (NERC) guidelines call for electromagnetic transient (EMT) studies of IBR and recommend that transmission operators collect distributed energy resource (DER) data to support such modeling. IEEE Standard P2800.1 is defining tests for model parameterization, so good model data should become available from inverter vendors. (EMT studies also apply to large power transformer reliability, and transformer vendors can provide EMT models.) Utilities don’t currently have the rest of the bulk system represented for EMT studies at large scale. An International Electrotechnical Commission (IEC) standard Common Information Model (CIM) provides a way of supporting these detailed models from physical asset data, e.g., conductors, towers, transformer data sheets, control block diagrams, while avoiding software vendor lock-in. CIM-for-EMT, with proposed schema extensions and open-source converters, provides a way to exchange EMT data between organizations and tools. This project leverages Office of Electricity (OE) funding of CIM-for-EMT code base through the GridAPPS-DTM project, and of GridPACKTM (parallelized transmission solver), for interoperability testing in CIM-for-EMT. The project also leverages partner PGSTech investments in EMTP® interoperability with CIM.
The application is considered of vortex lattice techniques to the problem of describing the aerodynamics and performance of statically thrusting propellers. A numerical lifting surface theory to predict the aerodynamic forces and power is performed. The chordwise and spanwise loading is modelled by bound vortices fixed to a twisted flat plate surface. In order to eliminate any apriori assumptions regarding the wake shape, it is assumed the propeller starts from rest. The wake is generated in time and allowed to deform under its own self-induced velocity field as the motion of the propeller progresses. The bound circulation distribution is then determined with time by applying the flow tangency boundary condition at certain selected control points on the blades. The aerodynamics of the infinite wing and finite wing are also considered. The details of wake formation and roll-up are investigated, particularly the localized induction effect. It is concluded that proper wake roll-up and roll-up rates can be established by considering the details of motion at the instant of start.
A study of Voyager 1 electric field measurements obtained by the plasma wave instrument in the Io plasma torus has been carried out. A survey of the data has revealed the presence of persistent peaks in electric field spectra in the frequency range 100-600 Hz consistent with their identification as lower hybrid noise for a heavy-ion plasma of sulfur and oxygen. Typical wave intensities are 0.1 mV/m, and the spectra also show significant Doppler broadening, Delta omega/omega approximately 1. A theoretical analysis of lower hybrid wave generation by a bump-on-tail ring distribution of ions is given. The model is appropriate for plasmas with a superthermal pickup ion population present. A general methodology is used to demonstrate that the maximum plasma heating rate possible through anomalous wave-particle heat exchange is less than approximately 10 to the -14th ergs per cu cm per s. Although insufficient to meet the power requirement of the EUV-emitting warm torus, the heating rate is large enough to maintain a low-density (0.01-0.1 percent) superthermal electron population of keV electrons, which may lead to a small but significant anomalous ionization effect.
Fourier analysis is applied to very high resolution image of coronal active regions obtained by the Normal Incidence X-Ray Telescope is used to find a broad isotropic power-law spectrum of the spatial distribution of soft X-ray intensities. Magnetic structures of all sizes are present down to the resolution limit of the instrument. Power spectra for the X-ray intensities of a sample of topologically different active regions are found which fall off with increasing wavenumber as 1/k-cubed. A model is presented that relates the basic features of coronal magnetic fluctuations to the subphotospheric hydrodynamic turbulence that generates them. The model is used to find a theoretical power spectrum for the X-ray intensity which falls off with increasing wavenumber as 1/k-cubed. The implications of a turbulent regime in active regions are discussed.
It is unlikely that the presence of chondrules, and thus their formation, within the protoplanetary nebula would be predicted if it were not for their ubiquitous presence in most chondritic meteorites. The study of these enigmatic, igneous objects has a direct influence on how meteoritic and solar system researchers model the processes operating and the materials present within our protoplanetary nebula. Key to understanding chondrule formation is a determination of constraints on their thermal histories. The three important variables in this history are their peak melting temperatures, the duration of their melting at peak temperatures, and the rate at which these object cool. Although these three variables are interdependent, it is cooling rate that provides the most powerful constraint. Cooling rate has a direct affect on the development of both crystal morphology and the elemental distributions within these grains. To date, experiments have indicated that chondrule cooling rates are in the range of 10's to 100's of degrees per hour for porphyritic chondrules (the most abundant type). The cooling rate for radial and barred chondrules is thought to be more rapid. To generate these cooling rates (rapid relative to the cooling of the nebula as a whole, but slow compared to simple black body radiation) the environment of chondrule formation must have been localized, and the abundance of solid materials must have been greatly enhanced above a gas of solar composition. Thus accurate determinations of chondrule cooling rates is critical in understanding both their formation and the nebular environment in which they formed. In a quest to more accurately determine the lower limit on cooling rates and to determine in more detail the effects of Fe loss from a molten sample to Pt wire loops, Weinbruch et al. have explored this issue experimentally and reevaluated the findings of Radomsky and Hewins in light of their new results. The basic conclusions of their paper are an important contribution to our understanding of how experimental techniques can affect established constraints on chondrule formation and are thus of interest to a wide audience. We do believe, however, that their methodology produces results that provide inappropriate impressions of the applicability of their study to chondrule formation and nebular processes. Furthermore, the extensive body of previous experimental work on chondrule bulk compositions cannot be invalidated, as they suggest, by the results of Weinbruch et al. It is for these reasons that within this comment we address to applicability of the results presented by Weinbruch et al. to previous studies, and illustrate how the experimental conditions chosen for their series of experiments introduced a significant bias in their results.
Rapid efficient access to Earth sciences data is fundamental to the Nation's efforts to understand the effects of global environmental changes and their implications for public policy. It becomes a bigger challenge in the future when data volumes increase further and missions with constellations of satellites start to appear. Demands on data storage, data access, network throughput, processing power, and database and information management are increased by orders of magnitude, while budgets remain constant and even shrink. The Global Change Data Center's (GCDC) mission is to provide systems, data products, and information management services to maximize the availability and utility of NASA's Earth science data. The specific objectives are (1) support Earth science missions be developing and operating systems to generate, archive, and distribute data products and information; (2) develop innovative information systems for processing, archiving, accessing, visualizing, and communicating Earth science data; and (3) develop value-added products and services to promote broader utilization of NASA Earth Sciences Enterprise (ESE) data and information. The ultimate product of GCDC activities is access to data and information to support research, education, and public policy.
To support the Revolutionary Aeropropulsion Concept Program, NASA Glenn Research Center' s Structural Mechanics and Dynamics Branch is developing a compact, nonpolluting, bearingless electric machine with electric power supplied by fuel cells for future more-electric aircraft. The use of such electric drives for propulsive fans or propellers depends on the successful development of ultra-high-power-density machines that can generate power densities of 50 hp/lb or more, whereas conventional electric machines generate usually 0.2 hp/lb. One possible candidate for such ultra-high-power-density machines, a round-rotor synchronous machine with an engineering current density as high as 20 000 A/cm2 was selected to investigate how much torque and power can be produced. A simple synchronous machine model that consists of rotor and stator windings and back-irons was considered first. The model had a sinusoidally distributed winding that produces a sinusoidal distribution of flux P poles. Excitation of the rotor winding produced P poles of rotor flux, which interacted with the P stator poles to produce torque.
The mantle components that represent the source region of ocean island basalts (OIB) and feed hotspot volcanism are predicted to contain 160 +/- 20 (2σ(sub m)) ng/g Th, a heat‐producing element. This critical model composition indicates that the OIB source region (OSR) comprises a significant amount of recycled oceanic crust and constitutes 19(sup +3)(sub −2) (2σ(sub m))% of the mantle by mass. The mass fraction of this reservoir supports a mantle architecture with a basal thermochemical layering at an average depth of 2000 +/- 100 (2σ(sub m)) km or two thermochemical piles that extend up to midmantle levels. The hotspot source described here generates 10 pW/kg of radiogenic heat and supplies 7.3 TW to the planet's total surface heat flux. Given that the silicate portion of the Earth produces some 20.4 TW of radiogenic power, with 7.2 TW derived from the continental crust, the mantle source responsible for mid‐ocean ridge volcanism provides only 5.9 TW of radiogenic power (or <2 pW/kg). As a result, the source of hotspots generates >5x more radiogenic heat than the source of mid‐ocean ridges, thus contributing to the energetics that drive mantle convection and potentially the formation of long‐lived plumes via bottom heating of the modern mantle. The potential for a sequestered or unsampled mantle reservoir would impact the relative mass fractions of the source regions of OIB and mid‐ocean ridge volcanism but not the compositional model of the OSR presented here.
CalderaCast is a user-friendly web-based tool for electrical-load forecast, providing stakeholders with a fully customizable decision-support framework that estimates the likely power draw from a possible future electric-vehicle (EV) charging station at a given location on a given day along an alternative fuel corridor (AFC). These EV charging profiles are accurately modeled in CalderaCast using the Caldera software framework developed by Idaho National Laboratories (INL), reflecting the realistic charging levels observed in actual charge events. This tool was developed as part of the National Electric Vehicle Infrastructure (NEVI) program, which is quickly generating substantial interest from would-be charging station operators (CSO), large and small electric utilities, and state transportation planners, some of whom had not seriously considered EV charging previously. All these entities—with or without background in EV infrastructure—must estimate the electricity load that a proposed charging station will generate. This load forecast is critically important for a utility to properly assess the capacity of their distribution network to support the proposed station or properly size grid upgrades for potential load growth due to future EV adoption, vehicle technology improvements, or station growth. This document describes each aspect of the CalderaCast tool and provides guidance to users who are interested in utilizing the tool for their work.
With large-scale distributed photovoltaics (PVs) being integrated into distribution networks (DNs), coordinated optimal power flow (OPF) of high voltage (HV) and medium voltage (MV) DNs should be investigated to optimally dispatch the distributed PVs and other network devices. Here, this paper presents a hierarchical distributed OPF method for HV and MV DNs with on-load tap changers, reactive power compensators, feeder switches and distributed PVs. A hierarchical master-slave control architecture is applied to implement coordinated OPF of two-layer DNs. The HV master problem and MV subproblems are transformed into mixed-integer convex problems respectively with second order cone programming and LinDistFlow approximation. Since there is no efficient distributed algorithm to solve such OPF models with integer subproblems, a novel distributed algorithm is proposed in this paper to efficiently solve the hierarchical coordinated OPF model with integer subproblems in a distributed manner. In the proposed algorithm, the coordinated OPF model is solved in a branch-and-bound framework, where in each branch node generalized Benders decomposition (GBD) algorithm is applied to decompose the coordinated OPF model into a master problem and relaxed subproblems and solves them iteratively to get optimal solution. The GBD optimal and feasible cutting planes generated in a branch node are proved to be valid for its descendants. Moreover, three acceleration techniques are introduced into the proposed algorithm to improve computational efficiency. Finally, the effectiveness and accuracy of the proposed method are verified via simulation tests in Jinzhai DNs of China.
The integrated hydrogen energy system incorporates hydrogen energy into the power grid, which has been recognized as a promising option for reaching a 100% renewable electricity supply. It can make a profit because the hydrogen produced can be sold as fuel or used to generate electricity for grid services. In this paper, we develop a planning model for the integrated hydrogen energy system that considers the uncertainty of the load demand, the renewable energy generation, and the market prices. To calculate the hydrogen load, we simulate the refueling operations at a hydrogen fueling station over the course of one day. Moreover, the long-term profitability of the integrated system under both cur-rent and future conditions is compared in the planning results.
In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.
The Fast Modular Reactor (FMR) is a 100-MW(thermal) gas-cooled fast reactor being developed by General Atomics Electromagnetic System with the goal of developing a FMR for flexible and dispatchable power to the U.S. electricity market in the mid-2030s. The conceptual design aims to develop and verify simplified design features. These include an inert helium gas coolant, pellet-loaded fuel rods, installations with air cooling as ultimate heat sink, and small and passive heat removal systems. The goal is to ensure the development of a safe, maintainable, cost-effective, and distributed nuclear energy-generating station. The baseline technologies selected to achieve this goal are a helium coolant that is an inert gas with no chemical reaction with structural components, not activated, single phase, enabling high-temperature operation and a high thermal efficiency Brayton cycle; conventional uranium dioxide (UO 2 ) fuel, which is the most widely used and well-known fuel material, capable of high burnup (100 MWd/kg) and a long fuel life; and silicon carbide composite (SiGA®) cladding and internal structures that are chemically inert in the helium environment, exceptionally radiation tolerant, and being derisked by accident tolerant fuel technology development. Further, the reactor was specifically designed with passive safety features, including high-temperature in-core materials and a reactor vessel cooling system consisting of cooling panels of naturally circulating water. The passive safety of the core was confirmed for the depressurized loss-of–forced cooling accident, which showed the peak cladding temperature at ~1600°C during the transient, which is below the current design limit of 1800°C. The conceptual design of the FMR has been conducted for the reactor system, vessel system, generator and turbomachine, instrumentation and control, residual heat removal system, plant service system, and containment, as well as pre-application licensing documents.