SEARCH · Search NASA

Results for “grid simulator testing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Unsteady Aerodynamic Force Sensing from Measured Strain

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm. A cantilevered rectangular wing built and tested at the NASA Langley Research Center (Hampton, Virginia, USA) in 1959 is used to validate the simple approach. Unsteady aerodynamic forces as well as wing deflections, velocities, accelerations, and strains are computed using the CFL3D computational fluid dynamics (CFD) code and an MSC/NASTRAN code (MSC Software Corporation, Newport Beach, California, USA), and these CFL3D-based results are assumed as measured quantities. Based on the measured strains, wing deflections, velocities, accelerations, and aerodynamic forces are computed using the proposed approach. These computed deflections, velocities, accelerations, and unsteady aerodynamic forces are compared with the CFL3D/NASTRAN-based results. In general, computed aerodynamic forces based on the lifting surface theory in subsonic speeds are in good agreement with the target aerodynamic forces generated using CFL3D code with the Euler equation. Excellent aeroelastic responses are obtained even with unsteady strain data under the signal to noise ratio of -9.8dB. The deflections, velocities, and accelerations at each sensor location are independent of structural and aerodynamic models. Therefore, the distributed strain data together with the current proposed approaches can be used as distributed deflection, velocity, and acceleration sensors. This research demonstrates the feasibility of obtaining induced drag and lift forces through the use of distributed sensor technology with measured strain data. An active induced drag control system thus can be designed using the two computed aerodynamic forces, induced drag and lift, to improve the fuel efficiency of an aircraft. Interpolation elements between structural finite element grids and the CFD grids and centroids are successfully incorporated with the unsteady aeroelastic computation scheme. The most critical technology for the success of the proposed approach is the robust on-line parameter estimator, since the least-squares curve fitting method depends heavily on aeroelastic system frequencies and damping factors.

shape sensing↗

DOE OSTI2022

The Application of a Genetic Algorithm to the Optimization of a Mesoscale Model for Emergency Response

Besides solving the equations of momentum, heat, and moisture transport on the model grid, mesoscale weather models must account for subgrid-scale processes that affect the resolved model variables. These are simulated with model parameterizations, which often rely on values preset by the user. Such “free” model parameters, along with others set to initialize the model, are often poorly constrained, requiring that a user select each from a range of plausible values. Finding the values to optimize any forecasting tool can be accomplished with a search algorithm, and one such process—the genetic algorithm (GA)—has become especially popular. As applied to modeling, GAs represent a Darwinian process: an ensemble of simulations is run with a different set of parameter values for each member, and the members subsequently judged to be most accurate are selected as “parents” who pass their parameters onto a new generation. At the U.S. Department of Energy’s Savannah River Site in South Carolina, we are applying a GA to the Regional Atmospheric Modeling System (RAMS) mesoscale weather model, which supplies input to a model to simulate the dispersion of an airborne contaminant as part of the site’s emergency response preparations. An ensemble of forecasts is run each day, weather data are used to “score” the individual members of the ensemble, and the parameters from the best members are used for the next day’s forecasts. As meteorological conditions change, the parameters change as well, maintaining a model configuration that is best adapted to atmospheric conditions. Significance Statement We wanted to develop a forecasting system by which a weather model is run over the Savannah River Site each day and repeatedly adjusted according to how well it performed the previous day. To run the model, a series of values (parameters) must be set to control how the model will calculate winds, temperatures, and other desired variables. Each day the model was run several times using different combinations of these parameters and later compared with observed meteorological conditions. Parameters that produced the most accurate forecasts were preferentially reused to create the forecasts for the next day. The process was tested for the summer of 2020 and exhibited lower errors than forecasts produced by the model using default values of the parameters.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗

NASA NTRS1993

Transient liquid-crystal technique used to produce high-resolution convective heat-transfer-coefficient maps

