SEARCH · Engineering Papers

Results for “large grids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Analytical Models of Frequency and Voltage in Large-Scale All-Inverter Power Systems

Low-order frequency response models for power systems have a decades-long history in optimization and control problems such as unit commitment, economic dispatch, and wide-area control. With a few exceptions, these models are built upon the Newtonian mechanics of synchronous generators, assuming that the frequency dynamics across a system are approximately homogeneous, and assume the dynamics of nodal voltages for most operating conditions are negligible, and thus are not directly computed at all buses. As a result, the use of system frequency models results in the systematic underestimation of frequency minimum nadir and maximum RoCoF, and provides no insight into the reactive power-voltage dynamics. This paper proposes a low-order model of both frequency and voltage response in grid-forming inverter-dominated power systems. The proposed model accounts for spatial-temporal variations in frequency and voltage behavior across a system and as a result, demonstrates the heterogeneity of frequency response in future renewable power systems. Electromagnetic transient (EMT) simulations are used to validate the utility, accuracy, and computational efficiency of these models, setting the basis for them to serve as fast, scalable alternatives to EMT simulation, especially when dealing with very large-scale systems, for both planning and operational studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Evaluation of Adaptive Volt-VAR to Mitigate PV Impacts [Slides]

Distributed generation (DG) sources like photovoltaic (PV) systems with advanced inverters are able to perform grid-support functions, like autonomous Volt-VAR that attempts to mitigate voltage issues by injecting or consuming reactive power. However, the Volt-VAR function operates with VAR priority, meaning real power may be curtailed to provide additional reactive power support. Since some locations on the grid may be more prone to higher voltages than others, PV systems installed at those locations may be forced to curtail more power, adversely impacting the value of that PV system. Adaptive Volt-VAR (AVV) could be implemented as an alternative, whereby the Volt-VAR reference voltage changes over time, but this functionality has not been well-explored in the literature. In this work, the potential benefits and grid impacts of AVV were investigated using yearlong quasi-static time-series (QSTS) simulations. After testing a variety of allowable AVV settings, we found that even with aggressive settings AVV resulted in <0.01% real power curtailment and significantly reduced the reactive power support required from the PV inverter compared to conventional Volt-VAR but did not provide much mitigation for extreme voltage conditions. The reactive power support provided by AVV was injected to oppose large deviations in voltage (in either direction), indicating that it could be useful for other applications like reducing voltage flicker or minimizing interactions with other voltage regulating devices.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

The Implementation of Framework for Improvement by Vertical Enhancement Into Energy Exascale Earth System Model

Abstract The low cloud bias in global climate models (GCMs) remains an unsolved problem. Coarse vertical resolution in GCMs has been suggested to be a significant cause of low cloud bias because planetary boundary layer parameterizations cannot resolve sharp temperature and moisture gradients often found at the top of subtropical stratocumulus layers. This work aims to ameliorate the low cloud problem by implementing a new computational method, the Framework for Improvement by Vertical Enhancement (FIVE), into the Energy Exascale Earth System Model (E3SM). Three physics schemes representing microphysics, radiation, and turbulence as well as vertical advection are interfaced to vertically enhanced physics (VEP), which allows for these processes to be computed on a higher vertical resolution grid compared to the rest of the E3SM model. We demonstrate the better representation of subtropical boundary layer clouds with FIVE while limiting additional computational cost from the increased number of levels. When the vertical resolution approaches the large eddy simulation‐like vertical resolution in VEP, the climatological low cloud amount shows a significant increase of more than 30% in the southeastern Pacific Ocean. Using FIVE to improve the representation of low‐level clouds does not come with any negative side effects associated with the simulation of mid‐ and high‐level cloud and precipitation, that can occur when running the full model at higher vertical resolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling and Simulation of Power-Electronic Inverters in Analog Electronic Circuit Simulators

This paper demonstrates how equivalent-circuit representations of grid-following power-electronic inverters can be realized within a SPICE-based development environment using common circuit components and VerilogA code. This facilitates computationally lean simulations of inverter networks leveraging the strengths of SPICE in large-scale simulations. We validate the approach with time-domain simulations for a modified version of the IEEE standard 118-bus system modeled in Virtuoso (a SPICE-based solver). Simulation results are compared-focusing on accuracy and computation speed-with results from a commonly used power-electronics simulation package. We note a significant decrease in simulation time with comparable signal resolution when simulating the network in SPICE.

Billmeyer, Ryan↗

DOE OSTI2022

Economic Impact of Permitting Timelines on Produced Geothermal Power

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Transmission Innovation Symposium: Modernizing the U.S. Electrical Grid

