Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Predictive performance of multi-model ensemble forecasts of COVID-19 across European nations

Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022. Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1–4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models’ predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models’ forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models’ past predictive performance. Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models’ forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models’ forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models’ forecasts of deaths (N=763 predictions from 20 models). Across a 1–4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models. Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks. Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 ( https://www.nfdi4health.de/task-force-covid-19-2 ) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. Additionally, we compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Human reliability analysis studies from simulator experiments using Bayesian inference

Probabilistic Safety Assessment (PSA) of complex facilities is performed to arrive at the risk posed by them. PSA also accounts for the contribution of the human errors towards the overall risk through Human Reliability Analysis (HRA) in terms of Human Error Probability (HEP). Human operators are part of the system and do not work in isolation. Their performance is influenced by the context in which the actions are performed. As a result, quantification of HEP requires operator performance data under the given context. Some good sources of operator performance data are plant‘s operation data, simulator data and expert judgement. The plant operation data pertaining to HRA is generally sparse. In this situation, a full scope plant simulator provides a good alternative for operator performance data generation. Many of the currently practiced HRA methods have been developed by combining the empirical evidence with expert judgement and contain a lot of uncertainty in their estimates. Bayesian inference is suitable for updating the prior HRA estimates with the simulator evidence to obtain the posterior HEP. Here, posterior HEP has been calculated for postulated accident scenarios in advanced reactor (first of its kind) at design stage, using plant simulator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The miniJPAS survey: AGN and host galaxy coevolution of X-ray-selected sources

Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass (MBH) and the stellar mass of their host galaxies (M * ). They even show similar histories across cosmic times of their differential terms: the star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We selected an X-ray sample of active galactic nuclei (AGN) up to z = 2.5 in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with the CIGALE code, constraining the emission to 68 bands, from which 54 are the narrow filters from the miniJPAS survey. For a final sample of 308 galaxies, we derived their physical properties, such as their M * , SFR, star formation history (SFH), and the luminosity produced by the accretion process of the central BH (L AGN ). For a subsample of 113 sources, we also fit their optical spectra to obtain the gas velocity dispersion from the broad emission lines and estimated the M BH . We calculated the BHAR in physical units depending on two radiative efficiency regimes. We find that the Eddington ratios (λ Edd ) and its popular proxy (L X /M * ) have a difference of 0.6 dex, on average, and a KS test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the MBH - M * scaling relation, which can explain the difference between λ Edd and its proxy. We also modeled three evolution scenarios for each source to recover the integral properties at z = 0. Using the SFR and BHAR, we show a notable diminution in the scattering between MBH - M * . For the last scenario, we considered the SFH and a simple energy budget for the AGN accretion, and we retrieved a relation similar to the calibrations known for the local Universe. Our study covers ~1 deg 2 in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between the differential values (SFR and BHAR) and their decoupling based on an energy limit is the key that leads to the local M BH - M * scaling relation. In the future, we plan to extend this methodology to a thousand degrees of the sky using JPAS with an X-ray selection from eROSITA, to obtain an unbiased distribution of BHAR and Eddington ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Signatures of linearized gravity in atom interferometers: A simplified computational framework

We develop a general framework for calculating the leading-order, general relativistic contributions to the gravitational phase shift in single-photon atom interferometers within the context of linearized gravity. We show that the atom gradiometer observable, which only depends on the atom interferometer propagation phase, can be written in terms of three distinct contributions: the Doppler phase shift, which accounts for the tidal displacement of atoms along the baseline, the Shapiro phase shift, which accounts for the delay in the arrival time of photons at atom-light interaction points, and the Einstein phase shift, which accounts for the gravitational redshift measured by the atoms. For specific atom gradiometer configurations, we derive the signal and response functions for two physically motivated scenarios: (i) transient gravitational waves in the transverse-traceless gauge and, for the first time, in the proper detector frame, and (ii) transient massive objects sourcing weak and slow-varying Newtonian potentials. Here, we find that the Doppler contribution of realistic Newtonian noise sources (e.g., a freight truck or a piece of space debris) at proposed atom gradiometer experiments, such as AION, MAGIS, and AEDGE, can exceed the shot noise level and thus affect physics searches if not properly subtracted. Note: 44 pages including appendices, 4 figures; v2: different formatting, fixed typos, 28 pages including appendices, 5 figures, agrees with published version

Atom interferometry↗
DOE OSTI2020

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗
DOE OSTI2022

Application of SCALE to Molten Salt Fueled Reactor Physics in Support of Severe Accident Analyses

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a molten salt-fueled reactor (MSR). SCALE capabilities for the modeling of MSR physics were demonstrated based on the Molten Salt Reactor Experiment (MSRE). Of primary interest were the determination of the system’s nuclide inventory, as well as the inventories in the various regions of the loop, considering that the fuel is continuously pumped through the system. This report contains discussions on the following: 1. Determination of the system-average fuel salt inventory considering fission gas removal in the off-gas system and noble metal removal through plating out at the heat exchanger using recent enhancements in SCALE’s depletion sequence TRITON, 2. Assessment of the nuclide spatial distribution throughout the loop using SCALE’s depletion solver ORIGEN, 3. Calculation of the core’s power profile, flux profile, temperature reactivity coefficients, and xenon reactivity using full-core calculations with SCALE’s Monte Carlo code KENO-VI. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the system, as well as the inventory and decay heat of individual regions in the loop, a zone-wise power profile, temperature feedback coefficients, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗
DOE OSTI2020

Excitation of Bloch Surface Waves in Zero-Admittance Multilayers for High-Sensitivity Sensor Applications

The growing need for classical as well as quantum optical sensing places increasingly stringent requirements upon the desired characteristics of the engendered fields. Specifically, achieving superior field enhancement plays a critical role in applications ranging from chem-bio sensing, Raman and infrared spectroscopies to ion trapping and qubit control in emerging quantum-information science. Due to their low optical losses and ability to exhibit resonant field enhancements, all dielectric multilayers are emerging as an optical material system not only useful to classical photonics and sensing but also of potential to be integrated with quantum materials and quantum sensing. The recently introduced concept of zero-admittance layers [1] within dielectric multilayer materials, enables the creation and control of resonant fields orders of magnitude larger than the exciting field. Here, invoking the zero-admittance concept, we design, fabricate, and characterize an all-dielectric nonabsorbing stack and demonstrate the engendered huge field enhancement. Describing the fields in terms of Bloch surface waves, we connect the surface field to the semiperiodicity in the dielectric domains of the stack. As a specific application of the resonant field, we propose and demonstrate refractive-index sensing for the detection of trace amounts of an analyte. The results include a quantification of the sensitivity of the device with respect to the profile of the exciting field. The experimental results are shown to be in good agreement with theoretical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2025

Rhenium-Selenido Corroles: Reflections on 5d Metalloporphyrins and Metallocorroles as Triplet Emitters and Photosensitizers

Following the successful synthesis of rhenium-oxido, rhenium-imido, and rhenium-sulfido corroles, a series of para-X-substituted rhenium-selenido triarylcorroles (X = OCH 3 , CH 3 , H, F, CF 3 ) have been prepared in 70–76% yields. A one-pot, two-step procedure was used, analogous to that used for ReS corroles, involving rhenium insertion via high-temperature reaction of Re 2 (CO) 10 and free-base corroles in refluxing decalin, followed by exposure to PCl 3 (which is thought to deoxygenate ReO corroles formed under the reaction) and elemental selenium. An analogous procedure, however, failed to yield rhenium-tellurido corroles, presumably reflecting, according to DFT calculations, the limited stability of these species. Unlike ReO corroles, ReS and ReSe corroles were found not to exhibit phosphorescence in the NIR-I regime (600–1000 nm); nor did they sensitize singlet oxygen formation. In the hope of obtaining a broader perspective of this negative result, we also examined ReN porphyrins and found them not to phosphoresce or to sensitize singlet oxygen formation. These results were explained by DFT and TDDFT calculations in terms of low-energy triplet states, which are not energetic enough to excite molecular oxygen to its lowest singlet state. Whether some of the new Re corroles exhibit phosphorescence in the NIR-II regime remains an interesting question for future investigation.

Aromatic compounds↗
DOE OSTI2025

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
DOE OSTI2026

Asynchronous Grid Connections Providing Fast-Frequency Response: System Integration Study

This paper presents an integration study for the recent power electronic-based fast-frequency response technology, "asynchronous grid connection" which operates as an aggregator for behind-the-meter resources and distributed generators. Both technical feasibility and techno-economic viability studies are presented. The fast-frequency response characteristics, validated against Power Hardware-in-the-Loop experiments, are integrated into an IEEE 9- bus system in DigSilent PowerFactory for system-level dynamic analysis. It demonstrates that droop-based control enhancements to local distributed generators allow their aggregation to provide grid-supporting functionalities and participate in the ancillary service markets. To this end, a long-term simulation embedding the system within the ancillary service market framework of PJM has been performed. The fast-frequency response regulation is subsequently used to calculate the potential revenue and project the results on a 15-year investment horizon. Finally, the techno-economic analysis provides recommendations for enhancements to access the full potential of distributed generators on a technical and regulatory level.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Nuclear Structure Studies for Schiff Moment Enhancement

One of the major open questions in modern physics is the origin of the matter-antimatter asymmetry observed in the universe. While the Standard Model contains sources of charge-parity (CP) violation, they appear insufficient to explain the observed dominance of matter over antimatter. As a result, searches for additional sources of CP violation are an important area of contemporary nuclear and particle physics. Atomic electric dipole moments (EDMs) provide sensitive probes of CP-violating interactions. In heavy nuclei, the observable atomic EDM can be enhanced through the presence of a nuclear Schiff moment, which arises from CP-violating nuclear forces and nuclear structure effects. The review paper Nuclear Schiff Moments and CP Violation by Jonathan Engel provides an overview of the theoretical framework connecting CP violation, nuclear structure, and Schiff moments. The long-term goal of this project is to contribute to the identification and characterization of nuclei that may exhibit enhanced Schiff moments. Such nuclei are often associated with strong octupole correlations or low-energy parity-doublet structures that can amplify CP-violating effects. This report summarizes progress made during the initial phase of the project, including a survey of candidate isotopes and preparation for nuclear-structure calculations using the HFODD code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Accuracy of HVAC Load Predictions: Validation of EnergyPlus and DOE-2 using FLEXLAB Measurements

The aim of the project reported here was to better understand the level of accuracy of three building energy simulation (BES) engines (‘engines’) — EnergyPlus™, DOE-2.1e, and DOE-2.2 — by identifying and investigating significant deviations between the performance predicted by these engines and actual performance as measured in the FLEXLAB® test facility at Lawrence Berkeley National Laboratory (LBNL). The specific test conditions included some of those prescribed in ANSI/ASHRAE Standard 140 - Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs. Detailed measurements of FLEXLAB performance, including indoor temperatures and heat fluxes and air-flow and water flow rates and temperatures in the Heating, Ventilating and Air Conditioning (HVAC) system, together with hourly weather data, were recorded and used in analyzing the simulation results from EnergyPlus v8.8, DOE-2.2 v3.65 and DOE-2.1e v127. These engines are commonly used in the United States for building energy code compliance, federal, state, and utility incentives programs, as well as energy efficient design of new buildings and energy retrofit of existing buildings. Seven conventional overhead mixing ventilation scenarios were tested and each engine was found to have a similar level of agreement with the measurements of space-level heating and sensible cooling loads. These results provide useful information regarding the accuracy of these engines in predicting the cooling and heating load elements of whole building energy performance. This information is intended for practitioners who are concerned about transitioning between simulation tools with different engines and for managers of utility programs leveraging these tools for evaluating and/or projecting measure savings to be incentivized under their programs. The results of the comparisons of simulated and measured performance indicate that the predictions from all three engines are not significantly different. The 24-hour average value of the absolute mean bias indicates the likely magnitude of the error in any particular case. The average mean bias is reduced by cancelation of overprediction in one case by underprediction in another. The daytime absolute mean biases, which may be more important for both energy performance and occupant comfort, are ~6%, presumably because of the greater complexity involved in simulating in the presence of solar radiation. EnergyPlus typically overpredicts the cooling load and/or underpredicts the heating load by ~1.5% and the DOE-2 engines typically underpredict the cooling load by approximately the same amount. The Root Mean Square Error is relatively more sensitive to shorter term variations in the difference between predicted and measured loads; the three engines have similar values, ~10%, suggesting that the uncertainties in their predictions of peak loads may also be similar in magnitude. The implication of these results is that users, both designers and program analysts, can use EnergyPlus, DOE-2.1e, or DOE-2.2 to model conventional commercial buildings equipped with overhead mixing ventilation with a similar level of confidence. Further work is required to better understand the variability in the level of agreement between the engine predictions and FLEXLAB measurements, where a particular engine will agree well with FLEXLAB in some cases and not so well in others and another engine will agree or disagree in different cases. As the sources of this variability are identified and eliminated or reduced significantly, it is recommended that the experimental capabilities and methods developed in the study reported here should be applied to validating heating and cooling load calculations for spaces with different types of furniture and miscellaneous loads. These methods should then be applied to low energy space conditioning systems in EnergyPlus including, in particular, radiant slab and radiant ceiling panel cooling and heating systems and ‘mixed mode’ systems that combine mechanical cooling and natural ventilation systems, focusing on controls, including control of thermal mass. The work reported here addresses the conventional method of heating and cooling occupied spaces; other methods, such as the use of radiant heating and cooling systems have the potential to provide equivalent occupant comfort, or better, with lower energy consumption. These systems are addressed more explicitly in EnergyPlus but there is a need for empirical validation to give users the same level of confidence in modeling these systems that they have, or should have, in modeling conventional systems, based on the results presented here.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗
DOE OSTI2022

Psilocybin: crystal structure solutions enable phase analysis of prior art and recently patented examples

Psilocybin {systematic name: 3-[2-(dimethylamino)ethyl]-1 H -indol-4-yl dihydrogen phosphate} is a zwitterionic tryptamine natural product found in numerous species of fungi known for their psychoactive properties. Following its structural elucidation and chemical synthesis in 1959, purified synthetic psilocybin has been evaluated in clinical trials and has shown promise in the treatment of various mental health disorders. In a recent process-scale crystallization investigation, three crystalline forms of psilocybin were repeatedly observed: Hydrate A, Polymorph A, and Polymorph B. The crystal structure for Hydrate A was solved previously by single-crystal X-ray diffraction. This article presents new crystal structure solutions for the two anhydrates, Polymorphs A and B, based on Rietveld refinement using laboratory and synchrotron X-ray diffraction data, and density functional theory (DFT) calculations. Utilizing the three solved structures, an investigation was conducted via Rietveld method (RM) based quantitative phase analysis (QPA) to estimate the contribution of the three different forms in powder X-ray diffraction (PXRD) patterns provided by different sources of bulk psilocybin produced between 1963 and 2021. Over the last 57 years, each of these samples quantitatively reflect one or more of the hydrate and anhydrate polymorphs. In addition to quantitatively evaluating the composition of each sample, this article evaluates correlations between the crystal forms present, corresponding process methods, sample age, and storage conditions. Furthermore, revision is recommended on characterizations in recently granted patents that include descriptions of crystalline psilocybin inappropriately reported as a single-phase `isostructural variant.' Rietveld refinement demonstrated that the claimed material was composed of approximately 81% Polymorph A and 19% Polymorph B, both of which have been identified in historical samples. In this article, we show conclusively that all published data can be explained in terms of three well-defined forms of psilocybin and that no additional forms are needed to explain the diffraction patterns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