In this transient technique the preheated isothermal model wall simulates the classic one-dimensional, semi-infinite wall heat transfer conduction problem. By knowing the temperature of the air flowing through the model, the initial temperature of the model wall, and the surface cooling rate measured at any location with time (using the fast-response liquid-crystal patterns recorded on video tape), the heat transfer coefficient can be calculated for the color isothermal pattern produced. Although the test was run transiently, the heat transfer coefficients are for the steady-state case. The upstream thermal boundary condition was considered to be isothermal. This transient liquid-crystal heat-transfer technique was used in a transient air tunnel in which a square-inlet, 3-to-1 exit transition duct was placed. The duct was preheated prior to allowing room temperature air to be suddenly drawn through it. The resulting isothermal contours on the duct surfaces were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the temperature and time data for all points on the duct surfaces during each test. The duct surfaces were uniformly heated using two heating systems: the first was an automatic temperature-controlled heater blanket completely surrounding the test duct like an oven, and the second was an internal hot-air loop through the inside of the test duct. The hot-air loop path was confined inside the test duct by insulated heat dams located at the inlet and exit ends of the test duct. A recirculating fan moved hot air into the duct inlet, through the duct, out of the duct exit, through the oven, and back to the duct inlet. The temperature nonuniformity of the test duct model wall was held very small. Test results are reported for two inlet Reynolds numbers of 200,000 and 1,150,000 (based on the square-inlet hydraulic diameter) and two free-stream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 20 percent (using a grid), which is typical of real engine conditions.

Hippensteele, Steven A.↗

NASA NTRS2011

NASA Tech Briefs, May 2011

Topics covered include: 1) Method to Estimate the Dissolved Air Content in Hydraulic Fluid; 2) Method for Measuring Collimator-Pointing Sensitivity to Temperature Changes; 3) High-Temperature Thermometer Using Cr-Doped GdAlO3 Broadband Luminescence; 4)Metrology Arrangement for Measuring the Positions of Mirrors of a Submillimeter Telescope; 5) On-Wafer S-Parameter Measurements in the 325-508-GHz Band; 6) Reconfigurable Microwave Phase Delay Element for Frequency Reference and Phase-Shifter Applications; 7) High-Speed Isolation Board for Flight Hardware Testing; 8) High-Throughput, Adaptive FFT Architecture for FPGA-Based Spaceborne Data Processors; 9) 3D Orbit Visualization for Earth-Observing Missions; 10) MaROS: Web Visualization of Mars Orbiting and Landed Assets; 11) RAPID: Collaborative Commanding and Monitoring of Lunar Assets; 12) Image Segmentation, Registration, Compression, and Matching; 13) Image Calibration; 14) Rapid ISS Power Availability Simulator; 15) A Method of Strengthening Composite/Metal Joints; 16) Pre-Finishing of SiC for Optical Applications; 17) Optimization of Indium Bump Morphology for Improved Flip Chip Devices; 18) Measuring Moisture Levels in Graphite Epoxy Composite Sandwich Structures; 19) Marshall Convergent Spray Formulation Improvement for High Temperatures; 20) Real-Time Deposition Monitor for Ultrathin Conductive Films; 21) Optimized Li-Ion Electrolytes Containing Triphenyl Phosphate as a Flame-Retardant Additive; 22) Radiation-Resistant Hybrid Lotus Effect for Achieving Photoelectrocatalytic Self-Cleaning Anticontamination Coatings; 23) Improved, Low-Stress Economical Submerged Pipeline; 24) Optical Fiber Array Assemblies for Space Flight on the Lunar Reconnaissance Orbiter; 25) Local Leak Detection and Health Monitoring of Pressurized Tanks; 26) Dielectric Covered Planar Antennas at Submillimeter Wavelengths for Terahertz Imaging; 27) Automated Cryocooler Monitor and Control System; 28) Broadband Achromatic Phase Shifter for a Nulling Interferometer; 29) Super Dwarf Wheat for Growth in Confined Spaces; 30) Fine Guidance Sensing for Coronagraphic Observatories; 31) Single-Antenna Temperature- and Humidity-Sounding Microwave Receiver; 32) Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping; 33) Optical Communications Link to Airborne Transceiver; 34) Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings; 35) Entanglement in Self-Supervised Dynamics; 36) Prioritized LT Codes; 37) Fast Image Texture Classification Using Decision Trees; 38) Constraint Embedding Technique for Multibody System Dynamics; 39) Improved Systematic Pointing Error Model for the DSN Antennas; 40) Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks; 41) More-Accurate Model of Flows in Rocket Injectors; 42) In-Orbit Instrument-Pointing Calibration Using the Moon as a Target; 43) Reliability of Ceramic Column Grid Array Interconnect Packages Under Extreme Temperatures; 44) Six Degrees-of-Freedom Ascent Control for Small-Body Touch and Go; and 45) Optical-Path-Difference Linear Mechanism for the Panchromatic Fourier Transform Spectrometer.

