Search NASA⌕ Search

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.

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

At least 685 records · Page 38

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

Measurements and Model Improvement: Insight into NWP Model Error Using Doppler Lidar and Other WFIP2 Measurement Systems

Abstract Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s −1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement. Significance Statement Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.

Meteorology & Atmospheric Sciences↗

Evaluation of Runtime Performance for xRage Equations of State Input Settings and Cropped EOS Tables

Using the xRage hydro code with its intrinsic equations of state (EOS) settings does not produce the most accurate representation of the problem due to xRage's internal alterations to the fundamental EOS data. Utilizing EOSPAC with default settings in the hydro code uses the unaltered EOS SESAME data to produce what is considered the most accurate results. However, this setup for the simulation is very computationally inefficient due to xRage's reliance on pressure and temperature as its look-up variables in combination with EOSPAC's use of a rational function interpolator. This report explores the results and performance of running a simulation with the option, within EOSPAC, to invert the table at setup, using a linear or rational interpolator, and using an EOS table cut down to only the region actually sampled in the problem. Air and Deuterium are used as the test materials for the simulations. While using EOS data directly with EOSPAC, the rational interpolator takes about twice as long as the linear one. Additionally, using a windowed EOS table does not show a large increase in performance. An EOS table that has had more energy grid points interpolated and added in by EOSPAC showed more accurate results, suggesting that a more dense EOS table would be more accurate.

97 MATHEMATICS AND COMPUTING↗

Comparisons of Two Spatial Implementations of a Crop Model Using Remotely Sensed Observations over Southeastern United States

Global food security is one of the most pressing issues of the current century, particularly for developing nations. Agricultural simulation models can be a key component in testing new technologies, seeds and cultivars etc. However, inaccurate input information, model related errors and the mode of implementation can also add to model uncertainties. In this study, the crop model is implemented in two separate fashions: a)gridded (GriDSSAT model) and b) using random spatial ensembles (RHEAS model). This is done in the Southeastern US to evaluate and understand the modelperformance over a region data availabilities. Once the model performance is assessed, multiple satellite based earth observation parameters such as soil moisture, vegetation index etc. can be assimilated into crop models to reduce input and model related uncertainties particularly in data limited regions. In this study, the National Agricultural Statistical Services (NASS) reported yield data at county levels are used for comparison andvalidation purposes. The GriDSSAT model estimation of corn yields in comparison with the reported NASS yields showed an overall RMSD of nearly 3720 (kg/ha) whereas RMSD for the RHEAS model implementation was 3550 (kg/ha). Overall the GriDSSAT model had negative bias of nearly 2400 kg/ha (except for 2013) while RHEAS had a slight positive bias of 400 kg/ha (approx.).

SERVIR↗

CFD Predictions of Fluctuating Pressure Environments on NASA’s SLS Ascent Unsteady Aerodynamics Wind Tunnel Test

This paper details work that is being conducted on predicting fluctuating pressure environments on NASA launch vehicles using computational fluid dynamics (CFD). The accurate characterization of these aeroacoustics environments is necessary in the analysis of vehicle structural health and aerodynamic performance and in the vehicle design process. High-fidelity simulations of the unsteady flow over NASA’s SLS Ascent Unsteady Aerodynamics wind tunnel Test (AUAT) using the hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulations (LES) methodologies in Mississippi States’ Loci/CHEM solver are presented. The two available methods, Nelson-Nichols and dynamic hybrid RANS-LES (DHRL), are compared head-to-head on identical grids at two Mach numbers. At the high-subsonic Mach number, locally supersonic flow expanding over a shoulder induces a separation-reattachment system that is predicted poorly by the Nelson-Nichols method due to an under-prediction in unsteady content. The DHRL solution predicts the spatial extent and frequency domain response of this system well but under-predicts the peak noise levels due to a delay in the RANS-to-LES transition. At the low-supersonic Mach number, compression corner dynamics dominate the unsteady flow. While the Nelson Nichols solution captures this phenomenon well, another solution on a refined grid indicates that some model stress depletion is occurring. The DHRL solution again performs well at this Mach number and matches the wind tunnel data fairly closely, but perhaps most impressively, is fairly agnostic to increasing grid resolution. The findings in this paper are expected to be generally applicable and will be used to guide the prediction of aeroacoustics environments of other NASA launch vehicles.

