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 1,279 records · Page 71

Modular Isolated Soft-Switching Medium Voltage String Inverter for Large-Scale PV Farm

This paper presents a single-stage modular isolated soft-switching medium voltage string inverter (MVSI) with tri-port configuration to interconnect PV and storage to medium voltage (MV) grid. The modular MVSI provides advanced functionalities such as energy storage elements integration without additional converters, connection to MV grid resulting in lower distribution losses in the PV farm, high-frequency galvanic isolation, independent control of PV, storage and grid power, high modularity and scalability. The topology of the MVSI, operation principle and model-predictive priority-switching (MPPS) control to achieve fast dynamic balancing across stacked modules are discussed. The independent power control is validated in PLECS simulation and the operation of stacked MVSI with MPPS control is verified through Hardware-in-the-loop (HIL) simulation. Here, the hardware prototype of 25 kVA MVSI building block is developed. 18 kW experimental results of single building block verifying tri-ports operation and 20 kW experimental results of two modules stacked operation are presented. A clean sinusoidal waveform with negligibly low THD in both tests is obtained. A low dv/dt of <; 1kV/us and soft-switching are achieved.

30 DIRECT ENERGY CONVERSION↗

Analytic Considerations and Design Basis for the IEEE Distribution Test Feeders

One of the greatest challenges for the analysis of distribution systems is the high level of diversity across the world. For nearly 20 years the Test Feeder Working Group of the Distribution System Analysis Subcommittee has been developing distribution test feeders. The purpose of these test feeders is to provide openly-available models of distribution systems that reflect the diversity in systems and the analytic challenges they present. Because of their utility and accessibility, the test feeders have been used for a wide range of research, some of which they were not originally intended to cover. The purpose of this paper is to provide an overview of the variety of existing distribution systems, and to clarify the specific analytic challenges that each test feeder was originally designed to examine. Additionally, the paper will provide guidance on which feeders are best suited for various types of analysis. The purpose of this paper is to provide the original intent of the Working Group and to inform researchers so that they may make an informed decision on which of the test feeders are most appropriate for their work.

Schneider, Kevin P.↗

THOR: A NEW AND FLEXIBLE GLOBAL CIRCULATION MODEL TO EXPLORE PLANETARY ATMOSPHERES