Source record↗

DOE OSTI2022

Thick Electrode Design for Facile Electron and Ion Transport: Architectures, Advanced Characterization, and Modeling

The demand for lithium ion batteries continues to expand for powering applications such as portable electronics, grid-scale energy storage, and electric vehicles. As the application requirements advance, the innovation of lithium ion batteries toward higher energy density and power output is required. Along with the investigation of new materials, an important strategy for increasing battery energy content is to design electrodes with high areal loading to minimize the fraction of nonactive materials such as current collectors, separators, and packaging components, resulting in significant gains in energy content and the reduction of the system-level cost. However, the adoption of thick high areal loading electrodes has been impeded by sluggish charge transport and mechanical instability. With conventional slurry cast electrodes, battery function significantly deteriorates with increases in electrode thickness due to high cell polarization and the incomplete utilization of active materials. Thus, a consideration of approaches that facilitate an understanding and eventual adoption of high-loading electrodes is warranted to enable the deliberate advancement of next-generation batteries. Herein, this Account considers three aspects critical to the science and technology of thick high-loading electrodes. The first discussion covers recent approaches to the design and fabrication of high-loading electrodes. Ensuring electrical contact throughout the electrode is accomplished through the manipulation of conductive additives or using a conductive scaffold within the electrode. Ion transport can be facilitated through electrode design and fabrication approaches that deliberately control the electrode porosity and tortuosity. Second, advanced characterization methodologies are presented as the ability to determine the origins of transport limitations provide the insight needed to deliberately approach future designs. Spectroscopic and diffraction methods have been used to characterize the 2D and 3D pore structure and composition of the electrodes. Furthermore, operando methods that yield spatially and temporally resolved information regarding the progression of the electrochemical reaction are highlighted. The third aspect considered is the utilization of modeling. Physically based continuum models linked with the results of experimental characterization have been demonstrated and then allow the rapid simulation of a variety of deliberate electrode designs and their impacts on functional electrochemistry. Variables relevant to the designs can be tested by the model under a series of use conditions to identify those of most promise for a specific application. Finally, an outlook on future opportunities for high-loading battery electrode research is provided to inform and entice practitioners in the field to pursue these important directions of inquiry.

25 ENERGY STORAGE↗

DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗

DOE OSTI2015

Fully iterative scatter corrected digital breast tomosynthesis using GPU-based fast Monte Carlo simulation and composition ratio update

Purpose: In digital breast tomosynthesis (DBT), scatter correction is highly desirable, as it improves image quality at low doses. Because the DBT detector panel is typically stationary during the source rotation, antiscatter grids are not generally compatible with DBT; thus, a software-based scatter correction is required. This work proposes a fully iterative scatter correction method that uses a novel fast Monte Carlo simulation (MCS) with a tissue-composition ratio estimation technique for DBT imaging. Methods: To apply MCS to scatter estimation, the material composition in each voxel should be known. To overcome the lack of prior accurate knowledge of tissue composition for DBT, a tissue-composition ratio is estimated based on the observation that the breast tissues are principally composed of adipose and glandular tissues. Using this approximation, the composition ratio can be estimated from the reconstructed attenuation coefficients, and the scatter distribution can then be estimated by MCS using the composition ratio. The scatter estimation and image reconstruction procedures can be performed iteratively until an acceptable accuracy is achieved. For practical use, (i) the authors have implemented a fast MCS using a graphics processing unit (GPU), (ii) the MCS is simplified to transport only x-rays in the energy range of 10–50 keV, modeling Rayleigh and Compton scattering and the photoelectric effect using the tissue-composition ratio of adipose and glandular tissues, and (iii) downsampling is used because the scatter distribution varies rather smoothly. Results: The authors have demonstrated that the proposed method can accurately estimate the scatter distribution, and that the contrast-to-noise ratio of the final reconstructed image is significantly improved. The authors validated the performance of the MCS by changing the tissue thickness, composition ratio, and x-ray energy. The authors confirmed that the tissue-composition ratio estimation was quite accurate under a variety of conditions. Our GPU-based fast MCS implementation took approximately 3 s to generate each angular projection for a 6 cm thick breast, which is believed to make this process acceptable for clinical applications. In addition, the clinical preferences of three radiologists were evaluated; the preference for the proposed method compared to the preference for the convolution-based method was statistically meaningful (p < 0.05, McNemar test). Conclusions: The proposed fully iterative scatter correction method and the GPU-based fast MCS using tissue-composition ratio estimation successfully improved the image quality within a reasonable computational time, which may potentially increase the clinical utility of DBT.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Numerical Study of Vortex Core in Turbulent Swirling Flow

