SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This viewgraph presentation contains information about (1) Framework for multi-body dynamics - Geometry Manipulation Protocol (GMP), (2) Simulation procedure using Chimera Grid Tools (CGT) and OVERFLOW-2 (3) Further recent developments in Chimera Grid Tools OVERGRID, Grid modules, Script library and (4) Future work.
When electrified transit systems make grid aware choices, improved social welfare is achieved by reducing grid stress, reducing system loss, and minimizing power quality issues. Electrifying transit fleet has numerous challenges like non availability of buses during charging, varying charging costs and so on, that are related the electric grid behavior. However, transit systems do not have access to the information about the co-evolution of the grid's power flow and therefore cannot account for the power grid's needs in its day-to-day operation. In this paper we propose a framework of transportation-grid co-simulation, analyzing the spatio-temporal interaction between the transit operations with electric buses and the power distribution grid. Real-world data for a day's traffic from Chattanooga city's transit system is simulated in SUMO and integrated with a realistic distribution grid simulation (using GridLAB-D) to understand the grid impact due to transit electrification. Charging information is obtained from the transportation simulation to feed into grid simulation to assess the impact of charging. We also discuss the impact to the grid with higher degree of transit electrification that further necessitates such an integrated transportation-grid co-simulation to operate the integrated system optimally. Our future work includes extending the platform for optimizing the charging and trip assignment operations.
Microgrids help facilitate the integration of renewable energy in distribution-level grids and increase the resilience of the electric grid to extreme weather, especially in rural areas. Compared to traditional microgrids, a networked microgrid leverages multiple grid-forming sources to form a potential meshed grid and is more flexible in operation. This paper demonstrates the simulation modeling of an actual networked microgrid located in Adjuntas, Puerto Rico. The model contains representations of power inverters that connect the battery energy storage systems and photovoltaic generation systems to the networked microgrid and is capable of simulating fast grid transients as well as long-term operation of the networked microgrid. The modeling technique for power inverters allows the time-efficient simulation of the microgrid with a minimal penalty on model accuracy.
This paper presents a mesoscopic model to populate a set of sub-fleets of battery electric vehicles (BEVs) that simulates how grid power demand changes with varied mixes of different charging strategies. The proposed model considers real driving activities of vehicles based on the National Household Travel Survey (NHTS) and various charging strategies while parked (at home and in other locations). To be easily scaled to any fleet size, a mesoscopic model is implemented, where BEVs with the same physical properties and driving schedules and clustered together into sub-fleets. The degree of granularity is explored for stable results. Monte Carlo simulations were run to demonstrate the potential of this model's use in grid loads under various charging strategies. Furthermore, an example cost function was optimized to demonstrate finding optimal allocations of charging strategies.
This project implements coupled simulation of electric power grid transmission and distribution (T&D) by instantiating the co-simulation platform HELICS using a commercial distribution planning software CYMDIST1 by Eaton Inc., a power management company, and by leveraging the open-source transmission grid simulator GridDyn and LLNL’s High Performance Computing (HPC). The project is funded by DOE’s Technology Commercialization Funding (TCF) and supported by DOE Office of Technology Transitions (OTT).
Sandia National Laboratories has an existing capability for hybrid control systems testing called SCEPTRE. This article proposes an architecture to add dynamic simulation capability for the underlying physical process (e.g. the power grid). Dynamic simulation for SCEPTRE will enable very accurate simulation, and allow the full integration of analog control systems hardware.
Implicit large-eddy simulations of the self-similar regions of two compressible mixing layers at high Reynolds number and convective Mach numbers of 0.381 and 0.690 were carried out. Experimental data was used as the input into a synthetic eddy method turbulent inflow boundary condition to initiate simulations. Computational grids with nearly isotropic cells were used and the grid counts ranged from 24 to 270 million points. Agreement with experimental velocity and Reynolds stress data was good given the simplifications of the simulation. The decrease in mixing layer growth rate with increasing compressibility compares well with the literature. Budgets for the Reynolds stress transport equation were extracted from the simulations. It was demonstrated that the grids used did not resolve the dissipation term at the experimental Reynolds number. However, low Reynolds number simulations indicate that the numerical dissipation replaces the physical dissipation and the dissipation term in the transport equation can be represented by the summation term that quantifies the imbalance in the budget. The budgets for the high Reynolds number mixing layers agree well with previous low Reynolds number studies; indicating that the primary mechanism for the reduced growth rate and increased anistropy of the Reynolds stress tensor due to compressibility is a reduction in the pressure-strain and production terms. A new scaling based on the magnitude of the Reynolds stress tensor is proposed. This scaling provides better relative comparisons between the data and shows that the decrease in pressure-strain in the shear component is the primary driver in reducing the transverse normal stress.
It is essential to design a propulsion powertrain real-time simulator using the hardware-in-the-loop (HIL) system that emulates an electrified aircraft propulsion (EAP) systems power grid. This simulator would enable us to facilitate in-depth understanding of the system principles, to validate system model analysis and performance prediction, and to demonstrate the proof-of-concept of the EAP electrical system. This paper describes how subscale electrical machines with their controllers can mimic the power components in an EAP powertrain. In particular, three powertrain emulations are presented to mimic 1) a gas turbo-=shaft engine driving a generator, consisting of two permanent magnet (PM) motors with brushless motor drives, coupled by a shaft, 2) a motor driving a propulsive fan, and 3) a turbo-shaft engine driven fan (turbofan engine) operation. As a first step towards the demonstration, experimental dynamic characterization of the two motor drive systems, coupled by a mechanical shaft, were performed. The previously developed analytical motor models1 were then replaced with the experimental motor models to perform the real-time demonstration in the predefined flight path profiles. This technique can convert the plain motor system into a unique EAP power grid emulator that enables rapid analysis and real-time simulation performance using hardware-in-the-loop (HIL).
We study stability of gas accretion in active galactic nuclei (AGNs). Our grid-based simulations cover a radial range from 0.1 to 200 pc, which may enable linking the galactic/cosmological simulations with small-scale black hole (BH) accretion models within a few hundreds of Schwarzschild radii. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive BH of 10{sup 8} M{sub Sun }. Our theoretical estimates and simulations show that for the X-ray luminosity, L{sub X} {approx} 0.008 L{sub Edd}, the gas is thermally and convectively unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then nonlinearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For L{sub X} = 0.015 L{sub Edd} or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For L{sub X} < 0.015 L{sub Edd}, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. who used smoothed particle hydrodynamic (SPH) simulations to tackle the same problem. Here, however, because we use a grid-based code to solve equations in one dimension and two dimensions, we are able to follow the gas dynamics at much higher spacial resolution and for longer time compared with the three-dimensional SPH simulations. One of the new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outward within the hot outflow. Therefore, these simulations may serve as an attractive model for the so-called narrow-line region in AGNs.
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.
Snowmass is a US long-term planning study for the high-energy community by the American Physical Society's Division of Particles and Fields. For its simulation studies, opportunistic resources are harnessed using the Open Science Grid infrastructure. Late binding grid technology, GlideinWMS, was used for distributed scheduling of the simulation jobs across many sites mainly in the US. The pilot infrastructure also uses the Parrot mechanism to dynamically access CvmFS in order to ascertain a homogeneous environment across the nodes. This report presents the resource usage and the storage model used for simulating large statistics Standard Model backgrounds needed for Snowmass Energy Frontier studies.
Power system dynamic simulation computes the system response to a sequence of large disturbance, such as sudden changes in generation or load, or a network short circuit followed by protective branch switching operation. It consists of a large set of differential and algebraic equations, which is computational intensive and challenging to solve using single-processor based dynamic simulation solution. High-performance computing (HPC) based parallel computing is a very promising technology to speed up the computation and facilitate the simulation process. This paper presents two different parallel implementations of power grid dynamic simulation using Open Multi-processing (OpenMP) on shared-memory platform, and Message Passing Interface (MPI) on distributed-memory clusters, respectively. The difference of the parallel simulation algorithms and architectures of the two HPC technologies are illustrated, and their performances for running parallel dynamic simulation are compared and demonstrated.
The scalar-filtered mass-density-function methodology is employed for large-eddy simulation of two swirl-stabilized nonpremixed flames. These are low-swirl (SM1) and high-swirl (SMA2) Sydney methane flames, both of which have been the subject of detailed laboratory measurements. Combustion chemistry is modeled via a flamelet model for the low-swirl flame, and a detailed finite-rate kinetics model for the high-swirl flame. The scalar-filtered mass-density-function is simulated by a Lagrangian Monte Carlo method on a domain discretized by unstructured grids. The simulated results are assessed via comparison with experimental data and show very good agreement. This demonstrates the capability of scalar-filtered mass-density-function for large-eddy simulation of complex flows and warrants future applications of the methodology for large-eddy simulation of practical combustor configurations.
Abstract In mesoscale and global atmospheric simulations with large horizontal domains, strong horizontal flow across the grid is often unavoidable, but its effects on cloud‐topped boundary layers have received comparatively little study. Here the effects of cross‐grid flow on large‐eddy simulations of stratocumulus and trade‐cumulus marine boundary layers are studied across a range of grid resolutions (horizontal × vertical) between 500 m × 20 m and 35 m × 5 m. Three cases are simulated: DYCOMS nocturnal stratocumulus, BOMEX trade cumulus, and a GCSS stratocumulus‐to‐trade cumulus case. Simulations are performed with a stationary grid (with 4–8 m s −1 horizontal winds blowing through the cyclic domain) and a moving grid (equivalent to subtracting off a fixed vertically uniform horizontal wind) approximately matching the mean boundary‐layer wind speed. For stratocumulus clouds, cross‐grid flow produces two primary effects on stratocumulus clouds: a filtering of fine‐scale resolved turbulent eddies, which reduces stratocumulus cloud‐top entrainment, and a vertical broadening of the stratocumulus‐top inversion which enhances cloud‐top entrainment. With a coarse (20 m) vertical grid, the former effect dominates and leads to strong increases in cloud cover and LWP, especially as horizontal resolution is coarsened. With a finer (5 m) vertical grid, the latter effect is stronger and leads to small reductions in cloud cover and LWP. For the BOMEX trade cumulus case, cross‐grid flow tends to produce fewer and larger clouds with higher LWP, especially for coarser vertical grid spacing. The results presented are robust to choice of scalar advection scheme and Courant number.
A state-of-the-art, grid-convergent simulation methodology was applied to three-dimensional calculations of a single-cylinder optical engine. A mesh resolution study on a sector-based version of the engine geometry further verified the RANS-based cell size recommendations previously presented by Senecal et al. (“Grid Convergent Spray Models for Internal Combustion Engine CFD Simulations,” ASME Paper No. ICEF2012-92043). Convergence of cylinder pressure, flame lift-off length, and emissions was achieved for an adaptive mesh refinement cell size of 0.35 mm. Furthermore, full geometry simulations, using mesh settings derived from the grid convergence study, resulted in excellent agreement with measurements of cylinder pressure, heat release rate, and NOx emissions. On the other hand, the full geometry simulations indicated that the flame lift-off length is not converged at 0.35 mm for jets not aligned with the computational mesh. Further simulations suggested that the flame lift-off lengths for both the nonaligned and aligned jets appear to be converged at 0.175 mm. With this increased mesh resolution, both the trends and magnitudes in flame lift-off length were well predicted with the current simulation methodology. Good agreement between the overall predicted flame behavior and the available chemiluminescence measurements was also achieved. Our present study indicates that cell size requirements for accurate prediction of full geometry flame lift-off lengths may be stricter than those for global combustion behavior. This may be important when accurate soot predictions are required.
Issues associated with analysis of 'icing effects' on airfoil and wing performances are discussed, along with accomplishments and efforts to overcome difficulties with ice. Because of infinite variations of ice shapes and their high degree of complexity, computational 'icing effects' studies using available software tools must address many difficulties in geometry acquisition and modeling, grid generation, and flow simulation. The value of each technology component needs to be weighed from the perspective of the entire analysis process, from geometry to flow simulation. Even though CFD codes are yet to be validated for flows over iced airfoils and wings, numerical simulation, when considered together with wind tunnel tests, can provide valuable insights into 'icing effects' and advance our understanding of the relationship between ice characteristics and their effects on performance degradation.
The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.
Climate change and increased electrification of space and water heating in buildings can significantly affect future electricity demand and hourly demand profiles, which has implications for electric grid greenhouse gas emissions and capacity requirements. We use EnergyPlus to quantify building energy demand under historical and under several climate change projections of 32 kinds of building prototypes in 16 different climate zones of California and imposed these impacts on a year 2050 electric grid configuration by simulation in the Holistic Grid Resource Integration and Deployment (HIGRID) model. We find that climate change only prompted modest increases in grid resource capacity and negligible difference in greenhouse gas emissions since the additional electric load generally occurred during times with available renewable generation. Heating electrification, however, prompted a 30-40% reduction in greenhouse gas emissions but required significant grid resource capacity increases, due to the higher magnitude of load increases and lack of readily available renewable generation during the times when electrified heating loads occurred. Overall, this study translates climate change and electrification impacts to system-wide endpoint impacts on future electric grid configurations and highlights the complexities associated with translating building-level impacts to electric system-wide impacts.