Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 217 records · Page 12

Validation of Hermes-3 turbulence simulations against the TCV-X21 diverted L-mode reference case

Electrostatic flux-driven turbulence simulations with the Hermes-3 code are performed in TCV L-mode conditions in forward and reversed toroidal field configurations, and compared to the TCV-X21 reference dataset (Oliveira et al 2022 Nucl. Fusion 62 096001) qualitatively and with a quantitative methodology. Using only the magnetic equilibrium, total power across the separatrix (120 kW) and total particle flux to the targets (3 x 10 21 s −1 ) as inputs, the simulations produce time-averaged plasma profiles in good agreement with experiment. Shifts in the target peak location when the toroidal field direction is reversed are reproduced in simulation, including the experimentally observed splitting of the outer strike point into two density peaks. The overall normalized discrepancy between simulation and observation is better than any previously reported in the reversed field configuration, and matches the best previously reported in forward field configuration. Differences between simulation and experiment include density profiles inside the separatrix and at the inner target in forward (favorable $\bigtriangledown B$) field configuration. These differences in target temperature in forward field configuration lead to differences in the balance of current to the inner and outer divertor in the private flux region. The cause of these differences is most likely the lack of neutral gas in these simulations, indicating that even in low recycling regimes neutral gas plays an important role in determining edge plasma profiles. These conclusions are consistent with findings in Oliveira et al (2022 https://github.com/SPCData/TCV-X21).

Physics - Plasma physics↗

Dynamic Equivalence of Large-Scale Power Systems Based on Boundary Measurements

Parallel computing helped speed up many tasks that can be done independently. The advance in gaming industry did not produce a technology that helps the power industry to reduce the computation time for dynamic simulation in individual cases. To perform faster than real-time simulation for the purpose of predicting power system dynamic trajectory, power industry continues to struggle to reduce the system size and simulation time of large-scale power systems while keeping its dynamic behavior under various disturbances. There is also an acute need for real-time dynamic model reduction as more renewables enter the generation mix with dramatic changes in the generation outputs. All existing model reduction solutions are based on having access to a detailed dynamic model of the system. The system changes by the minutes, dynamic models are only updated annually. To solve this problem, wide-area measurements obtained by the phasor measurement unit (PMU) at the boundaries between the reduced system and the study system is used to represent the external area. An artificial neuro-fuzzy inference system (ANFIS) is established to perform the mapping of the measurement to the external equivalent model. Here the external area is regarded as a black-box. Model reduction studies are conducted on the Northeast Power Coordinating Council (NPCC) under various types of contingencies, by using a co-simulation approach between PSS/E and MATLAB. The results look very promising and will be discussed in this paper.

Tong, Ning↗

Evaluation of coupled wind-wave model simulations of offshore winds in the Mid-Atlantic Bight using lidar-equipped buoys.

From 2014 to 2017, two Department of Energy buoys equipped with Doppler lidar were deployed off the U.S. East Coast to provide long term measurements of hub-height wind speed in the marine environment. In this study, we performed simulations of selected cases from the deployment using a 5-km configuration of the Weather Research and Forecasting (WRF) model, to see if simulated hub height speeds could produce closer agreement with the observations than existing reanalysis products. For each case we performed two additional simulations: one in which marine surface roughness height was one-way coupled to forecast wave parameters from a standalone WaveWatch III (WW3) simulation, and another in which WRF and WW3 were two-way coupled using the Coupled-Ocean-Atmosphere-Wave-Sediment-Transport (COAWST) framework. It was found that all the 5-km WRF simulations improved 90-m wind speed statistics for the tropical cyclone case of 08 May 2015 and the cold frontal case of 25 Mar 2016, but not the nor-easter of 18 Jan 2016. The impact of wave coupling on buoy-level (4 m) wind speed was modest and case dependent, but when present, the impact was typically seen at 90 m as well, being as large as 10% in stable conditions. One-way wave coupling consistently reduced wind speeds, improving biases for 25 Mar 2016 but worsening them for 08 May 2015. Two-way wave coupling mitigated these negative biases, improved wave field representation and statistics, and mostly improved 4-m wind field correlation coefficients, at least at the VA buoy, largely due to greater self-consistency between wind and wave fields.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulation of Vibrational Sum Frequency Generation Intensity for Non-Centrosymmetric Domains Interspersed in an Amorphous Matrix: A Case Study for Cellulose in Plant Cell Wall

Vibrational sum frequency generation (SFG) spectroscopy can specifically probe molecular species non-centrosymmetrically arranged in a centrosymmetric or isotropic medium. This capability has been extensively utilized to detect and study molecular species present at the two-dimensional (2D) interface at which the centrosymmetry or isotropy of bulk phases is naturally broken. The same principle has been demonstrated to be very effective for the selective detection of non-centrosymmetric crystalline nanodomains interspersed in three-dimensional (3D) amorphous phases. However, the full spectral interpretation of SFG features has been difficult due to the complexity associated with the theoretical calculation of SFG responses of such 3D systems. This paper describes a numerical method to predict the relative SFG intensities of non-centrosymmetric nanodomains in 3D systems as functions of their size and concentration as well as their assembly patterns, i.e., the distributions of tilt, azimuth, and rotation angles with respect to the lab coordinate. We applied the developed method to predict changes in the CH and OH stretch modes characteristic to crystalline cellulose microfibrils distributed with various orders, which are relevant to plant cell wall structures. As a result, the same algorithm can also be applied to any SFG-active nanodomains interspersed in 3D amorphous matrices.

36 MATERIALS SCIENCE↗

Process-oriented stochastic perturbations applied to the parametrization of turbulence and shallow convection for ensemble prediction

Weather forecasting nowadays often requires some estimation of uncertainties associated with the output of meteorological models, in order to better inform decision making, especially in the context of intense weather events. Ensemble prediction systems provide such information through sets of possible scenarios which are designed to represent various uncertainty sources, including model uncertainties. A wide variety of methods have been proposed to estimate model uncertainties, among which perturbation methods targeting uncertain processes are a promising research field. Here, we focus on the representation of small-scale variability by process-oriented perturbation schemes applied to two key physical processes, namely turbulence and shallow convection. The perturbations are applied to a single-column version of the convection-permitting AROME model, in three idealized boundary-layer cases. Large-eddy simulations (LESs) of the same cases serve as a reference for the subgrid variability that has to be represented, and the results are also compared to those given by the Stochastically Perturbed Parametrization Tendencies (SPPT) method, which is a method commonly used by weather forecast centres to represent model uncertainty. The spread produced by our process-oriented perturbations of turbulence and shallow convection does not represent all the small-scale variability implied by the LESs for temperature and humidity. However, it is of a similar order of magnitude for the wind, thanks to perturbations generated by the stochastic turbulence scheme. The dispersion is structurally different from what is obtained with SPPT. It is non-negligible in the lower levels, where SPPT perturbations are usually suppressed because of numerical instabilities, indicating a possible complementarity between the schemes.

54 ENVIRONMENTAL SCIENCES↗

Accurate and confident prediction of electron beam longitudinal properties using spectral virtual diagnostics

Abstract Longitudinal phase space (LPS) provides a critical information about electron beam dynamics for various scientific applications. For example, it can give insight into the high-brightness X-ray radiation from a free electron laser. Existing diagnostics are invasive, and often times cannot operate at the required resolution. In this work we present a machine learning-based Virtual Diagnostic (VD) tool to accurately predict the LPS for every shot using spectral information collected non-destructively from the radiation of relativistic electron beam. We demonstrate the tool’s accuracy for three different case studies with experimental or simulated data. For each case, we introduce a method to increase the confidence in the VD tool. We anticipate that spectral VD would improve the setup and understanding of experimental configurations at DOE’s user facilities as well as data sorting and analysis. The spectral VD can provide confident knowledge of the longitudinal bunch properties at the next generation of high-repetition rate linear accelerators while reducing the load on data storage, readout and streaming requirements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inverter-Based Operation of Maui: Electromagnetic Transient Simulations

As larger and larger power systems approach and reach 100% inverter-based resource (IBR) operation during some hours of the year, questions arise regarding the stability of such extremely-high IBR power systems, the potential need for grid-forming (GFM) inverter technology, and the potential need for synchronous condensers. Relatedly, questions also arise about the ability of conventional positive sequence power system modeling tools to capture high-IBR system dynamics. This presentation introduces electromagnetic transient (EMT) simulations in PSCAD of the near-future (year 2023) Maui power system at and near 100% IBR operation and compares those simulations to positive sequence (PSSE) simulations. The Maui PSCAD model is parallelized on 30 cores and includes the entire transmission system (>200 three-phase buses), >170 individual and aggregate IBR models, four wind plants, and three synchronous generators plus six synchronous condensers at two locations. We investigate system stability with varying levels of inertia using conventional grid-following IBR controls, and then we investigate the impact of GFM controls on stability. Results suggest that: 1) positive sequence simulations can miss key dynamics in extremely high IBR cases; 2) EMT simulations can also miss key dynamics if inverter inner control dynamics are not modeled; 3) synchronous condensers can stabilize a system in which 100% of the energy is supplied by IBRs, even conventional grid-following IBRs; 4) GFM controls on just some of the IBRs can stabilize a 100% IBR power system, even if that system has zero inertia (i.e. no synchronous condensers, though synchronous condensers may be needed for other purposes such as protection system operations).