Swirling flow has attracted a great deal of research interest. An application of swirling flow can be found in many engineering areas; such as gas cyclone and heat transfer enhancement. In nuclear engineering, swirling flow is used as heat transfer enhancement in rod bundle by installing grid spacer and mixing vane. Because of its extensive applications, there has been a great number of experimental and numerical investigations of swirling flow. McClusky et al., experimentally studied mapping of the lateral velocity in rod bundle using Laser Doppler Velocimetry (LDV) and Conner et al. predicted the flow field in rod bundle by CFD. To understand the characteristics of swirling flow such as tangential velocity and vortex core swirling flow, some researchers investigated swirling flow in a pipe instead of in rod bundle. This is due to the similarity between swirling flow in a pipe and rod bundle. In recent study, flow mapping of 2D velocity field of single-phase swirling flow had been conducted experimentally by Particle Image Velocimetry (PIV). However, this method needs calibration and complicated. The new measurement technique for flow mapping in swirling flow was done by phased array Ultrasonic Velocity Profiler (PAUVP). Hamdani et al., successfully plotted velocity field and determined the swirling core region (vortex core). This technique is nonintrusive and no need calibration. However, there was a limitation of the measurement because phased array UVP cannot capture 360 deg. velocity field inside the pipe. In addition, tangential velocity in CFD results was not agree well with experiment data due asymmetric problem. To evaluate this limitation of experimental method, the numerical simulations using Computational Fluid Dynamics (CFD) needs to be performed. The objective of this study is to use the commercial CFD code Fluent to model the developing turbulent swirling flow induced by twisted tape inside a straight pipe, and to gain an understanding of such a flow and the turbulence models involved in the simulation. Straight pipe is used instead of rod bundle for a simplification. In this work, the Renormalized Group (RNG) k-{epsilon} turbulence model and Reynolds Stress Model (RSM) turbulence model, are performed to study the phenomena of vortex core in a pipe with twisted tape inserts. Effects of the Reynolds number and asymmetric of twisted tape (inclination, {alpha} = 5 deg. and {alpha} = 20 deg.) on swirling core position are examined using water as the testing fluid. Furthermore, the Navier-Stokes equations will be solved to obtain results that can be validated with the previous experimental results. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Higher Efficiency, Demand Flexible Refrigerator with On-Demand Micro-Vibrational De-icing Technology

