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 paper presents experience with grid ready flexible transformer unit, which is in service for two years without any difficulty. Transformer unit is capable of changing short circuit impedance on load and it is equipped with state-of-the-art monitoring system.
Sequence 7: Shroud Operating (P) This test sequence used a downwind, rigid turbine with a 3.4° cone angle. The wind speed ranged from 5 m/s to 25 m/s. Turntable angles from 0° to 30° were achieved, but the yaw error angle was maintained at 0° by adjusting the nacelle yaw. The blade tip pitch was 3°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. An airfoil-shaped aluminum shroud was placed on the tower below the 34% span five-hole probe. It extended below the 91% span five-hole probe. The turntable angle was adjusted to position the shroud at various angles relative to the flow, but the turbine yaw error angle was maintained at 0°. Two additional points were collected where the turntable angle was 0°, and the turbine yaw angle was 20° (71000200) or -20° (710M0200).
Sequence R: Step AOA, No Probes (P) This sequence was designed to quantify the effect of the five-hole probes on the 3-D blade static angle-of-attack response in the presence of rotational influences by repeating Sequence K without five-hole probes. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speeds ranged from 6 m/s to 20 m/s, and data were collected at yaw angles of 0° and 30°. The rotor rotated at 72 RPM. Blade pressure measurements were collected. The five-hole probes were removed and the plugs were installed. Plastic tape 0.03-mm-thick was used to smooth the interface between the plugs and the blade. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. During post-processing, the probe channels were set to read –99999.99. The blade pitch angle ramped continuously at 0.18°/s over a wide range of increasing and decreasing pitch angles. A step sequence was also performed. The blade pitch was stepped 5°; the flow was allowed to stabilize; and the pitch angle was held for 5 seconds. Then the pitch angle step was repeated. Again, a wide range of pitch angles was obtained, both increasing and decreasing. The file lengths for this sequence varied from 96 seconds to 6 minutes, depending on the pitch angle range. Some short points were collected at 0° yaw and 3° pitch to ascertain the functionality of the instrumentation and repeatability over time. The file name convention used the initial letter R, followed by two digits specifying wind speed, followed by two digits for yaw angle, followed by RU, RD, or ST, followed by the repetition digit. The angle of attack motion was differentiated by RU (ramp up), RD (ramp down), and ST (step down, then step up). This sequence is related to Sequences K and L.
We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. Furthermore, the findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.
Reductive catalytic fractionation (RCF) is a promising lignin-first biorefining strategy that selectively extracts and depolymerizes lignin from whole biomass. Flow-through (FT)-RCF enables physical separation of the biomass and catalyst, but this process configuration typically operates at high solvent-to-biomass ratios. Here, we demonstrate multi-pass FT-RCF, wherein the lignin-enriched solvent obtained after an initial FT-RCF step is recycled and used in subsequent FT-RCF without intermediate lignin recovery. Multi-pass FT-RCF reduces the solvent-to-biomass ratio from 48 to 1.9 L/kg, which is a lower solvent loading than is accessible in batch reactors, without negative impacts on delignification and monomer yield with up to 12 wt % lignin in the solvent. Overall, this work demonstrates that solvent demand in RCF processes, which is a key cost and energy driver to enable this process at scale, can be reduced by recycling lignin oil between RCF reactions.
Recently, there has been considerable research interest in the potential for DC distribution systems in buildings instead of the traditional AC distribution systems. Due to the need for performing power conversions between DC and AC electricity, DC distribution may provide electrical efficiency advantages in some systems. To support comparative evaluations of AC-only, DC-only, and hybrid AC/DC distribution systems in buildings, a new modeling toolkit called the Building Electrical Efficiency Analysis Model (BEEAM) was developed and is described in this paper. To account for harmonics in currents or voltages arising from nonlinear devices, the toolkit implements harmonic power flow, along with nonlinear device behavioral descriptions derived from empirical measurements. This paper describes the framework, network equations, device representations, and an implementation of the toolkit in an open source software package, including a component library and graphical interface for creating circuits. Simulations of electrical behavior and device and system efficiencies using the toolkit are compared with experimental measurements of a small office environment in a variety of operating and load configurations. A detailed analysis of uncertainty estimation is also provided. Key findings were that a comparison of predicted versus measured efficiencies and power losses in the validation testbed using the initial toolkit implementation predicted device- and system-level efficiencies with reasonably good accuracy under both balanced and unbalanced AC scenarios. An uncertainty analysis also revealed that the maximum estimated error for system efficiency across all scenarios was 3%, and measured and modeled system efficiency agreed within the experimental uncertainty in approximately half of the scenarios. Based on the correspondence between simulation and measurement, the toolkit is proposed by the authors as a potentially useful tool for comparing efficiency in AC, DC, and hybrid AC/DC distribution systems in buildings.
In the existing operations, hydropower contributions to future system flexibility are generally modeled while maintaining traditional operating rules and constraints in supporting grid operation, such as the balancing of variable renewable energy production. Moreover, operation of large scale hydropower systems on major rivers has been investigated for decades, utilizing various systems engineering approaches, with the evolving electric grid, as the result of renewable resources integration, compounded by the changing climate (variability of river flows, intensification of hydrologic cycle resulting in more frequent extreme events) affecting water availability, the need for more advanced stochastic modeling and effective uncertainty analysis approaches have become necessary. The research results supported by this funding and presented in this report provide a new look at hydropower operational flexibility enforced by the changes identified above. Understanding how hydropower operates in response to the underlying uncertainties with respect to the system constraints is crucial in identifying its operational flexibility potentials. In this project, the flexibility of the operating hydropower facility is described by capturing uncertainties in both water and power system and formulating the operations as a multistage stochastic optimization problem. The proposed approach supports short- to seasonal-term operations and planning decision horizons.
One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.
In this report the plastic deformation behavior of single crystals of the quinary, equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy (HEA) with the face-centered cubic structure has been investigated in tension and compression as a function of crystal orientation and temperature from 10 K to 1373 K. The critical resolved shear stress (CRSS) for {111}<110> slip at room temperature is 42-45 MPa. It does not depend much on crystal orientation (i.e., the Schmid law holds true) and the sense (tension vs. compression) of the applied load. The CRSS for {111}<110> slip increases with the decrease in temperature, without showing any significant inertia effects at cryogenic temperatures below 77 K. Extrapolation from the measured yield stresses down to 10 K yields a CRSS value at 0 K of 168 MPa. At cryogenic temperatures, the measured strain-rate sensitivity of flow stress is consistent with a very small activation volume. The concept of stress equivalence holds true both for the temperature dependence of CRSS and the stress dependence of activation volume, indicating that solid-solution hardening is the major strengthening mechanism. Deformation twinning occurs at 77 K but not at room temperature, resulting in higher tensile elongation to failure at 77 K than at room temperature. Deformation twinning at 77 K occurs at a shear stress of 378 MPa on conjugate ($\bar1\bar11$) planes in the form of Lüders deformation after large plastic strain (about 85%) achieved by the stage I (easy glide) and stage II (linear work-hardening) deformation.
We demonstrated production of a superior performance biodiesel referred to here as fatty acid fusel alcohol esters (FAFE) – by reacting fusel alcohols (isobutanol, 3-methyl-1-butanol, and (S)-(-)-2-methyl-1-butanol) with oil (glyceryl trioleate) using lipase from Aspergillus oryzae. Reaction conditions corresponding to a molar ratio of 5:1 (fusel alcohols to oil), enzyme loading of 2% w/w, reaction temperature of 35 °C, shaking speed of 250 rpm, and reaction time of 24 h achieved >97% conversion to FAFE. Further, FAFE obtained from reacting a fusel alcohol mixture with corn oil were evaluated for use as a fuel for diesel engines. FAFE mixtures showed superior combustion and cold-flow properties, with the derived cetane numbers up to 4.8 points higher, cloud points up to -6 °C lower, and the heat of combustion up to 2.1% higher than the corresponding FAME samples, depending on the fusel mixture used. This represents a significant improvement for all three metrics, which are typically anti-correlated. Finally, FAFE provides a new opportunity for expanded usage of biodiesel by addressing feedstock limitations, fuel performance, and low temperature tolerance.
The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.
This paper develops a transactive energy system model to restore electricity to customers in an isolated distribution system after an outage. The model engages a variety of customer types -- prosumers, flexible loads, critical/noncritical customers, and distributed generators -- as active participants in the restoration process. Unlike many existing transactive approaches, the proposed model is developed for service restoration and accounts for various customer types and their \textit{autonomy} and \textit{privacy} through an iterative approach to determine the optimal market price, while maintaining system-level power flow and voltage constraints. The advantages of the proposed approach are numerically validated on a modified IEEE 123-bus test system.
Sequence P: Wake Flow Visualization, Upwind (P) This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 5 m/s to 15 m/s. Yaw angles of 0° to –60° were achieved. The blade tip pitch was 3°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The aluminum blade tip designed to contain a smoke generator was installed, and counterweights were installed in the non-instrumented blade tip to compensate. The turbine was positioned at the appropriate yaw angle, and the smoke generator was ignited remotely. The campaign duration was 3 minutes for all tests except P1000000, which was 2 minutes. The file name convention was the standard format except for P10000A0, which indicated a 3° pitch angle. File P1000000 used a 12° pitch angle. After these two campaigns were collected, it was determined that all subsequent data should be collected with a 3° pitch angle. Pressure data were not acquired during this sequence, so all associated data values were flagged as not applicable by setting the measured values in the data file to 0.000 Pa. Corresponding pressure data are available from Sequence H for the 3° pitch angle test points. Flow visualization data obtained from wall- and ceiling-mounted video cameras were recorded to videotape. The camera locations and calibration procedures are described in Appendix J.
This project is a component of a broader effort focused on geothermal heating and cooling (GHC) with the aim of illustrating the numerous benefits of incorporating GHC and geothermal heat exchange (GHX) into community energy planning and national decarbonization strategies. To better assist private sector investment, it is currently necessary to define and assess the potential of low-temperature geothermal resources. For shallow GHC/GHX fields, there is no formal compilation of subsurface characteristics shared among industry practitioners that can improve system design and operations. Alaska is specifically noted in this work, because heretofore, it has not received a similar focus in geothermal potential evaluations as the contiguous United States. The methodology consists of leveraging relevant data to generate a baseline geospatial dataset of low-temperature resources (less than 150 degrees C) to compare and analyze information accessible to anyone trying to understand the potential of GHC/GHX and small-scale low-temperature geothermal power in Alaska (e.g., energy modelers, communities, planners, and policymakers). Importantly, this project identifies data related to (1) the evaluation of GHC/GHX in the shallow subsurface, and (2) the evaluation of low-temperature geothermal resource availability. Additionally, data is being compiled to assess repurposing of oil and gas wells to contribute co-produced fluids toward the geothermal direct use and heating and cooling resource potential. In this work we identified new data from three different datasets of isolated geothermal systems in Alaska and bottom-hole temperature data from oil and gas wells that can be leveraged for evaluation of low-temperature geothermal resource potential. The goal of this project is to facilitate future deployment of GHC/GHX analysis and community-led programs and update the low-temperature geothermal resources assessment of Alaska. A better understanding of shallow potential for GHX will improve design and operations of highly efficient GHC systems. The deployment and impact that can be achieved for low-temperature geothermal resources will contribute to decarbonization goals and facilitate widespread electrification by shaving and shifting grid loads.
Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.
Cross-seam high voltage DC (HVDC) transmission in the United States can provide an efficient and flexible bulk power delivery highway between the interconnections. Recent studies have shown significant economic benefits of cross-seam HVDC transmission systems by taking the advantages of load diversity, frequency response, renewable energy diversity, and energy arbitrage in the different interconnections. This paper introduces a station-hybrid system scheme for cross-seam interconnections. Here, the proposed station-hybrid system combines the advantages of line commutated converter (LCC) and voltage source converter (VSC) technologies, thus enabling reliable and flexible bidirectional power flows across the interconnections. Control strategies of the station-hybrid system are investigated with a focus on power flow reversal control under normal operation conditions and grid support control under emergency conditions. The feasibility and effectiveness of the proposed station-hybrid system and corresponding control strategies are verified by simulations in PSCAD/EMTDC.
This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.
To curb the energy consumption of buildings and their related CO2 emissions, Oak Ridge National Laboratory (ORNL) has developed the thermally anisotropic building envelope (TABE) —a multi-layer design comprising insulation materials and metal foils connected to thermal loops. In this study, a machine learning-assisted framework was developed to control the TABE in residential buildings to reduce the computation load for future optimal rule-based control and application. First, a 2D finite element model was established in COMSOL to calculate the hourly heat flux through exterior walls installed with the TABE. Then, TABE wall heat fluxes were simulated for various given indoor and outdoor boundary conditions, as well as thermal loops fluid temperatures and flow rates. Since the finite element simulations are computationally expensive, an artificial neural network (ANN) was then trained to use as a proxy of the finite element (COMSOL) modeling. Finally, the trained ANN model was coupled with the EnergyPlus model to predict the energy consumption of a US Department of Energy prototype single-family house installed with the TABE. An optimal simple rule-based control was determined from predefined rules for a case study. The results demonstrate that the developed machine learning–assisted framework can reduce 99.9% of the computation time while efficiently managing residential building energy for installed TABE walls.