Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional expectation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Forest reorganisation effects on fuel moisture content can exceed changes due to climate warming in wet temperate forests

Abstract The distributions of vegetation and fire activity are changing rapidly in response to climate warming. In many regions, climate effects on dead fuel moisture content (FMC) are expected to increase future wildfire activity. However, forest FMC is largely driven by microclimate conditions, which are moderated from open weather by vegetation canopies. As shifts in vegetation increase under climate warming, the extent to which future fire activity will be driven by climate directly or associated vegetation shifts remains unresolved. Here, we present a study aimed at quantifying the relative magnitudes of (i) direct climate warming, and (ii) vegetation change, on FMC. Field sites to evaluate these effects were established in a natural laboratory of altered forest states to mature wet temperate forest in south‐eastern Australia. FMC was estimated using a process‐based model and 48 years of reconstructed climate data. Canopy effects on microclimate were captured by transferring inputs from climate to microclimate using models parameterised with field observations. To evaluate the relative magnitude of climate and vegetation effects, we calculated the maximum difference in mean annual FMC across annual climate replicates and compared this to FMC differences across reorganising forest sites. Our results show vegetation effects on FMC can exceed those related to expected climate change. Changes to forest structure and composition increased (+15.7%) and decreased (−12.3%) mean annual FMC, with a larger negative effect when forest cover was completely removed (−18.5%). In contrast, the largest climate effect on FMC was −6.6% across 48‐years of data. Our study demonstrates that the magnitude of vegetation effects on FMC can exceed expected climate change effects. Models of future fire activity that do not account for changing vegetation effects on microclimate are omitting a key biophysical control on FMC and therefore may not be accurately predicting future fire activity.

Brown, Tegan P.↗
DOE OSTI2021

Atomic Force Microscopy and Infrared Nanospectroscopy of COVID-19 Spike Protein for the Quantification of Adhesion to Common Surfaces

The COVID-19 pandemic has claimed millions of lives worldwide, sickened many more, and has resulted in severe socio-economic consequences. As society returns to normal, understanding the spread and persistence of SARS CoV-2 on commonplace surfaces can help to mitigate future outbreaks of coronaviruses and other pathogens. We hypothesize that such an understanding can be aided from studying the binding and interaction of viral proteins with non-biological surfaces. Here, we propose a methodology for investigating the adhesion of the SARS CoV-2 spike glycoprotein on common inorganic surfaces such as aluminum-, copper-, iron-, silica-, and ceria-oxides as well as metallic gold. Quantitative adhesion was obtained from the analysis of measured forces at the nanoscale using an atomic force microscope operated under ambient conditions. Without imposing further constraints on the measurement conditions, our preliminary findings suggest that spike glycoproteins interact with similar adhesion forces across the majority of the metal oxides tested with the exception to gold, for which attraction forces ~10 times stronger than all other materials studied were observed. Ferritin, which was used as a reference protein, was found to exhibit similar adhesion forces as SARS CoV-2 spike protein. Further, this study results show that glycoprotein adhesion forces for similar ambient humidity, tip shape, and contact surface are non-specific to the properties of metal oxide surfaces, which are expected to be covered by a thin water film. The findings suggest that at ambient conditions glycoprotein adhesion to metal oxides is primarily controlled by the water capillary forces and they depend on the surface tension of the liquid (water). We discuss further strategies warranted to decipher the intricate nanoscale forces for improved quantification of the adhesion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2023

RPC based tracking system at CERN GIF++ facility

With the HL-LHC upgrade of the LHC machine, an increase of the instantaneous luminosity by a factor of five is expected and the current detection systems need to be validated for such working conditions to ensure stable data taking. At the CERN Gamma Irradiation Facility (GIF++) many muon detectors undergo such studies, but the high gamma background can pose a challenge to the muon trigger system which is exposed to many fake hits from the gamma background. A tracking system using RPCs is implemented to clean the fake hits, taking profit of the high muon efficiency of these chambers. This work will present the tracking system configuration, used detector analysis algorithm and results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Wildfire Smoke Demonstrates Significant and Predictable Black Carbon Light Absorption Enhancements