Refrigerator technology has advanced significantly over the last couple of decades. Today’s refrigerators use only about 25% of the energy that was required to power models built in 1975. Even as they continually improve efficiency to meet standards, refrigerators have increased in size by almost 20%, added energy-consuming features such as through-the-door ice, and provide more benefits than ever before. However, a few challenges and technology gaps are preventing further improvement of the demand responsiveness and efficiency of the refrigerators. One of the major technology gaps in existing refrigerators is their outdated de-icing process. When the evaporator generates frost, an old-fashioned resistive heating element melts the ice. Most refrigerators have a timed defrost cycle, rather than an active system that could monitor the state of the frost. In these systems, not only is the precious electricity used at its least efficient form of conversion (direct conversion of electricity to heat), but also all the latent heat associated with the ice is wasted during the melting process. On top of that, the refrigerator needs to work harder to pull the temperature down after defrosting, and, last but not least, the food quality is severely impacted by the temperature swings during the defrost cycle. According to a study, the EU alone wastes 89 million tons of food in the supply chain every year. Any temperature swing during defrosting (about 6F according to Emerson for low-temperature cases) can negatively impact the shelf life of meat and other products for multiple days. All these issues can happen during the peak demand time of the electric grid. Unlike the conventional systems, the proposed novel advanced micro-vibrational deicing process uses no heat for defrosting. Instead, it uses the micro vibrations generated by a piezoelectric or vibration-generating module to mechanically break ice from the heat exchanger almost instantaneously. The project titled “Higher Efficiency, Demand Flexible Refrigerator with On-Demand Micro-Vibrational De-icing Technology, performed by Ultrasonic Technology Solutions, LLC (UTS) of Knoxville, TN, in collaboration with Emerson (now Copeland), represents the final phase of a multi-year effort funded under the U.S. Department of Energy’s Building Technologies Office (BTO) BENEFIT FOA 2020. Initiated on October 1, 2021, and completed after a nine-month no-cost extension ending September 30, 2025, this project aimed to develop and validate a novel micro-vibrational mechanical defrosting system, achieving more than 25% improvement in defrosting energy efficiency over conventional baseline defrosting technologies. Over sixteen quarters, the project advanced from fundamental ice-mechanical characterization and prototype development to full-scale system integration and validation. Initial efforts established project management infrastructure and characterized ice adhesion properties, followed by the design and fabrication of early aluminum-based prototypes for resonance frequency testing. Subsequent quarters saw rapid technical progression, including the identification of optimal piezoelectric and motor-based vibration mechanisms, the demonstration of effective de-icing over 6x6-inch aluminum surfaces. The team achieved its Go/No-Go milestone by exceeding the 25% energy-efficiency improvement target—reaching up to 3,340% under optimized conditions—and later confirmed that motor-driven systems offered superior performance and energy efficiency compared to piezoelectric alternatives. Continued refinement led to the development of amplifier systems on printed circuit boards, improved control and instrumentation hardware, and integration into full-scale heat exchanger (HX) prototypes at both UTS and Copeland facilities. Multiple vibration-mounting studies and frost-growth experiments guided mechanical optimization and noise-mitigation strategies, achieving a 17.5 dB reduction in sound pressure level and verifying robust mechanical performance. Advanced analyses, including modal and harmonic simulations, established a quantitative understanding of vibrational behavior and de-icing efficiency across >1000 cm² systems. The final project phase successfully demonstrated scalable integration within reach-in and chest freezer prototypes, confirmed >25% efficiency improvements in large-area systems, and completed a comprehensive business model and scale-up strategy identifying electric defrost systems as the primary beachhead market. The culmination of this DOE-supported effort establishes micro-vibrational defrosting as a viable, high-efficiency, low-noise, and demand-flexible de-icing technology, paving the way for commercial deployment and broader application in next-generation refrigeration systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

WISP: Watching grid Infrastructure Stealthily through Proxies (Final Technical Report)

The complex interdependencies of cyber systems (sensors and communications), physical grids and associated electricity market operations make protecting electric power grids a significant challenge. The energy sector is constantly under new, targeted, advanced and dangerous cyber-attacks that have the potential to result in the loss of human life. These threats are further exacerbated by our need to modernize the grid. One focus of cyber security research in smart grids is the securing of the SCADA system through advanced intrusion detection systems (IDS) and bad data detection algorithms in state estimation. These methods either require full knowledge of the system topology and parameters or fail to understand the physical behaviors under attack. WISP (Watching grid Infrastructure Stealthily through Proxies) is designed to provide additional protection to the power grid using only publicly available data. In particular, WISP exploits the spatio-temporal nature of the real time locational marginal prices (LMPs), in conjunction with other information such as bids, weather, outages and load data to analyze anomalous power pricing behaviors and then correlate those observations to localize regions of interest and identify potential cyber events. WISP is non-intrusive as the tool is deployed as a service in the Cloud or on premise and provides reliable information to system operators for enhanced situational awareness, without impeding energy delivery functions. The WISP technology comprises three modules: the data-driven anomaly detection core, the vulnerability and risk analysis and the root cause analysis. The data-driven anomaly detection core performs the tasks of feature selection, anomaly detection and attack region localization. The vulnerability and risk analysis module provides system level information of the vulnerable variables and times, assisting the operators in selecting monitoring and protection nodes. The root cause analysis module takes the detection results and identifies potential operational conditions that contribute to the detected anomalies. In Phase I, we have demonstrated the feasibility and effectiveness of WISP. We developed a realistic electricity market simulator capable of generating normal and attack market data under various operational conditions. We developed a series of cyber-attack detection and analysis algorithms and evaluated them under multiple data sources. Finally, we integrated all modules into an end-to-end software, providing functions for data management, data analytics and visualization. Specifically, we have achieved: (i) real-time data acceptance from external utility interfaces with >99% acceptance rate; (ii) high performance anomaly detection algorithms with >98% detection accuracy and <0.1% false alarm rate; and (iii) ultra-low computing delay <50 milliseconds. Additionally, our team developed algorithms to identify the vulnerable variables in electricity market operations and root cause analysis functions to identify major contributors to the price spikes. These ancillary modules are necessary when deploying WISP in real world industry environment. In Phase II, we have demonstrated the effectiveness of WISP software on realistic largescale power systems. We performed red team testing for the Phase I WISP software and identified software vulnerabilities and implemented corresponding mitigation solutions. We adapted the electricity market simulator for the Texas synthetic 2000-bus system and generated datasets for the false data injection attacks. We created database and visualization interfaces for the Texas system and the ISO New England system. We performed software optimization in terms of operation efficiency, computing speed and detection accuracy. Finally, we tested the software on the Texas system and the ISO New England system and evaluated the detection performance. Overall, we achieved above 89% detection rate, below 3% false alarm rate and below 37 seconds of end-to-end detection delay.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS2012

Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion

Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Leffell, Joshua Isaac↗

