Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Simulation”

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.

At least 199 records · Page 11

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

Accelerating geostatistical modeling using geostatistics-informed machine Learning

Ordinary Kriging (OK) is a popular geostatistical algorithm for spatial interpolation and estimation. The computational complexity of OK changes quadratically and cubically for memory and speed, respectively, given the number of data. Therefore, it is computationally intensive and also challenging to process a large set of data, especially in three-dimensional (3D) cases. This paper develops a geostatistics-informed machine learning (GIML) model to improve the efficiency of OK by reducing the number of points required to be estimated using OK. Specifically, only a very few of the unknown points are estimated by OK to get the weights and estimations, which are used as the training dataset. Moreover, the governing equations of OK are used to guide our proposed machine learning to better reproduce the spatial distributions. Our results show that the proposed GIML can reduce the computational time of OK by at least one order of magnitude. The effectiveness of the GIML is evaluated and compared using a 2D case. Furthermore, we demonstrate its efficiency and robustness by considering a different number of training samples on various 3D simulation grids.

58 GEOSCIENCES↗

Modeling Analysis of the Swell and Wind-Sea Climate in the Salish Sea

Accurate model hindcast of wave climate in complex estuarine systems is challenging because high-resolution wave model and wind forcing are required. In this study, a modeling approach using the unstructured-grid Simulating WAves Nearshore and wind product from a high-resolution regional Weather Research and Forecasting hindcast was used to simulate the swell and wind-sea climate in the Salish Sea, a large estuary with many interconnected waterways on the Pacific Northwest coast of North America. The model hindcast was validated with observed data at four wave buoys in the Salish Sea. Spatial distribution and seasonal variations of wave climate in the Salish Sea were analyzed. Of the three major basins in the Salish Sea, the Strait of Juan de Fuca has the largest waves and is dominated by swells propagated from the Pacific Ocean. Significant wave heights in the Strait of Georgia have spatial and seasonal distribution patterns similar to those found in the Strait of Juan de Fuca. Waves in Puget Sound are small and primarily dominated by the wind-sea climate. Strong seasonal variations are observed in the Strait of Juan de Fuca and Strait of Georgia, but there is little seasonality of wave climate in Puget Sound. The high-resolution wave hindcast conducted in this study provides a comprehensive and important data set for better understanding the role of wave climate in coastal processes and natural hazards assessment in the Salish Sea.

wave climate, salish sea, modeling, swell, sea sta↗

Wave resource characterization at regional and nearshore scales for the U.S. Alaska coast based on a 32-year high-resolution hindcast

A wave resource characterization was performed for the southern coast of Alaska based on a 32-year hindcast covering the period from 1979 to 2010. The characterization closely followed International Electrotechnical Commission (IEC) Technical Specifications. An unstructured-grid Simulating WAves Nearshore (SWAN) model, which had an approximate spatial resolution of 300 m within 30 km from the nearest shoreline, was developed. Additionally, extensive model validation and error characterization was performed. The model was found to perform well with an average absolute percent error of 8.6% in significant wave height, averaged over 18 buoys. Statistics for the six IEC wave resource parameters were calculated and aggregated at 20 km from shore to quantify the incident wave power and its variability at a regional scale. The southern coast of the Aleutian Archipelago was found to have the most available wave energy in the region. A nearshore resource assessment was performed by evaluating resource hotspots located 1 km from shore. More than 900 nearshore stations had an Optimum Hotspot Identifier value of 5 (kW/m) at diverse water depths, thereby positioning Alaska as a promising location for wave energy development.

54 ENVIRONMENTAL SCIENCES↗

Source Contributions to Fine Particulate Matter and Attributable Mortality in India and the Surrounding Region

Fine particulate matter (PM 2.5 ) is a leading mortality risk factor in India and the surrounding region of South Asia. This study evaluates the contribution of emission sectors and fuels to PM 2.5 mass for 29 states in India and 6 surrounding countries (Pakistan, Bangladesh, Nepal, Bhutan, Sri Lanka, and Myanmar) by combining source-specific emission estimates, stretched grid simulations from a chemical transport model, high resolution hybrid PM 2.5 , and disease-specific mortality estimates. We find that 1.02 (95% Confidence Interval (CI): 0.78-1.26) million deaths in South Asia attributable to ambient PM 2.5 in 2019 were primarily from three leading sectors: residential combustion (28%), industry (15%), and power generation (12%). Solid biofuel is the leading combustible fuel contributing to the PM 2.5 -attributable mortality (31%), followed by coal (17%), and oil and gas (14%). State-level analyses reveal higher residential combustion contributions (35%-39%) in states (Delhi, Uttar-Pradesh, Haryana) with high ambient PM 2.5 (>95 µg/m 3 ). The combined mortality burden associated with residential combustion (ambient) and household air pollution (HAP) in India is 0.72 million (95% CI:0.54-0.89) (68% attributable to HAP, 32% attributable to residential combustion). Our results illustrate the potential to reduce PM 2.5 mass and improve population health by reducing emissions from traditional energy sources across multiple sources in South Asia.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Selective bond formation triggered by short optical pulses: quantum dynamics of a four-center ring closure