Abstract Black carbon (BC) is estimated to have the second largest anthropogenic radiative forcing in earth‐systems models (ESMs), but there is significant uncertainty in its impact due to complex mixing with organics. Laboratory‐generated particles show that co‐mixed non‐absorbing material enhances absorption by BC by a factor of 2–3.5 as predicted by optical models. However, weak or no enhancements are often reported for field studies. The cause of lower‐than‐expected absorption is not well understood and implies a lower radiative impact of BC compared to how many ESMs currently treat aerosols. By analyzing BC aerosol particle‐by‐particle we reconcile observed and expected absorption for ambient smoke plumes varying in geographic origin, fuel types, burn conditions, atmospheric age and transport. Although particle‐by‐particle tracking is computationally prohibitive for sophisticated ESMs we show that realistic BC absorption is reliably estimated by bulk properties of the plume providing a suitable parameterization to constrain black carbon radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A shared-mobility-based framework for evacuation planning and operations under forecast uncertainty

To meet evacuation needs from carless populations who need personalized assistance to evacuate safely, in this article we propose a ridesharing-based evacuation program that recruits volunteer drivers before a disaster strikes, and then matches volunteer drivers with evacuees once demand is realized. Here we optimize resource planning and evacuation operations under uncertain spatiotemporal demand, and construct a two-stage stochastic mixed-integer program to ensure high demand fulfillment rates. We consider three formulations to improve the number of evacuees served, by minimizing an expected penalty cost, imposing a probabilistic constraint, and enforcing a constraint on the conditional value at risk of the total number of unserved evacuees, respectively. We discuss the benefits and disadvantages of the different risk measures used in the three formulations, given certain carless population sizes and the variety of evacuation modes available. We also develop a heuristic approach to provide quick, dynamic and conservative solutions. We demonstrate the performance of our approaches using five different networks of varying sizes based on regions of Charleston County, South Carolina, an area that experienced a mandatory evacuation order during Hurricane Florence, and utilize real demographic data and hourly traffic count data to estimate the demand distribution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Isolation and conversion of electrolyte components into a value-added product

The overall objective of this work was to develop a process for effectively (i) recycling the electrolyte components from the end-of-life batteries and reduce the environmental hazard, (ii) understand the fundamental differences between pristine and recycled electrolyte components, (iii) electrochemically and analytically evaluate the recycled electrolyte components and (iv) establish deviations from the pristine electrolyte, and effectively separation of value-added product. We focused on determining scaling up the electrolyte recovery process by distillation of end-of-life batteries to get yield enough to use it for electrolyte formulation of reuse in the battery cell. The used cells from the ORNL battery manufacture facility (BMF)as well as received from the industry partner, Austin Elements Inc., was scaled to 10 used multilayered pouch cells in each batch of distillation and systematically separated and identified the electrolyte components of the used cells with FTIR and NMR studies. We were able to recovery 10g of solvent mass from the used pouch cells. NMR spectroscopy showed that the solvent was pure EMC,DMC and the sample was used to make a new battery cell. The end-of-life battery, it was noted that some cells were visually drier than the other cells and were more advanced in aged. In that case, even though the same distillation conditions were utilized, no usable mass of electrolyte solvent was recovered. It is expected that most of the solvent was decomposed during the battery charge-discharge cycling. The solid salt residue was separated by chemical and water treatment processes. Components of resulting solid product was analyzed by XRD and XPS measurements. The solid products are very much related to the salts used in the electrolyte solution.

25 ENERGY STORAGE↗
DOE OSTI2025

Isolation and Conversion of Electrolyte Components into a Value-Added Product

The overall objective of this work was to develop a process for effectively (i) recycling the electrolyte components from the end-of-life batteries and reduce the environmental hazard, (ii) understand the fundamental differences between pristine and recycled electrolyte components, (iii) electrochemically and analytically evaluate the recycled electrolyte components and (iv) establish deviations from the pristine electrolyte, and effectively separation of value-added product. We focused on determining scaling up the electrolyte recovery process by distillation of end-of-life batteries to get yield enough to use it for electrolyte formulation of reuse in the battery cell. The used cells from the ORNL battery manufacture facility (BMF) as well as received from the industry partner, Austin Elements Inc., was scaled to 10 used multilayered pouch cells in each batch of distillation and systematically separated and identified the electrolyte components of the used cells with FTIR and NMR studies. We were able to recovery 10g of solvent mass from the used pouch cells. NMR spectroscopy showed that the solvent was pure EMC, DMC and the sample was used to make a new battery cell. The end-of-life battery, it was noted that some cells were visually drier than the other cells and were more advanced in aged. In that case, even though the same distillation conditions were utilized, no usable mass of electrolyte solvent was recovered. It is expected that most of the solvent was decomposed during the battery charge-discharge cycling. The solid salt residue was separated by chemical and water treatment processes. Components of resulting solid product was analyzed by XRD and XPS measurements. The solid products are very much related to the salts used in the electrolyte solution.