We have designed and developed, from scratch, a global circulation model (GCM) named THOR that solves the three-dimensional nonhydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the “pole problem” (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing an icosahedral grid. Irregularities in the grid, which lead to grid imprinting, are smoothed using the “spring dynamics” technique. We validate our implementation of spring dynamics by examining calculations of the divergence and gradient of test functions. To prevent the computational time step from being bottlenecked by having to resolve sound waves, we implement a split-explicit method together with a horizontally explicit and vertically implicit integration. We validate our GCM by reproducing the Earth and hot-Jupiter-like benchmark tests. THOR was designed to run on graphics processing units (GPUs), which allows for physics modules (radiative transfer, clouds, chemistry) to be added in the future, and is part of the open-source Exoclimes Simulation Platform (www.exoclime.org).

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of Computed and Measured Vortex Evolution for a UH-60A Rotor in Forward Flight

A Computational Fluid Dynamics (CFD) simulation using the Navier-Stokes equations was performed to determine the evolutionary and dynamical characteristics of the vortex flowfield for a highly flexible aeroelastic UH-60A rotor in forward flight. The experimental wake data were acquired using Particle Image Velocimetry (PIV) during a test of the fullscale UH-60A rotor in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel. The PIV measurements were made in a stationary cross-flow plane at 90 deg rotor azimuth. The CFD simulation was performed using the OVERFLOW CFD solver loosely coupled with the rotorcraft comprehensive code CAMRAD II. Characteristics of vortices captured in the PIV plane from different blades are compared with CFD calculations. The blade airloads were calculated using two different turbulence models. A limited spatial, temporal, and CFD/comprehensive-code coupling sensitivity analysis was performed in order to verify the unsteady helicopter simulations with a moving rotor grid system.

Evolution↗

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster↗

Quantifying Power Distribution System Resiliency Using Code-Based Metric

It is essential to improve the resiliency of power distribution systems (PDS) with the increase in extreme weather events, number of malicious threats and consumers’ need for higher reliability. Here, this paper provides a formal approach to evaluate the operational resiliency of PDS, and quantify the resiliency of a system using a code-based metric. A combination of steady state and dynamic simulation tools is used to determine the resiliency metric. Dynamic simulation tools help with analyzing impact of short-term events, which might affect operational resiliency in long term. A dynamic optimization algorithm for changing operating criteria to increase the sustainability of the most critical loads has been proposed. The proposed theoretical approach is validated using a simple power distribution system model and simulation results demonstrate the ability to quantify the resiliency using the proposed code-based metric. The time-dependent quantification of resiliency has been demonstrated on a test system of two connected CERTS microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗

Particle Morphology and Size Results from the Smoke Aerosol Measurement Experiment-2

Results are presented from the Reflight of the Smoke Aerosol Measurement Experiment (SAME-2) which was conducted during Expedition 24 (July-September 2010). The reflight experiment built upon the results of the original flight during Expedition 15 by adding diagnostic measurements and expanding the test matrix. Five different materials representative of those found in spacecraft (Teflon, Kapton, cotton, silicone rubber and Pyrell) were heated to temperatures below the ignition point with conditions controlled to provide repeatable sample surface temperatures and air flow. The air flow past the sample during the heating period ranged from quiescent to 8 cm/s. The smoke was initially collected in an aging chamber to simulate the transport time from the smoke source to the detector. This effective transport time was varied by holding the smoke in the aging chamber for times ranging from 11 to 1800 s. Smoke particle samples were collected on Transmission Electron Microscope (TEM) grids for post-flight analysis. The TEM grids were analyzed to observe the particle morphology and size parameters. The diagnostics included a prototype two-moment smoke detector and three different measures of moments of the particle size distribution. These moment diagnostics were used to determine the particle number concentration (zeroth moment), the diameter concentration (first moment), and the mass concentration (third moment). These statistics were combined to determine the diameter of average mass and the count mean diameter and, by assuming a log-normal distribution, the geometric mean diameter and the geometric standard deviations can also be calculated. Overall the majority of the average smoke particle sizes were found to be in the 200 nm to 400 nm range with the quiescent cases producing some cases with substantially larger particles.

Urban, David L.↗

Computational Assessment of a 3-Stage Axial Compressor Which Provides Airflow to the NASA 11- by 11-Foot Transonic Wind Tunnel, Including Design Changes for Increased Performance

A 24 foot diameter 3-stage axial compressor powered by variable-speed induction motors provides the airflow in the closed-return 11- by 11-Foot Transonic Wind Tunnel (11-Foot TWT) Facility at NASA Ames Research Center at Moffett Field, California. The facility is part of the Unitary Plan Wind Tunnel, which was completed in 1955. Since then, upgrades made to the 11-Foot TWT such as flow conditioning devices and instrumentation have increased blockage and pressure loss in the tunnel, somewhat reducing the peak Mach number capability of the test section. Due to erosion effects on the existing aluminum alloy rotor blades, fabrication of new steel rotor blades is planned. This presents an opportunity to increase the Mach number capability of the tunnel by redesigning the compressor for increased pressure ratio. Challenging design constraints exist for any proposed design, demanding the use of the existing driveline, rotor disks, stator vanes, and hub and casing flow paths, so as to minimize cost and installation time. The current effort was undertaken to characterize the performance of the existing compressor design using available design tools and computational fluid dynamics (CFD) codes and subsequently recommend a new compressor design to achieve higher pressure ratio, which directly correlates with increased test section Mach number. The constant cross-sectional area of the compressor leads to highly diffusion factors, which presents a challenge in simulating the existing design. The CFD code APNASA was used to simulate the aerodynamic performance of the existing compressor. The simulations were compared to performance predictions from the HT0300 turbomachinery design and analysis code, and to compressor performance data taken during a 1997 facility test. It was found that the CFD simulations were sensitive to endwall leakages associated with stator buttons, and to a lesser degree, under-stator-platform flow recirculation at the hub. When stator button leakages were modeled, pumping capability increased by over 20 of pressure rise at design point due to a large reduction in aerodynamic blockage at the hub. Incorporating the stator button leakages was crucial to matching test data. Under-stator-platform flow recirculation was thought to be large due to a lack of seals. The effect of this recirculation was assessed with APNASA simulations recirculating 0.5, 1, and 2 of inlet flow about stators 1 and 2, modeled as axisymmetric mass flux boundary conditions on the hub before and after the vanes. The injection of flow ahead of the stators tended to re-energize the boundary layer and reduce hub separations, resulting in about 3 increased stall margin per 1 of inlet flow recirculated. In order to assess the value of the flow recirculation, a mixing plane simulation of the compressor which gridded the under-stator cavities was generated using the ADPAC CFD code. This simulation indicated that about 0.65 of the inlet flow is recirculated around each shrouded stator. This collective information was applied during the redesign of the compressor. A potential design was identified using HT0300 which improved overall pressure ratio by removing pre-swirl into rotor 1, replacing existing NASA 65 series rotors with double circular arc sections, and re-staggering rotors and the existing stators. The performance of the new design predicted by APNASA and HT0300 is compared to the existing design.

Turbomachinery↗

A Method for Flow Simulation About Complex Geometries Using Both Structured and Unstructured Grids

A computational fluid dynamics code which utilizes both structured and unstructured grids was developed. The objective of this study was to develop and demonstrate the ability of such a code to achieve solutions about complex geometries in two dimensions. An unstructured grid generator and flow solver were incorporated into the PARC2D structured flow solver. This new unstructured grid generator capability allows for easier generation and manipulation of complex grids. Several examples of the grid generation capabilities are provided. The coupling of different grid topologies and the manipulation of individual grids is shown. Also, grids for realistic geometries, a NACA 0012 airfoil and a wing/nacelle installation, were created. The flow over a NACA 0012 airfoil was used as a test case for the flow solver. Eight separate cases were run. They were both the inviscid and viscous solutions for two freestream Mach numbers and airfoil angle of attacks of 0 to 3.86 degrees. The Mach numbers chosen were for a subsonic case, Mach 0.6, and a case where supersonic regions and a shock wave exists, Mach 0.8. These test case conditions were selected to match experimentally obtained data for code comparison. The results show that the code accurately predicts the flow field for all cases.

Debonis, James R.↗

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗

Formulation of Three-Phase State Estimation Problem Using a Virtual Reference

State Estimation (SE) is the backbone of modern Energy Management System (EMS) due its capability of processing real time measurements and provide reliable information to system operators. Since its introduction to power systems in the 70’s, SE has been widely used in transmission networks. Distribution grids on the other hand lack sufficient number of real time measurements, and for that reason SE has not been widely implemented on these systems. The recent increase in the number of renewable energy sources connected to the grid at lower voltage levels, the advent of Distribution Automation (DA) and Smart Grids necessitate closer monitoring of distribution networks and thus forcing utilities to upgrade their operations and deploy Advanced Distribution Management Sytems (ADMS). Therefore, Distribution System State Estimation (DSE) is paramount to provide real time monitoring of active distribution grids. This papers investigates the formulation of the three-phase distribution state estimation where in the most general case, none of the system buses may have balanced three phase voltages. In such a case, the choice and definition of the slack bus need to be revisited. Most of the existing work arbitrarily assign the root bus of the distribution feeder as the balanced three-phase reference which may not be realistic for today’s active distribution systems. This paper addresses this shortcoming of the existing SE formulations by introducing a virtual reference bus. Here, the proposed approach is effectively validated by simulations carried out on distribution test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impact of Cyber Attacks on High Voltage DC Transmission Damping Control

Electric grids have evolved to being a huge cyberphysical system (CPS) that combines intelligent electronic devices (IEDs), computing resources, communication networks and electricity infrastructures. The wide attack surface and increasing dependence on communication make power systems vulnerable to cyber attacks. This paper studies the impact of cyber attacks on high voltage DC (HVDC) transmission oscillation damping control. Three kinds of cyber attack models are considered: timing attack , replay attack, and false data injection attack. The attacks are applied to the HVDC damping controller, which requires widearea measurements and long-distance communication. A modified IEEE New England 39-Bus AC-HVDC system is used as the test system. Simulation results have shown that all three kinds of attacks are capable of driving the AC-HVDC system into swing or unstable conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mapping Surface Cover Parameters Using Aggregation Rules and Remotely Sensed Cover Classes

A coupled model, which combines the Biosphere-Atmosphere Transfer Scheme (BATS) with an advanced atmospheric boundary-layer model, was used to validate hypothetical aggregation rules for BATS-specific surface cover parameters. The model was initialized and tested with observations from the Anglo-Brazilian Amazonian Climate Observational Study and used to simulate surface fluxes for rain forest and pasture mixes at a site near Manaus in Brazil. The aggregation rules are shown to estimate parameters which give area-average surface fluxes similar to those calculated with explicit representation of forest and pasture patches for a range of meteorological and surface conditions relevant to this site, but the agreement deteriorates somewhat when there are large patch-to-patch differences in soil moisture. The aggregation rules, validated as above, were then applied to remotely sensed 1 km land cover data set to obtain grid-average values of BATS vegetation parameters for 2.8 deg x 2.8 deg and 1 deg x 1 deg grids within the conterminous United States. There are significant differences in key vegetation parameters (aerodynamic roughness length, albedo, leaf area index, and stomatal resistance) when aggregate parameters are compared to parameters for the single, dominant cover within the grid. However, the surface energy fluxes calculated by stand-alone BATS with the 2-year forcing, data from the International Satellite Land Surface Climatology Project (ISLSCP) CDROM were reasonably similar using aggregate-vegetation parameters and dominant-cover parameters, but there were some significant differences, particularly in the western USA.

Arain, Altaf M.↗

DSMC Simulation of Separated Flows About Flared Bodies at Hypersonic Conditions

This paper describes the results of a numerical study of interacting hypersonic flows at conditions that can be produced in ground-based test facilities. The computations are made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is on Mach 10 flows about flared axisymmetric configurations, both hollow cylinder flares and double cones. The flow conditions are those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel. The range of flow conditions, model configurations, and model sizes provides a significant range of shock/shock and shock/boundary layer interactions at low Reynolds number conditions. Results presented will highlight the sensitivity of the calculations to grid resolution, contrast the differences in flow structure for hypersonic cold flows and those of more energetic but still low enthalpy flows, and compare the present results with experimental measurements for surface heating, pressure, and extent of separation.

Moss, James N.↗

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗

Efficient Simulation of Cascading Outages Using an Energy Function-Embedded Quasi-Steady-State Model

Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.

Guo, Zhenping [Univ. of Tennessee, Knoxville, TN (↗

Verification assessment of piston boundary conditions for Lagrangian simulation of compressible flow similarity solutions

This work is concerned with the use of similarity solutions of the compressible flow equations as benchmarks or verification test problems for finite-volume compressible flow simulation software. In practice, this effort can be complicated by the infinite spatial/temporal extent of many candidate solutions or “test problems.” Methods can be devised with the intention of ameliorating this inconsistency with the finite nature of computational simulation; the exact strategy will depend on the code and problem archetypes under investigation. For example, self-similar shock wave propagation can be represented in Lagrangian compressible flow simulations as rigid boundary-driven flow, even if no such “piston” is present in the counterpart mathematical similarity solution. The purpose of this work is to investigate in detail the methodology of representing self-similar shock wave propagation as a piston-driven flow in the context of various test problems featuring simple closed-form solutions of infinite spatial/temporal extent. The closed-form solutions allow for the derivation of similarly closed-form piston boundary conditions (BCs) for use in Lagrangian compressible flow solvers. Finally, the consequences of utilizing these BCs (as opposed to directly initializing the self-similar solution in a computational spatial grid) are investigated in terms of common code verification analysis metrics (e.g., shock strength/position errors and global convergence rates).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