Simulating buoyancy-driven airflow in buildings by coarse-grid fast fluid dynamics
Not Available
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.
Not Available
This article proposes an approach to address the current and the aggregated active power control challenge for large-scale inverter-based resources subjected to partially or completely loss of inverters and grid voltage variations. To address this problem, a distributed active power mechanism is proposed which generates desired inverters current. Then an adaptive mechanism distributes the voltage control input automatically between inverters in response to the partial or complete loss of inverters. An L2-gain-based controller is designed for each inverter to track the desired current and rejects the grid voltage disturbance. Here, simulation results show a significant robust tracking for collective active power and current.
Thermal security region (TSR) is a powerful tool for monitoring and controlling the thermal security of bulk power grid with high penetration of renewable energy. One of the major challenges of TSR in practical application is to obtain the exact analytical expression of the TSR boundary (TSRB) because TSRB is a nonlinear hypersurface described by an implicit function. For this reason, the TSRB is usually approximated by hyperplane instead of getting the exact analytical expression. Traditionally, the coefficients of the hyperplane are estimated with a set of original TSRB points using least square estimation (LSE). However, LSE is a point estimation method, which is incapable of evaluating the quality of the estimation. In addition, the accuracy and computational efficiency of the hyperplane approximation are influenced by the number of TSRB points significantly. Here, a bootstrap based confidence interval estimation is proposed to estimate not only the coefficients of TSRB approximation hyperplane, but also the standard deviations and confidence intervals of the coefficients for evaluating the quality and reliability of the approximation results. First, empirical distribution functions (EDFs) of the coefficients of TSRB approximation hyperplane are approximated from a set of original TSRB points by using residual resampling bootstrap method. Then, the EDFs of the coefficients are employed to estimate the coefficients of the approximation hyperplane. Meanwhile, the standard deviations and confidence intervals of the estimated coefficients of the hyperplanes are also calculated for evaluating the quality and reliability of the approximation. The proposed approach is tested on the China Southern Power Grid (CSG). Results of simulations validate the accuracy and efficiency of the proposed approach in approximating the TSRB.
We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.
Hailstones and extreme precipitation generate substantial economic losses across the United States (US) and the globe. Their strong association with short-lived, intense convective storms poses a great challenge in predicting their future changes. Here, we conducted model simulations at 1.2 km grid spacing for severe convective storms with large hail and heavy precipitation that occurred in two typical types of synoptic-scale environments in spring seasons over the central US under both current and future climate conditions. We find that the responses of large hail (diameters >2.5 cm) to anthropogenic climate change (ACC) are markedly different between the hailstorms developed in the two types of synoptic-scale environments, with over 110% increase in large hail occurrences for the frontal systems, whereas less than 30% increase for the Great Plains low-level jet (GPLLJ) systems. This is explained by the larger increase in convective intensity and updraft width and a smaller increase in warm cloud depth in the frontal storms compared with the GPLLJ storms. Interestingly, the occurrences and intensity of heavy precipitation (rain rate >20 mm hr-1) in both types of systems are similarly sensitive to ACC (e.g., 40% and 33% increases in the occurrences for the frontal and GPLLJ systems, respectively). These results advance our knowledge of hail projection and have important implications for managing risks for future hail.
Lattices remain an attractive class of structures due to their design versatility; however, rapidly designing lattice structures with tailored or optimal mechanical properties remains a significant challenge. With each added design variable, the design space quickly becomes intractable. To address this challenge, research efforts have sought to combine computational approaches with machine learning (ML)-based approaches to reduce the computational cost of the design process and accelerate mechanical design. While these efforts have made substantial progress, significant challenges remain in (1) building and interpreting the ML-based surrogate models and (2) iteratively and efficiently curating training datasets for optimization tasks. Here, we address the first challenge by combining ML-based surrogate modeling and Shapley additive explanation (SHAP) analysis to interpret the impact of each design variable. We find that our ML-based surrogate models achieve excellent prediction capabilities (R 2 > 0.95) and SHAP values aid in uncovering design variables influencing performance. We address the second challenge by utilizing active learning-based methods, such as Bayesian optimization, to explore the design space and report a 5 × reduction in simulations relative to grid-based search. Collectively, these results underscore the value of building intelligent design systems that leverage ML-based methods for uncovering key design variables and accelerating design.
Simulation codes for galaxy formation and evolution take on board as many physical processes as possible beyond the standard gravitational and hydrodynamical physics. Most of this extra physics takes place below the resolution level of the simulations and is added in a ''sub-grid'' fashion. However, these sub-grid processes affect the macroscopic hydrodynamical properties of the gas and thus couple to the ''on-grid'' physics that is explicitly integrated during the simulation. In this paper, we focus on the link between partial ionization and the hydrodynamical equations. We show that the energy stored in ions and free electrons constitutes a potential energy term which breaks the linear dependence of the internal energy on temperature. Correctly taking into account ionization hence requires modifying both the equation of state and the energy-temperature relation. We implemented these changes in the cosmological simulation code GADGET2. As an example of the effects of these changes, we study the propagation of Sedov-Taylor shock waves through an ionizing medium. This serves as a proxy for the absorption of supernova feedback energy by the interstellar medium. Depending on the density and temperature of the surrounding gas, we find that up to 50% of the feedback energy is spent ionizing the gas rather than heating it. Thus, it can be expected that properly taking into account ionization effects in galaxy evolution simulations will drastically reduce the effects of thermal feedback. To the best of our knowledge, this potential energy term is not used in current simulations of galaxy formation and evolution.
Turbulence in the interstellar medium is anisotropic due to the ubiquitous magnetic fields. This anisotropy depends on the strength of the magnetic field and leaves an imprint on observations of spectral line maps. We use a grid of ideal magnetohydrodynamic simulations of driven turbulence and produce synthetic position–position–velocity maps to study the turbulence anisotropy in velocity channels of various resolutions. We found that the average structure function of velocity channels is aligned with the projection of the magnetic field on the plane of the sky. We also found that the degree of such anisotropy increases with the magnitude of the magnetic field. For thick velocity channels (low velocity resolution), the anisotropy is dominated by density, and the degree of anisotropy in these maps allows one to distinguish sub-Alfvénic and super-Alfvénic turbulence regimes, but it also depends strongly on the sonic Mach number. For thin channels (high velocity resolution), we find that the anisotropy depends less on the sonic Mach number. An important limitation of this technique is that it only gives a lower limit on the magnetic field strength because the anisotropy is related only to the magnetic field component on the plane of the sky. It can, and should, be used in combination with other techniques to estimate the magnetic field, such as the Fermi-Chandrasekhar method, anisotropies in centroids, Faraday rotation measurements, or direct line-of-sight determinations of the field from Zeeman effect observations.
ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.
ABSTRACT The reverse shock in the ejecta of core-collapse supernovae is potentially able to destroy newly formed dust material. In order to determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code astrobear in order to model a shock wave interacting with clumpy supernova ejecta. Dust motions and destruction rates were computed using our newly developed external, post-processing code paperboats, which includes gas drag, grain charging, sputtering, and grain–grain collisions. We have determined dust destruction rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes and clump gas density. We found that up to $30\,\mathrm{{{\ \rm per\ cent}}}$ of the carbon dust mass is able to survive the passage of the reverse shock if the initial grain size distribution is narrow with radii around ∼10–50 nm for high gas densities, or with radii around $\sim 0.5\!-\!1.5\,\mathrm{\mu m}$ for low and medium gas densities. Silicate grains with initial radii around 10–30 nm show survival rates of up to $40\,\mathrm{{{\ \rm per\ cent}}}$ for medium- and high-density contrasts, while silicate material with micron-sized distributions is mostly destroyed. For both materials, the surviving dust mass is rearranged into a new size distribution that can be approximated by two components: a power-law distribution of small grains and a lognormal distribution of grains having the same size range as the initial distribution. Our results show that grain–grain collisions and sputtering are synergistic and that grain–grain collisions can play a crucial role in determining the surviving dust budget in supernova remnants.
We study intermittency of circulation moments in turbulent superfluid helium by using experimental grid turbulence and numerical simulations of the Hall-Vinen-Bekarevich-Khalatnikov model. More precisely, we compute the velocity circulation Γ r in loops of size r laying in the inertial range. For both experimental and numerical data, the circulation variance shows a clear Kolmogorov scaling < $Γ$$^2_r$ > ~ r 8/3 in the inertial range, independently of the temperature. Scaling exponents of high-order moments are comparable, within error bars, to previously reported anomalous circulation exponents in classical turbulence and low-temperature quantum turbulence numerical simulations.
Not Available
With high-penetration levels of renewable generating sources being integrated into the existing electric power grid, conventional generators are being replaced and grid inertial response is deteriorating. This technical challenge is more severe with photovoltaic (PV) generation than with wind generation because PV generation systems cannot provide inertial response unless special countermeasures are adopted. To enhance the inertial response, this paper proposes to synthesize a virtual inertia emulator (VIE) by using a battery energy storage system (BESS) and a three-phase grid-tied inverter to simulate a traditional generator such that the inertial response can be appropriately enhanced without sacrificing energy efficiency. Control systems for the VIE are presented in this paper, along with simulation results from PSCAD and MATLAB to validate the effectiveness of the proposed scheme.
This paper develops a hierarchical control frame-work to aggregate and to manage behind-the-meter distributed energy resources (DERs), which will be ubiquitous in future distribution systems. In the proposed framework, firstly, each controller in the hierarchy determines the flexibility of the DERs such that the obtained flexibility is feasible with respect to its operational purview. For example, the operational purview of a home energy management system may only consider consumer comfort preferences, while that for an aggregator or a grid controller may consider network voltage management as well. Based on the feasible flexibility, optimal setpoints for the DERs is then determined by the hierarchical controllers to help the distribution power network in voltage regulation, coordination issues with existing transmission-level conventional generators, etc. Therefore, the proposed strategy, which is based on model-predictive control, can be effectively utilized by the distribution network to coordinate several DERs to provide grid regulation services. Numerical simulations performed on the IEEE 37-bus test system demonstrate the efficacy of the proposed approach.
As the coordination of distributed energy re-sources becomes more necessary to provide grid services through a distributed energy resource management system (DERMS), studying the communication requirements for suc-cessful real-world, low-cost implementation becomes increasingly important. Specifically, this paper studies the necessary communication time resolutions between the different system components for a primal-dual-based DERMS, a highly de-veloped DERMS framework. We design a metric to evaluate the functionality of a DERMS with respect to providing grid services. Using numerical simulations based on a real-world feeder in Colorado, we show that the upper bound on the time resolutions is on the order of minutes instead of the previously assumed order of seconds.
Simulation of modern distribution system powerflow increasingly requires capturing the impact of advanced PV inverter voltage regulation on powerflow. With Volt/var control, the inverter adjusts its reactive power flow as a function of the point of common coupling (PCC) voltage. Similarly, Volt/watt control curtails active power production as a function of PCC voltage. However, with larger systems and higher penetrations of PV, this active/reactive power flow itself can cause significant changes to the PCC voltage potentially introducing oscillations that slow the convergence of system simulations. Improper treatment of these advanced inverter functions could potentially lead to incorrect results. This paper explores a simple approach to speed such convergence by blending in the previous iteration's reactive power estimate to dampen these oscillations. Results with a single large (5MW) PV system and with multiple 500kW advanced inverters show dramatic improvements using this approach.
To help bridge the gap between traditional power system engineering instruction and emerging industry needs, a new set of classroom and research tools are needed. Real-time simulation tools emulating power system control room software present an opportunity to introduce students to the array of operational considerations, technical challenges, and decision-making associated with power system operations. The GridAPPS-D platform is proposed to support coursework and academic research as it provides an open-source platform for simulation, application development, and software integration. The GridAPPS-D platform, simulation capabilities, development environment, and interactive training are discussed in the context of lessons-learned from implementation for undergraduate and graduate students in the US and India. A series of potential GridAPPS-D supported academic capabilities are introduced.
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.