41 EE - Solar Energy Technologies Office (EE-4S)↗

2D radial-azimuthal particle-in-cell benchmark for E × B discharges

In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural uncertainty assessment for low-Mach wall-resolved large-eddy simulations: Plane Channel and Periodic Hill Use Cases

A structural uncertainty assessment is provided for a set of wall-resolved large-eddy simulations (WRLES) using a turbulent channel flow (𝑅𝑒 𝜏 395) and a periodic hill in the turbulent regime (𝑅𝑒 𝑏 2800) validation suite. The standard one-equation subgrid scale turbulent kinetic energy model (𝑘 𝑠𝑔𝑠 ) using both a static and dynamic coefficient approach along with Wall Adapting Local Eddy Viscosity (WALE) are exercised. The 𝑘 𝑠𝑔𝑠 model activated in this study also allows for a correction that supports theoretical turbulent viscosity and total kinetic energy scaling (cubit and quadratic, respectively) in the near-wall regime. Channel and hill results showcase very good agreement between prediction and simulation for the WALE and the dynamic coefficient 𝑘 𝑠𝑔𝑠 model. The study also demonstrates that near-wall damping improves predictivity for the static coefficient 𝑘 𝑠𝑔𝑠 model, while proving an unnecessary additional contribution for the dynamic coefficient 𝑘 𝑠𝑔𝑠 model as near-wall scaling is naturally captured. Overall conclusions are that this suite of WRLES models are well suited for production usage in the Sierra-based low-Mach codes.