CFD↗

Common Envelope Light Curves. I. Grid-code Module Calibration

The common envelope (CE) binary interaction occurs when a star transfers mass onto a companion that cannot fully accrete it. The interaction can lead to a merger of the two objects or to a close binary. The CE interaction is the gateway of all evolved compact binaries, all stellar mergers, and likely many of the stellar transients witnessed to date. CE simulations are needed to understand this interaction and to interpret stars and binaries thought to be the byproduct of this stage. At this time, simulations are unable to reproduce the few observational data available and several ideas have been put forward to address their shortcomings. The need for more definitive simulation validation is pressing and is already being fulfilled by observations from time-domain surveys. In this article, we present an initial method and its implementation for post-processing grid-based CE simulations to produce the light curve so as to compare simulations with upcoming observations. Here we implemented a zeroth order method to calculate the light emitted from CE hydrodynamic simulations carried out with the 3D hydrodynamic code Enzo used in unigrid mode. The code implements an approach for the computation of luminosity in both optically thick and optically thin regimes and is tested using the first 135 days of the CE simulation of Passy et al., where a 0.8 M {sub ⊙} red giant branch star interacts with a 0.6 M {sub ⊙} companion. This code is used to highlight two large obstacles that need to be overcome before realistic light curves can be calculated. We explain the nature of these problems and the attempted solutions and approximations in full detail to enable the next step to be identified and implemented. We also discuss our simulation in relation to recent data of transients identified as CE interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

FLISR in the Presence of DERs ADMS Test Bed Use Case

Fault location, isolation, and restoration (FLISR) is one of the distribution automation applications that utilities are most interested in. It operates groups of switches on distribution feeders to improve the reliability of power delivery after localizing outages. FLISR is an essential function for enabling a self-healing grid, and it directly affects grid reliability and resilience. The presence of DERs brings both new opportunities and challenges for FLISR applications. For instance, some DERs can work as backup generation sources and help energize local networks so more switching options can be achieved with the coordination of FLISR operations; however, the existence of DERs also brings new challenges. The bidirectional power flow resulting from DERs might require dynamic changes to protective devices' fault current settings, and thus the FLISR application will need to dynamically adjust its operations in coordination with the dynamic protection settings. Also, some FLISR applications make their service restoration decisions and close tie-switches based on predefined rules; however, the presence of DERs and their impact on altering net loads on the local feeder network are not considered. For this use case, the project team is using the FLISR application of the Survalent ADMS operating on a real-time simulation of a feeder from Central Georgia EMC, a cooperative electric utility. The performance of the FLISR application is being evaluated for different locations and combinations of DERs, including grid-forming battery energy storage systems.

ADMS↗

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗

Survivability of Autonomous Microgrid during Overload Events

Grid-forming sources are voltage sources that draw necessary currents to meet any load changes. A load step can cause part or all of these sources to become overloaded in a microgrid. This paper presents an overload mitigation controller that addresses the two overload issues in a microgrid by actively controlling the sources’ frequency. When part of the sources in a microgrid is overloaded, the controller autonomously transfers the extra load to other sources by rapidly reducing its frequency. The frequency difference between sources during transient results in a change of phase angle, which redistributes the power flow. When all sources in a microgrid are overloaded, each source keeps dropping the frequency. Therefore, under frequency load shedding can be used to trip the non-critical loads resulting in the survival of microgrid. The advantages of these concepts are that communications between sources are not needed during transient, and the robust voltage control is maintained. Lastly, simulation and field tests from CERTS/AEP microgrid test site verify that the control strategy is effective in both purely inverter-based microgrids and inverter & generator mixed microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Fabrication, Modeling, and Testing of a Prototype Thermal Energy Storage Containment