The foundation of the United States Department of Energy (DOE) Transmission Reliability research program was established 20 years ago by a series of commissioned white papers. Those white papers described the dramatic institutional and regulatory changes that the U.S. electricity transmission grid was undergoing at the time and articulated the technical challenges that these changes created. The challenges outlined in the white papers were the basis for the initial research goals of the DOE Transmission Reliability program. To a large extent, the reliability research needs outlined in the original white papers have now been met. As a result, now is an appropriate time to step back and review the technical challenges that the industry currently faces and to use those challenges as the basis for identifying the next set of targets for DOE’s transmission-related research and development (R&D) programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Assessment of Economic Impact of Permitting Timelines on Produced Geothermal Power in Imperial County, California

Despite having a large geothermal power potential in the United States, only a small fraction has been developed for power generation. Various barriers, including technical, financial, and regulatory permit delays, are attributed to lower contribution of geothermal energy in the national grid. Unpredictable environmental reviews and permitting timelines are some of the non-technical barriers that can cause delays in geothermal exploration and utilization plans. Our study shows that the geothermal permitting timelines can vary from six months to several years, depending on the presence or absence of biological resources, cultural resources, and sensitive environmental issues at the project site. The potential impacts of these permit barriers can range from investors abandoning geothermal development to making the product (i.e., electricity) more expensive and uncompetitive. In this study, we conducted economic analysis to assess the impact of permitting timelines on cost of produced electricity from geothermal resources using data from existing geothermal plants as well as prospective sites. In this paper, we present collected timelines data, approach, and results of economic impact of permitting timelines on geothermal power. We evaluated the various environmental management and permit review processes by considering a hypothetical geothermal project in the Salton Sea Known Geothermal Resource Area. Because of the variety of the biological and environmental issues and the involvement of local, state, and federal agencies with overlapping jurisdictions, this project could go through one of the many California Environmental Quality Act (CEQA) and National Environmental Policy Act (NEPA) review scenarios that range from the least to the most complex in its circumstances. The fastest CEQA/NEPA review timelines would have the project completed in six years. In contrast, the project would substantially need longer time to complete if it were located in an area with significant environmental resources or cultural issues that required permitting from various agencies. With increasing project completion timelines, the simplified levelized cost of electricity (sLCOE) can be 4 to 11% higher with longer CEQA/NEPA review timelines than the sLCOE value with the fastest CEQA/NEPA review timeline. Lengthier CEQA/NEPA review timelines could also result in loss of $64 million to $227 million in potential revenue. Such significant economic impacts could determine the success of a geothermal project.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Why is My Zero Energy Home Not a Zero Carbon Home?

For years, carbon calculations were done very simply. The method of calculation was to take annual totals of energy consumption and multiply by an average emission factor, either for the grid serving a project or for a larger region (e.g. an EPA eGRID sub region). The level of accuracy of this approximation was reasonably good, although the issue of accuracy was not, to our knowledge, tested. And the data required were minimal – just a year’s worth of bills for each fuel and one lookup factor. But this method assures that a net zero energy home is automatically a net zero carbon home because zero times any possible emission factor is still zero. Starting in the early 2010s, things changed – grids were starting to rely more and more heavily on renewables, and the difference was showing up on aggregate load curves. This was perhaps noticed first in California, where aggressive renewable policies led to significant renewable power generation large enough to affect the overall shape of the diurnal load curve for the Independent Systems Operator.

14 SOLAR ENERGY↗

DOE OSTI2022

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗

DOE OSTI2024

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗

DOE OSTI2023

Investigation on Tuning Power–Frequency Droop for Improved Grid–Forming Inverter and Synchronous Generator Transient Load Sharing

Poor transient load sharing has been observed during the parallel operation of synchronous generators (SGs) and droop-controlled grid-forming (GFM) inverter-based resources (IBRs) in islanded microgrids. This can result in overcurrent protections, causing the inverters to trip, which leads to the collapse of voltage and frequency. In this article, small-signal analysis supported by electromagnetic transients (EMTs) simulation is performed on a detailed model of a microgrid containing IBRs and an SG to identify underdamped modes and their sensitivities to the power–frequency (P-F) droop parameter. Further EMT simulation of a more complex microgrid model with multiple GFM IBRs and SGs is performed to determine whether trends identified from the modal analysis extend to cases where the number of SGs in operation may vary. The objective is to investigate the effectiveness of tuning this parameter to reduce the likelihood of overcurrents in response to load and generation disturbances. The work uncovers both fast and slow modes of concern. Primarily, the findings show that increasing the P-F droop value improves the damping of the slower mode while degrading the damping of the faster mode, leading to an upper limit on the damping improvement of the GFM IBR output current. In EMT simulations of the more complex microgrid, the greatest damping improvement still exhibits significant maximum overshoot. Furthermore, the optimal value of the P-F droop parameter is sensitive to the number of SGs in operation. The conclusion is that tuning the P-F droop parameter is not an effective means to avoid the overcurrent tripping of GFM IBRs during response to large load or generation disturbances when operating in parallel to SGs in islanded microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Application of Numerical Atmospheric Models