NASA NTRS1994

Applications Performance on NAS Intel Paragon XP/S - 15#

The Numerical Aerodynamic Simulation (NAS) Systems Division received an Intel Touchstone Sigma prototype model Paragon XP/S- 15 in February, 1993. The i860 XP microprocessor with an integrated floating point unit and operating in dual -instruction mode gives peak performance of 75 million floating point operations (NIFLOPS) per second for 64 bit floating point arithmetic. It is used in the Paragon XP/S-15 which has been installed at NAS, NASA Ames Research Center. The NAS Paragon has 208 nodes and its peak performance is 15.6 GFLOPS. Here, we will report on early experience using the Paragon XP/S- 15. We have tested its performance using both kernels and applications of interest to NAS. We have measured the performance of BLAS 1, 2 and 3 both assembly-coded and Fortran coded on NAS Paragon XP/S- 15. Furthermore, we have investigated the performance of a single node one-dimensional FFT, a distributed two-dimensional FFT and a distributed three-dimensional FFT Finally, we measured the performance of NAS Parallel Benchmarks (NPB) on the Paragon and compare it with the performance obtained on other highly parallel machines, such as CM-5, CRAY T3D, IBM SP I, etc. In particular, we investigated the following issues, which can strongly affect the performance of the Paragon: a. Impact of the operating system: Intel currently uses as a default an operating system OSF/1 AD from the Open Software Foundation. The paging of Open Software Foundation (OSF) server at 22 MB to make more memory available for the application degrades the performance. We found that when the limit of 26 NIB per node out of 32 MB available is reached, the application is paged out of main memory using virtual memory. When the application starts paging, the performance is considerably reduced. We found that dynamic memory allocation can help applications performance under certain circumstances. b. Impact of data cache on the i860/XP: We measured the performance of the BLAS both assembly coded and Fortran coded. We found that the measured performance of assembly-coded BLAS is much less than what memory bandwidth limitation would predict. The influence of data cache on different sizes of vectors is also investigated using one-dimensional FFTs. c. Impact of processor layout: There are several different ways processors can be laid out within the two-dimensional grid of processors on the Paragon. We have used the FFT example to investigate performance differences based on processors layout.

Saini, Subhash↗

DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

NASA NTRS2015

Exercise Equipment Usability Assessment for a Deep Space Concept Vehicle