Here we report bond formation induced by an ultrashort UV pulse. The photochemical process is described by quantum dynamics as coherent electronic and nuclear motions during the ultrashort pulse induced ring closure of norbornadiene to quadricyclane. Norbornadiene consists of two ethylene moieties connected by a rigid (CH 2 ) 3 bridge. Upon photoexcitation, two new sigma bonds are formed, resulting in the closure of a four-atom ring. As a medium-sized polyatomic molecule, norbornadiene exhibits a high density of strongly coupled electronic states from about 6 eV above the ground state. We report on inducing the formation of the new bonds using a short femtosecond UV pulse to pump a non-equilibrium electronic density in the open form that evolves towards the closed ring form. As the coherent electronic-nuclear coupled dynamics unfold, the excited states change character through non-adiabatic interactions and become valence states for the two new C–C bonds of quadricyclane. Our three-dimensional fully quantum dynamical grid simulations during the first 200 fs show that short UV pulses of different polarization initiate markedly different initial non-equilibrium electronic densities that follow different dynamical paths to the S 0 /S 1 conical intersection. They lead to different initial relative yields of quadricyclane, thereby opening the way to controlling bond-making with attopulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A NEW DENSITY VARIANCE-MACH NUMBER RELATION FOR SUBSONIC AND SUPERSONIC ISOTHERMAL TURBULENCE

The probability density function of the gas density in subsonic and supersonic, isothermal, driven turbulence is analyzed using a systematic set of hydrodynamical grid simulations with resolutions of up to 1024{sup 3} cells. We perform a series of numerical experiments with root-mean-square (rms) Mach number M ranging from the nearly incompressible, subsonic (M=0.1) to the highly compressible, supersonic (M=15) regime. We study the influence of two extreme cases for the driving mechanism by applying a purely solenoidal (divergence-free) and a purely compressive (curl-free) forcing field to drive the turbulence. We find that our measurements fit the linear relation between the rms Mach number and the standard deviation (std. dev.) of the density distribution in a wide range of Mach numbers, where the proportionality constant depends on the type of forcing. In addition, we propose a new linear relation between the std. dev. of the density distribution {sigma}{sub {rho}} and that of the velocity in compressible modes, i.e., the compressible component of the rms Mach number, M{sub comp}. In this relation the influence of the forcing is significantly reduced, suggesting a linear relation between {sigma}{sub {rho}} and M{sub comp}, independent of the forcing, and ranging from the subsonic to the supersonic regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Divertor leg filaments in NSTX-U

This study presents the first observation of divertor localized turbulence in the NSTX-Upgrade spherical tokamak. Previous work in NSTX discharges described fluctuations on the divertor target due to blobs generated on the low field side midplane. In the region that shows disconnection from upstream turbulence (private flux region and proximity of the strike point), divertor-localized fluctuations are observed via imaging of C III and D-α emission in diverted L-mode discharges and are characterized in this paper. Field-aligned filaments connected to the divertor target plate are radially localized at the separatrix on the outer divertor leg and in the private flux region on the inner divertor leg. These are limited to the region below the X-point with fluctuation levels up to 10%–20%. The filaments have comparable poloidal and radial correlation lengths (10–100 ion gyroradii ρ i ) and parallel correlation lengths of several meters. Toroidal mode numbers are in the range of 10–30 and 2–10 for outer and inner leg filaments. The poloidal motion of outer leg filaments has the same direction and comparable magnitude to advection due to E x B drift calculated by the multi-fluid edge transport code UEDGE with inclusion of cross field drifts. Disconnection between inner and outer leg filaments as well as the absence of correlation with upstream turbulence support the hypothesis of X-point disconnection for divertor leg filaments. Simulations with the ArbiTER linear eigenvalue solver and a resistive-balooning model were performed with a simulation grid limited to the divertor legs. Instabilities localized to the bad curvature side of the divertor legs were observed, in qualitative agreement with the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Joint Matrix Completion and Compressed Sensing for State Estimation in Low-observable Distribution System