In this study, a numerical modeling framework for simulating extreme storm events was established using the Weather Research and Forecasting (WRF) model. Such a framework is necessary for the derivation of engineering parameters such as probable maximum precipitation that are the cornerstone of large water-management infrastructure design. Here, this framework was built based on a heavy storm that occurred in Nashville, Tennessee (U.S.), in 2010, and verified using two other extreme storms. To achieve the optimal setup, several combinations of model resolutions, initial/boundary conditions (IC/BC), cloud microphysics, and cumulus parameterization schemes were evaluated using multiple metrics of precipitation characteristics. The evaluation suggests that WRF is most sensitive to the IC/BC option. Simulation generally benefits from finer resolutions up to 5 km. At the 15 km level, NCEP2 IC/BC produces better results, whereas NAM IC/BC performs best at the 5 km level. The recommended model configuration from this study is: NAM or NCEP2 IC/BC (depending on data availability), 15 km or 15–5 km nested grids, Morrison microphysics, and Kain-Fritsch cumulus schemes. Validation of the optimal framework suggests that these options are good starting choices for modeling extreme events similar to the test cases. This optimal framework is proposed in response to emerging engineering demands of extreme storm event forecasting and analyses for design, operations, and risk assessment of large water infrastructures.

Chen, Xiaodong↗

DOE OSTI2024

A collision-based hybrid method for the BGK equation

In this article, we apply the collision-based hybrid method introduced by Hauck and McClarren to the Boltzmann equation with the BGK operator and a hyperbolic scaling. An implicit treatment of the source term is used to handle stiffness associated with the BGK operator. Although it helps the numerical scheme become stable with a large time step size, it is still not obvious to achieve the desired order of accuracy due to the relationship between the size of the spatial cell and the mean free path. Without asymptotic preserving property, a very restricted grid size is required to resolve the mean free path, which is not practical. Our approaches are based on the noncollision-collision decomposition of the BGK equation. We introduce the arbitrary order of nodal discontinuous Galerkin (DG) discretization in space with a semi-implicit time-stepping method; we employ the backward Euler time integration for the uncollided equation and the 2nd order predictor-corrector scheme for the collided equation, i.e., both source terms in uncollided and collided equations are treated implicitly and only streaming term in the collided equation is solved explicitly. This improves the computational efficiency without the complexity of the numerical implementation. Numerical results are presented for various Knudsen numbers to present the effectiveness and accuracy of our hybrid method. Also, we compare the solutions of the hybrid and non-hybrid schemes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Enabling kilometer-scale E3SM land model simulation over North America: A new integrated framework solution

This study introduces a novel framework designed to enhance the performance, scalability, and portability of the kilometer-scale E3SM Land Model (km-ELM) within the E3SM modeling infrastructure. By seamlessly integrating cutting-edge data tools, we address existing challenges such as slow performance, limited scalability, and difficulties in software integration in current data-driven ELM simulation over large geographic areas. Our innovative approach leverages the KiloCraft data toolkit to generate unified inputs for simulations ranging from a single-cite case, to a 72,083-cell regional case to a continental configuration encompassing 21.6 million land grid cells at a 1 km × 1 km resolution. We conduct extensive strong- and weak-scaling experiments on three state-of-the-art supercomputers, utilizing up to 100,800 CPU cores across 2400 compute nodes to evaluate end-to-end metrics including wall-clock time, simulation-years-per-day (SYPD), initialization costs, and I/O throughput. Our results reveal the land (LND) component’s efficient scaling, demonstrating near-ideal weak scaling and strong-scaling parallel efficiencies reaching up to 87% at 50,400 cores. We confirm portability and reproducibility through bitwise-equivalent outputs across different machines using identical inputs over supported machines. Notably, at extreme scales, we identify I/O as a critical bottleneck and that leads to effective solution with the SCORPIO/ADIOS stack. Collectively, these findings validate the deployment of km-ELM at a continental scale with high parallel efficiency and provide essential guidance on configuration, decomposition, and I/O settings for optimized kilometer-scale land simulations in E3SM. This work emphasizes the innovative design and practical solutions that enhance the operational capabilities of km-ELM, focusing on software performance and scalability while leaving detailed scientific evaluations of simulated land processes for future investigations.

E3SM land model (ELM), km-ELM, scalability, perfor↗

DOE OSTI2021

Storage Futures Study: The Challenge of Defining Long-Duration Energy Storage

This document explores the definition of "long duration" as applied to energy storage. Given the growing use of this term, a uniform definition could aid in communication and consistency among various stakeholders. There is large and growing use of the ARPA-E definition of greater than 10 hours. However, the term "long-duration energy storage" is often used as shorthand for storage with sufficient duration to provide firm capacity and support grid resource adequacy. The actual duration needed for this application varies significantly from as little as a few hours to potentially multiple days. This leads to the unsatisfying conclusion that there cannot be a simple, uniform, and static definition of long duration storage that captures its ability to provide firm capacity and also aids consistent communication.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