Increasing penetration of variable renewable energy resources requires the deployment of energy storage at a range of durations. Long-duration energy storage (LDES) technologies will fulfill the need to firm variable renewable energy resource output year round; lithium-ion batteries are uneconomical at these durations. Thermal energy storage (TES) is one promising technology for LDES applications because of its siting flexibility and ease of scaling. Particle-based TES systems use low-cost solid particles that have higher temperature limits than the molten salts used in traditional concentrated solar power systems. A key component in particle-based TES systems is the containment silo for the high-temperature (>1100 degrees C) particles. This study combined experimental testing and computational modeling methods to design and characterize the performance of a particle containment silo for LDES applications. A laboratory-scale silo prototype was built and validated the congruent transient finite element analysis (FEA) model. The performance of a commercial-scale silo was then characterized using the validated model. The commercial-scale model predicted a storage efficiency above 95% after 5 days of storage with a design storage temperature of 1200 degrees C. Insulation material and concrete temperature limits were considered as well. The validation of the methodology means the FEA model can simulate a range of scenarios for future applications. This work supports the development of a promising LDES technology with implications for grid-scale electrical energy storage, but also for thermal energy storage for industrial process heating applications.

clean energy↗

Black-Starting Microgrids Using Inverter-Based Distributed Energy Resources

Recent trends in power systems show increasing deployment of inverter-based distributed energy resources (DERs) at the distribution level. Efforts have also been made to use inverter-based DERs for system restoration of critical sections of distribution feeders. In the past few years, there has been active consideration of using local grid-forming (GFM) DERs, which can form their own voltage and frequency, to black-start microgrids instead of waiting for the bulk system for supply restoration. However, using current-limited inverter-based DERs for black-start operations with high inrush currents can pose some technical challenges. This paper presents detailed electromagnetic transient studies on a digital real-time simulator to assess the operational feasibility of using local inverter-based GFM DERs for microgrid black-start operations as a bottom-up approach instead of a top-down approach where the substation is used to restore critical sections of the feeder. The paper presents key modeling details for distribution system components for studying GFM DER-based microgrid black-start operations. The paper presents simulation results comparing different black-start approaches using GFM DERs for two test systems: i) a simple test microgrid and ii) a real-world utility microgrid currently being designed and deployed by Holy Cross Energy in Colorado.

black start↗

Mapping unstructured grid computations to massively parallel computers

Investigated here is this mapping problem: assign the tasks of a parallel program to the processors of a parallel computer such that the execution time is minimized. First, a taxonomy of objective functions and heuristics used to solve the mapping problem is presented. Next, we develop a highly parallel heuristic mapping algorithm, called Cyclic Pairwise Exchange (CPE), and discuss its place in the taxonomy. CPE uses local pairwise exchanges of processor assignments to iteratively improve an initial mapping. A variety of initial mapping schemes are tested and recursive spectral bipartitioning (RSB) followed by CPE is shown to result in the best mappings. For the test cases studied here, problems arising in computational fluid dynamics and structural mechanics on unstructured triangular and tetrahedral meshes, RSB and CPE outperform methods based on simulated annealing. Much less time is required to do the mapping and the results obtained are better. Compared with random and naive mappings, RSB and CPE reduce the communication time two fold for the test problems used. Finally, we use CPE in two applications on a CM-2. The first application is a data parallel mesh-vertex upwind finite volume scheme for solving the Euler equations on 2-D triangular unstructured meshes. CPE is used to map grid points to processors. The performance of this code is compared with a similar code on a Cray-YMP and an Intel iPSC/860. The second application is parallel sparse matrix-vector multiplication used in the iterative solution of large sparse linear systems of equations. We map rows of the matrix to processors and use an inner-product based matrix-vector multiplication. We demonstrate that this method is an order of magnitude faster than methods based on scan operations for our test cases.

Hammond, Steven Warren↗

Microgrid Energy Management System Integration with Advanced Distribution Management System