With international aspirations to send astronauts to deep space, the world is now faced with the complex problem of keeping astronauts healthy in unexplored hostile environments for durations of time never before attempted by humans. The great physical demands imparted by space exploration compound the problem of astronaut health, as the astronauts must not only be healthy, but physically fit upon destination arrival in order to perform the scientific tasks required of them. Additionally, future deep space exploration necessitates the development of environments conducive to long-duration habitation that would supplement propulsive vehicles. Space Launch System (SLS) core stage barrel sections present large volumes of robust structure that can be recycled and used for long duration habitation. This assessment will focus on one such conceptual craft, referred to as the SLS Derived Habitat (SLS-DH). Marshall Space Flight Center's (MSFC) Advanced Concepts Office (ACO) has formulated a high-level layout of this SLS-DH with parameters such as floor number and orientation, floor designations, grid dimensions, wall placement, etc. Yet to be determined, however, is the layout of the exercise area. Currently the SLS-DH features three floors laid out longitudinally, leaving 2m of height between the floor and ceilings. This short distance between levels introduces challenges for proper placement of exercise equipment such as treadmills and stationary bicycles, as the dynamic envelope for the 95th percentile male astronauts is greater than 2m. This study aims to assess the optimal equipment layout and sizing for the exercise area of this habitat. Figure 1 illustrates the layout of the DSH concept demonstrator located at MSFC. The exercise area is located on the lower level, seen here as the front half of the level occupied by a crew member. This small volume does not allow for numerous or bulky exercise machines, so the conceptual equipment has been limited to a treadmill and stationary bicycle. With the most current treadmill aboard the International Space Station (ISS), the Combined Operational Load-Bearing External Resistance Treadmill (COLBERT), being located in an International Standard Payload Rack (ISPR), the bottom of the conceptual treadmill features a height of 38in. Making the treadmill flush with the floor would be impossible in this rack configuration, as the distance from the outer wall of the spacecraft to the bottom floor would be too shallow. From preliminary sizing, the 38in required for the bottom of the treadmill combined with a 78in operational envelope for a 95th percentile may not be accommodated in the exercise area in a vertical orientation. Figure 2 demonstrates the volume required (in maroon) for an ISPR-bound treadmill in the concept demonstrator. Early indications as seen in this figure indicate that the crew members would contact the ceiling in such an arrangement. An assessment will be conducted to evaluate various orientations of exercise equipment in the concept demonstrator. Orientations to be tested include putting the bottom of the treadmill on the wall, having the treadmill at an angle in the floor both horizontally and vertically, and having a shorter (non-rack bound) treadmill in a vertical orientation on the floor. This assessment will yield findings regarding sizing of the area and how well participants feel they could exercise in such an environment. Due to the restrictions of assessing a microgravity vehicle in a normal-gravity environment, simulations in MSFC's Virtual Environments Lab (VEL) may be necessary. Final deliverables will include recommendations regarding the location and size of possible exercise equipment aboard the SLS-Derived DSH.

Rhodes, Brooke M.↗

DOE OSTI2018

Assessment of CO 2 Storage Resources in Depleted Oil and Gas Fields in the Ship Shoal Area, Gulf of Mexico