Limited measurement availability at the distribution grid presents challenges for state estimation and situational awareness. This paper combines the advantages of two sparsity-based state estimation approaches (matrix completion and compressive sensing) that have been proposed recently to address the challenge of unobservability. The proposed approach exploits both the low rank structure and a suitable transform domain representation to leverage the correlation structure of the spatio-temporal data matrix while incorporating the powerflow constraints of the distribution grid. Simulations are carried out on three phase unbalanced IEEE 37 test system to verify the effectiveness of the proposed approach. The performance results reveal - (1) the superiority over traditional matrix completion and (2) very low state estimation errors for high compression ratios representing very low observability.

Dahale, Shweta↗

Integrated Transmission and Distribution Systems Restoration with Distributed Generation Scheduling

The self-healing smart grids enable system operators to quickly and efficiently recover transmission and distribution systems from outages and blackouts. The current restoration procedure solves the transmission and distribution system restoration problems separately without considering their impacts on each other, which may prolong the restoration process. However, with the emerging active distribution networks, distributed generators can contribute into the bottom-up restoration strategy. The interaction between the restoration of two interconnected systems should be taken into account. In this paper, a hybrid model of transmission and distribution systems restoration (HTDSR) is proposed by exchanging data of power and voltage between two separate optimization problems. The effectiveness of the HTDSR algorithm is demonstrated through case studies of the integrated IEEE 13-bus distribution network and IEEE 39-bus transmission grid. Simulation results demonstrate that the integrated restoration strategy with optimal scheduling of distributed generation sources can expedite the restoration process.

Distributed generation↗

Using a critical infrastructure game to provide realistic observation of the human in the loop by criminal justice students

Understanding human behavior is crucial in anticipating adversarial actions during cyberattacks. The Criminal Justice (CJ) discipline offers the necessary frameworks to unpack the complex facets of adversarial behavior and movement, and should therefore be leveraged for their possible contributions to the area of proactive cybersecurity. Yet the discipline remains weak at training current and future CJ workforce on these matters in a hands-on manner. This paper presents a cybersecurity training exercise where a power grid simulator is used to educate CJ students via experiential learning about concepts of cyberattacks and cybersecurity as well as exposing them to doing hands-on cybersecurity field research. The paper reports on Game use as an important opportunity to observe humans put under additional stress in operating conditions. The paper discusses what CJ students learn from multidisciplinary simulation-based exercises, the challenges and limitations they face, and how training this workforce could help contribute towards proactive cyberdefense of critical infrastructure.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cyber Framework for Steering and Measurements Collection Over Instrument-Computing Ecosystems

We propose a framework to develop cyber solutions to support the remote steering of science instruments and measurements collection over instrument-computing ecosystems. It is based on provisioning separate data and control connections at the network level, and developing software modules consisting of Python wrappers for instrument commands and Pyro server-client codes that make them available across the ecosystem network. We demonstrate automated measurement transfers and remote steering operations in a microscopy use case for materials research over an ecosystem of Nion microscopes and computing platforms connected over site networks. The proposed framework is currently under further refinement and being adopted to science workflows with automated remote experiments steering for autonomous chemistry laboratories and smart energy grid simulations.

Al Najjar, Anees↗

Sequence Impedance Measurement of Utility-Scale Wind Turbines and Inverters – Reference Frame, Frequency Coupling, and MIMO/SISO Forms

Sequence impedance responses with or without considering frequency coupling in both MIMO and SISO forms are increasingly used for the stability analysis of three-phase power electronic systems; however, many aspects of sequence impedance measurement are not fully explored. It is not clear if the sequence impedance has a reference frame similar to the dq impedance. If so, the role of the grid voltage angle estimation in aligning the sequence impedance reference frame has not been discussed. Additionally, existing methods for measuring the sequence impedance with frequency coupling are complicated, are not feasible for large wind turbines and inverters, and provide the sequence impedance responses in only either MIMO or SISO form. This paper presents a sequence impedance measurement method that considers the frequency coupling, performs reference frame alignment, demonstrates the impact of the grid voltage angle estimation, and obtains the sequence impedance response in both MIMO and SISO forms. This paper demonstrates the proposed method and practical problems associated with the sequence impedance measurement of utility-scale wind turbines and inverters on a 1.9-MW Type III wind turbine and a 2.2-MVA inverter using an impedance measurement system built around a 7-MW/13.8-kV grid simulator and a 5-MW dynamometer.

17 WIND ENERGY↗

smoothG