25 ENERGY STORAGE↗
DOE OSTI2023

Sensitivity study on the fin effect of thermocouple mounted on the heated surface under film boiling condition

In the reactivity-initiated accident (RIA) scenario, right after the power burst with the occurrence of departure from nucleate boiling (DNB), the cladding temperature quickly reaches its maximum and is maintained in a steady film boiling regime. Then, under limited heat transfer performance of film boiling, the local cooling effect by the thermocouple wire may not be neglected. A computational model is developed to estimate the local cooling effect of the thermocouple wire mounted on the cladding surface under film boiling condition. The model predicts the temperature profile along the cladding surface as physically expected and is verified against the benchmark problem. Sensitivity studies reveal the significant effects of the thickness of vapor film layer and the diameter of thermocouple on the temperature drop near the thermocouple, for they increase the heat transfer area to the liquid water region. The impact of film boiling heat transfer coefficient, however, is not significant, as the primary heat transfer mechanism near the thermocouple is the heat conduction to the wire. Future study will include the contribution of radiative heat transfer between the cladding surface and the vapor film layer, as well as a transient power history which represents an actual RIA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Influence of Phase-Change Materials and Operating Periods on Load Shifting for Thermal Storage Heat Pumps in European Climate Zones

As the reliance on heat pumps (HPs) for space conditioning increases worldwide, ex-tensive grid strains during utility peak hours are expected. Significant research efforts focus on investigating the use of integrated thermal energy storage (TES) systems to reduce peak energy demands. Conventionally, TES systems are often integrated to store and provide either cooling or heating, but not both. In this work, a single room-temperature phase-change material (PCM)-TES was integrated for dual-mode opera-tion with an R290 air-to-water 10.5 kW (3-Ton) HP. Two PCM melting temperatures (17°C, 22°C) were assessed to investigate their impact on system performance and cost reductions for Milan, Italy (4A); Barcelona, Spain (3A); and Oslo, Norway (6A). HP and HP-TES performance maps were generated using Modelica-based transient models and co-simulated with a prototype residential building using the Spawn of En-ergyPlusTM with simple rule-based controls. Our simulations indicate that PCM-17°C provided annual cost savings at higher peak-to-off-peak cost ratios, as the demand re-duction dominated the savings potential. The recharge energy costs and requirements became more dominant at lower cost ratios, making PCM-22°C more beneficial at lower cost ratios. Finally, despite Oslo's heating-dominated climate, the utility peak hours favor cooling demand reductions, making PCM-17°C more appealing. In all cases, the cost savings potential stabilized between 4 – 10% regardless of cost ratio, suggesting that system-level component and operation control optimization is required to maximize cost savings further.

25 ENERGY STORAGE↗
DOE OSTI2026

Enhancing the mass resolving power of FRIB’s proposed high-voltage MR-ToF mass separator and spectrometer: Addressing non-ideal conditions

Multi-reflection time-of-flight mass separators and spectrometers (MR-ToF MSs) are indispensable tools at radioactive ion beam (RIB) facilities. These electrostatic ion beam traps act as highly selective mass separators and high-precision mass spectrometers for rare and exotic nuclei. When well-tuned and designed to minimize higher-order flight-time aberrations, state-of-the-art MR-ToF MSs approach, and slightly exceed, mass resolving powers of m/Δm = 10 6 . Achieving m/Δm ≥ 3 ⋅ 10 6 would provide the ability to resolve > 90 % of all known isomeric states with half-lives above 10 ms. However, the ability to mass separate in all practical setups is limited by non-ideal conditions which place such resolving powers out of reach. Here, to this end, we present a simulated analysis of these conditions in the newly proposed high-voltage MR-ToF MS for the Facility for Rare Isotope Beams (FRIB). It is expected to store ions at 30 keV beam energy and increase ion throughput by two orders of magnitude compared to current devices. Existing efforts to mitigate the effects of non-ideal conditions employed for current MR-ToF devices storing ions at < 3 keV beam energy will already enable mass resolving powers approaching 10 6 for FRIB’s high-voltage MR-ToF device. Simulations of newly proposed mitigation strategies show that even mass resolving powers approaching 10 7 might become feasible.

Electrostatic ion beam traps↗
DOE OSTI2025

Anatomy of a Turtle Confiscation: What to Expect When You're Not Expecting and Recommendations for Improving Outcomes