The Gulf of Mexico is one of the most important regions in the United States for energy resources and infrastructure. Gulf of Mexico federal offshore oil and gas production accounts for 17% of total U.S. crude oil production and 5% of total U.S. gas production (EIA Gulf of Mexico Fact Sheet). This region presents an excellent combination of high need and significant opportunity for large scale geologic storage of CO 2 . The Ship Shoal Area is located about 20 miles offshore Louisiana within the Gulf Coast federal waters. Miocene and Pliocene sediments in the Ship Shoal Area are proven to provide excellent and secure traps for oil and gas. The Ship Shoal Area contains a large number of depleted oil and gas fields either currently abandoned, or planned for abandonment by 2025, which may provide very significant potential CO 2 storage capacity. GeoMechanics Technologies conducted a comprehensive research project to better characterize Neogene sediments in the Ship Shoal Area for high volume CO 2 storage. The data generated from literature, well data, well log, and formation evaluation were input into Rockwork 16 geologic software to create a detailed geologic model for the Ship Shoal study area. The geologic grids were then fed into TOUGH2, the gas migration model software and FLAC 3D, the geomechanical model software. The fluid flow modeling for both the Miocene and Pliocene simulation from both Block 107 and Block 84 indicated that Pliocene and Miocene are a good reservoir for the CO 2 sequestration. Thirty years of CO 2 injection (at a rate of 1 million metric tons CO 2 per year) and 30 years of observation simulations were run with the fluid flow modeling. The 3D geomechanical modeling for Ship Shoal Block 107 and Block 84 were developed to evaluate induced surface deformation and potential fault reactivation after 30 years CO 2 injection at base of Pliocene and at upper Miocene injection locations. Evaluation results indicate low to no risks for fault slips or fault reactivation. Twelve well bores and 76 well bores for Ship Shoal Blocks 84 and 107 respectively were reviewed for their cement history. Most wells (57 out of 76 and 7 out of 12 wells in Ship Shoal Block 107 and 84 respectively) have good integrity. Nineteen and 5 wells (in Ship Shoal Block 107 and 84 respectively) with no top plug, incomplete cement or Plug and Abandonment information are given yellow cautionary indicators. These cautionary wells may provide leakage paths of CO 2 through the well bores to the USDWs. Using our Quantitative Risk & Decision Analysis Tool (QRDAT) for caprock integrity evaluation, we compared Ship Shoal’s risk to that of In Salah, Sleipner, Kevin Dome, Loudon, Illinois Industrial CCS and Wilmington Graben. We found the risk at the Ship Shoal Blocks 84 and 107 fields are similar to the known CO 2 active sequestration sites, but lower than the Wilmington Graben turbidities offshore California studied site. The risk of natural seismicity in the Gulf of Mexico is relatively low. Geomechanics Technologies has documented the top 25 CO 2 emission sources within the close proximity of the Ship Shoal Block 84 and 107 fields. All the offshore and onshore pipelines have been digitized and can be viewed in an interactive website (http://www.geomechanicstech.com/shipshoal.html). We also performed a feasibility study on the potential for converting existing oil or gas pipelines for CO 2 transport. There are abandoned, idled and retired onshore and offshore pipelines. However, there is no standard specification for maximum pipeline pressures needed for CO 2 transport. It is a function of design, materials, and testing provided; thus it will be the responsibility of the pipeline operator to correctly determine, maintain and operate within the limits of the pipelines. There are a few transit corridors extending from onshore Louisiana (Cailou Bay) to offshore trunk-lines. The cost for constructing a pipeline has increased about 46% from 2015 to an average cost of $5,064,046 per mile in 2016. We use the NETL CO 2 storage resource mass estimate to calculate the potential resources for Ship Shoal Block 84 and 107 fields using the depleted oil and gas reservoir volume and sand volume generated from our geologic modeling. The estimated storage resource results are greater when using the sand volumes derived through geologic modeling than the oil and gas reservoir data. The difference is due to the depleted oil and gas reservoir information not accounting for the water-flooded sand located below the oil/gas-water contact or unproductive sand units. Resource calculation using the sand volume obtained through geologic modeling overestimate the storage capacity as the model accounts for all sand within the formation, not just the interconnected sand.

02 PETROLEUM↗

DOE OSTI2022

Using Energy Storage-Based Grid Forming Inverters for Operational Reserve in Hybrid Diesel Microgrids

In remote arctic communities, where access to a bulk electrical grid interconnection is not available, the implementation of islanded microgrids has been the most viable way to produce and distribute electricity services to their inhabitants. Historically, these islanded grids have relied primarily on diesel generators or hydropower resources to supply the baseload. However, this practice can result in increased expense due to the high costs associated with fuel transportation and the significant amounts of on-site storage necessary when fuel transportation is unavailable during winter months. In order to mitigate this problem, arctic microgrids have started to transition to a hybrid-source operational mode by incorporating renewable energy sources that are inherently variable in nature, such as wind or solar. Due to their highly stochastic behavior, these hybrid-source islanded microgrids can pose potential issues related to power quality due to introduction of rapid net load fluctuations and inability of diesel generators to respond rapidly. In addition, non-firm stochastic sources can require significant idling diesel generator resources to serve as spinning reserves, which is inefficient and wasteful. This work studies the problems that may arise in the transient dynamics of a real-world hybrid diesel microgrid when subjected to a loss of wind generation. Moreover, this work proposes a transition from a diesel spinning reserve to a battery energy-storage system (BESS) operating reserve scheme. The study of the proposed transition is important in establishing the fundamental implication of transient dynamics and the potential benefits of integrating a BESS as a spinning reserve in terms of stability, frequency nadir, and transient voltage deviation. The methods to investigate and validate the transient dynamics relied on both electromagnetic simulation models of GFMIs and a commercially available GFMI in an experimental power hardware-in-the-loop setup. The simulation results showed that the proposed operating reserve scheme improves the power quality of the system in terms of voltage deviation and frequency nadir when the microgrid is subjected to a loss of wind generation scenario. Depending on the simulation cases, adding a GFMI reduced the frequency nadir between 65.3% and 86.7%. Moreover, the reductions in the voltage deviations improved between 3.6% and 23.0%. From these results it can be concluded that the integration of a GFMI can reduce the frequency nadir in a hybrid diesel microgrid, and in turn, reduce diesel consumption, which improves system reliability while reducing fuel expenses. Furthermore, the novelty of this work relies on the fact that the offline simulation results were validated using a power hardware-in-the-loop platform that incorporated a 100 kVA commercially available GFMI as the device under test.

Hernandez-Alvidrez, Javier↗