The Integrated Distribution Management System (IDMS) project was initiated to demonstrate the interactive operation of microgrid systems and the distribution systems with which they interconnect. The key technologies for this are the microgrid management system and the utility distribution management system. The IDMS project successfully demonstrated that a utility’s advanced distribution management system/distributed energy resources management system (ADMS/DERMS) could effectively manage microgrids to provide visibility and control functionalities as well as use the microgrid as a dispatchable resource to support the utility grid. The DERMS accomplishes this by determining active and reactive power needs at the point of common coupling (PCC) using advanced applications like volt/VAR watt optimization (VVWO) and load relief (LR) with the underlying core applications state estimation (SE) and load flow (LF). The ADMS/DERMS can use the microgrid as a resource to resolve and prevent violations in the grid and to optimize the operational working state of the grid. The IDMS project demonstrated that a utility-operated ADMS with embedded DERMS functionality can flexibly manage a variety of microgrids and other aggregated distributed energy resources (DER) in concert with the wider distribution grid. Microgrids can provide grid services in any number of different ways to meet the operational needs of the distribution utility. The manner of aggregation — microgrid or virtual power plant — is not necessarily relevant to the utility as long as the grid services from the aggregated DER are available and can be managed by its ADMS/DERMS for the stability and reliability of the grid. The IDMS project demonstrated this integrated ADMS concept by combining hardware/software-in-the-loop testing with commercial products from different vendors and a utility’s network model. The project integrated the Schneider Electric EcoStruxure™ ADMS with DERMS with the Schweitzer Engineering Laboratories (SEL) POWERMAX ® Microgrid Control System, and the simulated resources with energy company PECO’s model of a utility-owned microgrid, establishing an operational relationship in which the utility manages the operational functionality of a microgrid at the PCC and that provides the utility with the capability to control the comprehensive power system, inclusive of the macrogrid and microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A simple atmospheric model on the sphere with 100% parallelism

A simple shallow water equations model used to test numerical weather prediction schemes is presented. The spherical geometry is transformed into a formally doubly periodic configuration which can be implemented on a massively parallel computer at 100% parallelism. A nonstaggered grid with all variables defined at the grid points (circles) is used. The parallel structure was simulated by performing step 1 in a FORTRAN array, on an Amdahl serial machine, changing the signs of the winds on the dark side.

Kalnay, E.↗

Parallel Implementation of the Discontinuous Galerkin Method

This paper describes a parallel implementation of the discontinuous Galerkin method. Discontinuous Galerkin is a spatially compact method that retains its accuracy and robustness on non-smooth unstructured grids and is well suited for time dependent simulations. Several parallelization approaches are studied and evaluated. The most natural and symmetric of the approaches has been implemented in all object-oriented code used to simulate aeroacoustic scattering. The parallel implementation is MPI-based and has been tested on various parallel platforms such as the SGI Origin, IBM SP2, and clusters of SGI and Sun workstations. The scalability results presented for the SGI Origin show slightly superlinear speedup on a fixed-size problem due to cache effects.

Baggag, Abdalkader↗

Reinforcement Learning Applied to Cognitive Space Communications

The future of space exploration depends on robust, reliable communication systems. As the number of such communication systems increase, automation is fast becoming a requirement to achieve this goal. A reinforcement learning solution can be employed as a possible automation method for such systems. The goal of this study is to build a reinforcement learning algorithm which optimizes data throughput of a single actor. A training environment was created to simulate a link within the NASA Space Communication and Navigation (SCaN) infrastructure, using state of the art simulation tools developed by the SCaN Center for Engineering, Networks, Integration, and Communications (SCENIC) laboratory at NASA Glenn Research Center to obtain the closest possible representation of the real operating environment. Reinforcement learning was then used to train an agent inside this environment to maximize data throughput. The simulation environment contained a single actor in low earth orbit capable of communicating with twenty-five ground stations that compose the Near-Earth Network (NEN). Initial experiments showed promising training results, so additional complexity was added by augmenting simulation data with link fading profiles obtained from real communication events with the International Space Station. A grid search was performed to find the optimal hyperparameters and model architecture for the agent. Using the results of the grid search, an agent was trained on the augmented training data. Testing shows that the agent performs well inside the training environment and can be used as a foundation for future studies with added complexity and eventually tested in the real space environment.

Schubert, Carson D.↗