Turtles are being subjected to unprecedented collection for the illegal wildlife trade, but only a portion of the trade is detected. When confiscations do happen, wildlife agencies must decide what to do with those animals—euthanize them, place them in permanent captivity, or release them back into the wild. Here, we present a case study of a confiscation of .200 eastern box turtles (Terrapene carolina carolina) and our efforts to repatriate them back to the wild. Twentyfive percent of turtles died in captivity, and at least another 33% died during the soft-release phase. Approximately half of the confiscated turtles survived until they were released from their soft-release pen 9 months post-confiscation. For each phase, from the time turtles were seized until released, we describe our objectives, the challenges we encountered, and our recommendations for improving future turtle confiscations. Given the extended stressful conditions that confiscated turtles often experience before being seized, it is important to recalibrate our expectations regarding future confiscation outcomes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Modeling Air Handling Units to Create a Diverse Fault Dataset for FDD Innovation: Lessons Learned and Recommendations

As energy management and information systems (e.g., automated fault detection and diagnostics [AFDD] tools) become more prevalent in the commercial building stock, it is important to determine the effectiveness of these technologies by benchmarking their performance. The authors have been working to develop the largest publicly available dataset of HVAC fault datasets for performance benchmarking applications, covering the most common HVAC systems and designs including chiller plants, rooftop packaged units, dual duct air handling unit and single duct air handling units. This study covers the development, modeling, and validation of a synthetic fault dataset for the air handling unit (AHU), one of the most common HVAC configurations found in the commercial building stock. Despite this being a common system, real-world time series data are scarce and usually do not span a wide range of weather conditions. Due to this limitation, two detailed AHU models, which included the single duct AHU and dual duct AHU developed in the Modelica language and HVACSIM+ were employed to carry out annual simulations of numerous common sensor faults, mechanical faults, and control sequence faults. The fault inclusive data were then validated by comparing fault effects on system performance to expected symptoms. We summarize the nature of each fault and their impacts under different weather and operation conditions. We report some lessons learnt during the efforts of validating the high volumes of the FDD data sets. Finally, we highlight considerations for FDD developers that may want to use this dataset to assess their algorithms’ performance and their improvement over time.

Casillas, Armando↗
DOE OSTI2022

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10 15 W/cm 2 ), counterstreaming plasma flows (V≲ 2000 km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ~5–6 linear growth times into the evolution. Our new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (~300–500 ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (~300 μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ~10–30 T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Correcting the PFNS for more consistent fission modeling

For FY20, we had a deliverable to write a report detailing efforts to simultaneously evaluate both the prompt neutron multiplicity, $\overline{ν}$, and the prompt neutron fission neutron spectrum, PFNS, using CGMF. CGMF is the LANL-developed fission fragment decay code that consistently evaporates prompt neutrons and γ rays using the Hauser-Feshbach statistical theory of compound reactions. The decay begins by constructing the initial conditions of the fission fragments, then decaying each one from the excited state by neutrons and γ rays, conserving energy, momentum, spin, and parity in each step of the emission. The initial conditions of the fragments, along with the multiplicity, energy, and direction of each emitted neutron or γ ray, are recorded, allowing for the full reconstruction of the fission event. These event histories allow us to reconstruct average quantities, as well as correlations between observables, that can be compared with experimental or evaluated data. In that initial report, although there was already a favorable comparison between $\overline{ν}$ from CGMF, experiment, and the current ENDF/B-VIII.0 evaluation, we showed that there was still significant work to be done to improve the PFNS from CGMF. Historically, the PFNS is calculated too soft by Hauser-Feshbach fission models, and CGMF is no exception. The incorrect shape presents a significant challenge in fission modeling, including for our understanding of the fission process and for our ability to consistently calculate and predict a variety of prompt fission observables (such as fission fragment initial conditions, neutron and γ-ray multiplicities and energies, and the correlations between all observables). In our companion report, we detail our success in using CGMF to evaluate $\overline{ν}$. Although not included in the optimization explicitly, we also keep the initial conditions of the fission fragments physical, along with reproducing reasonably well the neutron multiplicity distribution. As we would expect from the sensitivities calculations from, the average neutron energies change very little from the $\overline{ν}$ optimization along with the PFNS (as will be shown in Sec. 2.6). The conclusion was that the global and statistical models would have to be investigated instead of just the fission fragment initial conditions (as is sufficient for $\overline{ν}$). This report details those efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