smoothG is a collection of parallel C++ classes/functions that algebraically constructs reduced models of different resolutions from a given high-fidelity graph model. In addition, smoothG also provides efficient linear solvers for the reduced models. Other than pure graph problem, the software finds its application in subsurface flow and power grid simulations in which graph Laplacians are found

Barker, AndrewT.↗

Fundamental Drivers of the Cost and Price of Operating Reserves

Operating reserves impose a cost on the electric power system by forcing system operators to keep partially loaded spinning generators available for responding to system contingencies variable demand. In many regions of the United States, thermal power plants provide a large fraction of the operating reserve requirement. Alternative sources of operating reserves, such as demand response and energy storage, may provide more efficient sources of these reserves. However, to estimate the potential value of these services, the cost of reserve services under various grid conditions must first be established. This analysis used a commercial grid simulation tool to evaluate the cost and price of several operating reserve services, including spinning contingency reserves and upward regulation reserves. These reserve products were evaluated in a utility system in the western United States, considering different system flexibilities, renewable energy penetration, and other sensitivities. The analysis demonstrates that the price of operating reserves depend highly on many assumptions regarding the operational flexibility of the generation fleet, including ramp rates and the fraction of fleet available to provide reserves.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cost-Benefit Analysis of Flexibility Retrofits for Coal and Gas-Fueled Power Plants: August 2012 - December 2013

High penetrations of wind and solar power plants can induce on/off cycling and ramping of fossil-fueled generators. This can lead to wear-and-tear costs and changes in emissions for fossil-fueled generators. Phase 2 of the Western Wind and Solar Integration Study (WWSIS-2) determined these costs and emissions and simulated grid operations to investigate the full impact of wind and solar on the fossil-fueled fleet. This report studies the costs and benefits of retrofitting existing units for improved operational flexibility (i.e., capability to turndown lower, start and stop faster, and ramp faster between load set-points).

02 PETROLEUM↗

A Community-Based Approach to Leading the Nation in Smart Energy Use

Project Objectives The AEP Ohio gridSMART® Demonstration Project (Project) achieved the following objectives: • Built a secure, interoperable, and integrated smart grid infrastructure in northeast central Ohio that demonstrated the ability to maximize distribution system efficiency and reliability and consumer use of demand response programs that reduced energy consumption, peak demand, and fossil fuel emissions. • Actively attracted, educated, enlisted, and retained consumers in innovative business models that provided tools and information reducing consumption and peak demand. • Provided the U.S. Department of Energy (DOE) information to evaluate technologies and preferred smart grid business models to be extended nationally. Project Description Ohio Power Company (the surviving company of a merger with Columbus Southern Power Company), doing business as AEP Ohio (AEP Ohio), took a community-based approach and incorporated a full suite of advanced smart grid technologies for 110,000 consumers in an area selected for its concentration and diversity of distribution infrastructure and consumers. It was organized and aligned around: • Technology, implementation, and operations • Consumer and stakeholder acceptance • Data management and benefit assessment Combined, these functional areas served as the foundation of the Project to integrate commercially available products, innovative technologies, and new consumer products and services within a secure two-way communication network between the utility and consumers. The Project included Advanced Metering Infrastructure (AMI), Distribution Management System (DMS), Distribution Automation Circuit Reconfiguration (DACR), Volt VAR Optimization (VVO), and Consumer Programs (CP). These technologies were combined with two-way consumer communication and information sharing, demand response, dynamic pricing, and consumer products, such as plug-in electric vehicles and smart appliances. In addition, the Project incorporated comprehensive cyber security capabilities, interoperability, and a data assessment that, with grid simulation capabilities, made the demonstration results an adaptable, integrated solution for AEP Ohio and the nation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Fundamental Drivers of the Cost and Price of Operating Reserves

Operating reserves impose a cost on the electric power system by forcing system operators to keep partially loaded spinning generators available for responding to system contingencies variable demand. In many regions of the United States, thermal power plants provide a large fraction of the operating reserve requirement. Alternative sources of operating reserves, such as demand response and energy storage, may provide more efficient sources of these reserves. However, to estimate the potential value of these services, the cost of reserve services under various grid conditions must first be established. This analysis used a commercial grid simulation tool to evaluate the cost and price of several operating reserve services, including spinning contingency reserves and upward regulation reserves. These reserve products were evaluated in a utility system in the western United States, considering different system flexibilities, renewable energy penetration, and other sensitivities. The analysis demonstrates that the price of operating reserves depend highly on many assumptions regarding the operational flexibility of the generation fleet, including ramp rates and the fraction of fleet available to provide reserves.

24 POWER TRANSMISSION AND DISTRIBUTION↗