97 MATHEMATICS AND COMPUTING↗

Challenges for density functional theory in simulating metal–metal singlet bonding: A case study of dimerized VO2

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.

Chemistry↗

Updraft dynamics and microphysics: on the added value of the cumulus thermal reference frame in simulations of aerosol–deep convection interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.

42 ENGINEERING↗

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Simulating current flow through a well casing and an induced fracture

A system includes a receiver component that receives a model of a well system, the model comprising a representation of an electrically energized well casing and/or a fracture as a transmission line that leaks electric current as the current traverses the well casing and/or fracture. The receiver component also receives a value that indicates an amount of the current that is applied to the well casing and a location of a source of the current on the well casing. The system further includes an electromagnetic field calculator component that calculates an estimated electromagnetic field at at least one location on the surface of the earth based at least in part upon the representation of the electrically energized well casing and/or fracture, the value that indicates the amount of current that is applied to the well casing, and the location of the source of the current on the well casing.

02 PETROLEUM↗

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗

Analysis of the neutral fluxes in the divertor region of Wendelstein 7-X under attached and detached conditions using EMC3-EIRENE

Abstract This paper analyzes the neutral fluxes in the divertor region of the W7-X standard configuration for different input powers, both under attached and detached conditions. The performed analysis is conducted through EMC3-EIRENE simulations. They show the importance of the horizontal divertor to generate neutrals, and resolve the neutral plugging in the divertor region. Simulations of detached cases show a decrease in the number of generated neutrals compared to the attached simulations, in addition to a higher fraction of the ion flux arriving on the baffles during detachment. As the ionization takes place further inside the plasma during detachment, a larger percentage of the generated neutral particles leave the divertor as neutrals. The leakage in the poloidal and toroidal direction increases, just as the fraction of collected particles at the pumping gap. The fraction of pumped particles increases with a factor two, but stays below one percent. This demonstrates that detachment with the current target geometry, although it improves the power exhaust, is not yet leading to an increased particle exhaust.

Physics↗

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z dr and larger specific differential phase K dp above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z dr and K dp observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Z dr and K dp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

54 ENVIRONMENTAL SCIENCES↗

BESTEST-GSR (Building Energy Simulation Test - Generation Simulation and Reporting) 2023 [SWR 18-23]

Building Energy Simulation Test (BESTEST) is an NREL-developed method to validate the qualitative performance of different whole building simulations engines relative to each other. https://www.nrel.gov/docs/legosti/old/6231.pdf The purpose of this repository is to generate BESTEST test case models, run simulations, and populate data for ASHRAE Standard 140 reporting spreadsheets for EnergyPlus® based whole building simulation tools. It was originally setup for 2014 version of Standard 140. In May of 2022 it was updated to the 2020 version of Standard 140. This update included updates and additions to existing test suites, the bulk of which was in Section 5.2 (Building Thermal Envelope and Fabric Load Tests). We did not add Section 5.5 Airside HVAC Equipment Performance, but we hope to add that later in 2022. At some point we also hope to add Section 5.2.4 ground modeling, which is currently excluded. Supported Tools The default IDF generation is based on the OpenStudio® CLI, but the workflow supports a 'Bring your own IDF' use case. Additionally, for non-EnergyPlus® based tools the post processing scripts can be used if simulation results are provided as a CSV file.The scripts on this repository should work on Mac, Windows, and Linux. Dependencies Install OpenStudio® 3.4.0 make sure command line can recognize the 'openstudio' command This includes EnergyPlus® 22.1 Install Ruby on your system if it isn't already setup. 2.7 is used for development but other versions may work Since OpenStudio has its own embedded Ruby, which is used for running measures, you don't necessarily have to use a version of Ruby supported by OpenStudio. Install RubyXL Ruby gem This is used to modify Microsoft Excel spreadsheets Install Parallel Ruby gem This allows the CLI to run simulations in parallel

Goldwasser, David↗