Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V CODES”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

QR-code optical encryption in the scheme with spatially incoherent illumination based on two micromirror light modulators

An optical encryption scheme with spatially incoherent illumination based on two micromirror light modulators has been experimentally implemented for the first time. Currently, such modulators are the fastest tools for spatio-temporal light modulation; their high frame rates provide possibilities for developing optical encryption systems with a bandwidth of several gigabites per second. (optical encryption)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Benchmark Specification for FFTF LOFWOS Test #13

The Fast Flux Test Facility (FFTF) at the Hanford site in Washington was designed by the Westinghouse Electric Corporation for the U.S. Department of Energy. FFTF was a 400 MW thermal, oxide-fueled, liquid sodium cooled test reactor, built to assist development and testing of advanced fuels and materials for fast breeder reactors. After reaching criticality in 1980, FFTF operated until 1992, providing the U.S. Department of Energy (DOE) with the means to test fuels, materials, and other components in a fast neutron flux environment. In July 1986, a series of unprotected transients (with the plant protection system intentionally disabled) were performed in FFTF as part of the passive safety demonstration program. Among these were thirteen loss of flow without scram (LOFWOS) tests. The goals of this program included confirming the liquid metal reactor safety margins, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of specific design features. The test defined in this benchmark is LOFWOS Test #13, which was initiated at 50% power and 100% flow with the pump pony motors turned off. This benchmark specification is intended to support collaborative efforts within international partnerships on the validation of simulation tools and models in the area of Sodium-cooled Fast Reactor (SFR) safety. Validated tools and models are needed to evaluate SFR inherent safety characteristics and assess the impact of passive design features in response to accident initiators. Comparisons with experimental data and the results of safety analyses from other groups create unique opportunities to improve predictive capabilities of computational codes and methods for SFR modeling and simulation. The conditions of the LOFWOS test along with the feedback from FFTF’s limited free bow core restraint system and the novel passive reactivity control Gas-Expansion Modules (GEMs) pose a very challenging and uniquely valuable benchmark exercise.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Physics↗
DOE OSTI2022

Spinor $GW$ Bethe-Salpeter calculations in BerkeleyGW: Implementation, symmetries, benchmarking, and performance

Computing the GW quasiparticle band structure and Bethe-Salpeter equation (BSE) absorption spectra for materials with spin-orbit coupling have commonly been done by treating GW corrections and spin-orbit coupling (SOC) as separate perturbations to density-functional theory. However, accurate treatment of materials with strong spin-orbit coupling (such as many topological materials of recent interest, and thermoelectrics) often requires a nonperturbative approach using spinor wave functions in the Kohn-Sham equation and GW/BSE. Such calculations have only recently become available, in particular for the BSE. Here, we have implemented this approach in the plane-wave pseudopotential GW/BSE code BerkeleyGW, which is highly parallelized and widely used in the electronic-structure community. We present reference results for quasiparticle band structures and optical absorption spectra of solids with different strengths of spin-orbit coupling, including Si, Ge, GaAs, GaSb, CdSe, Au, and Bi 2 Se 3 . The calculated quasiparticle band gaps of these systems are found to agree with experiment to within a few tens of meV. SOC splittings are found to be generally in better agreement with experiment, including quasiparticle corrections to band energies. The absorption spectrum of GaAs is not significantly impacted by the inclusion of spin-orbit coupling due to its relatively small value (0.2 eV) in the Λ direction, while the absorption spectrum of GaSb calculated with the spinor GW/BSE captures the large spin-orbit splitting of peaks in the spectrum. For the prototypical topological insulator Bi 2 Se 3 , we find a drastic change in the low-energy band structure compared to that of DFT, with the spinorial treatment of the GW approximation correctly capturing the parabolic nature of the valence and conduction bands after including off-diagonal self-energy matrix elements. We present the detailed methodology, approach to spatial symmetries for spinors, comparison against other codes, and performance compared to spinless GW/BSE calculations and perturbative approaches to SOC. This work aims to spur further development of spinor GW/BSE methodology in excited-state research software and enables a more accurate and detailed exploration of electronic and optical properties of materials containing elements with large atomic numbers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Feasibility and performance of the staged Z-pinch: A one-dimensional study with FLASH and MACH2

Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (>300) but achieve fuel ion temperatures of many kilo-electron volts. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

JETSCAPE framework: p+p results

The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in p + p collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and dijet mass cross sections. These observables in p+p collisions provide the baseline for their counterparts in nuclear collisions. Quantifying the level of agreement of JETSCAPE results with p+p data is thus necessary for meaningful applications of JETSCAPE to A+A collisions. The calculations use the JETSCAPE PP19 tune, defined in this paper, based on version 1.0 of the JETSCAPE framework. For the observables discussed in this work calculations using JETSCAPE PP19 agree with data over a wide range of collision energies at a level comparable to standard Monte Carlo codes. These results demonstrate the physics capabilities of the JETSCAPE framework and provide benchmarks for JETSCAPE users.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

RSE-Ops

If DevOps is the intersection of "Developer" and "Operations," then how does this concept map to high performance computing? The answer is RSE-ops, or the intersection of Research Software Engineering and Operations. Research Software Engineers (RSEs) are the individuals writing code for scientific software, and supporting researchers to use codes on high performance computing systems. RSE-ops, then, generally refers to best practices for development and operations of scientific software, which typically happens in a high performance computing environment. A comparison can be made to Cloud Native Devops, which is a similar term that has emerged around Cloud Native Computing. It refers to the application of DevOps practices to the cloud.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Time-Energy Uncertainty Relation for Noisy Quantum Metrology

Detection of very weak forces and precise measurement of time are two of the many applications of quantum metrology to science and technology. To sense an unknown physical parameter, one prepares an initial state of a probe system, allows the probe to evolve as governed by a Hamiltonian 𝐻 for some time 𝑡, and then measures the probe. If 𝐻 is known, we can estimate 𝑡 by this method; if 𝑡 is known, we can estimate classical parameters on which 𝐻 depends. The accuracy of a quantum sensor can be limited by either intrinsic quantum noise or by noise arising from the interactions of the probe with its environment. In this work, we introduce and study a fundamental trade-off, which relates the amount by which noise reduces the accuracy of a quantum clock to the amount of information about the energy of the clock that leaks to the environment. Specifically, we consider an idealized scenario in which a party Alice prepares an initial pure state of the clock, allows the clock to evolve for a time that is not precisely known, and then transmits the clock through a noisy channel to a party Bob. Meanwhile, the environment (Eve) receives any information about the clock that is lost during transmission. We prove that Bob’s loss of quantum Fisher information about the elapsed time is equal to Eve’s gain of quantum Fisher information about a complementary energy parameter. We also prove a similar, but more general, trade-off that applies when Bob and Eve wish to estimate the values of parameters associated with two noncommuting observables. We derive the necessary and sufficient conditions for the accuracy of the clock to be unaffected by the noise, which form a subset of the Knill-Laflamme error-correction conditions. A state and its local time-evolution direction, if they satisfy these conditions, are said to form a metrological code. We provide a scheme to construct metrological codes in the stabilizer formalism. We show that there are metrological codes that cannot be written as a quantum error-correcting code with similar distance in which the Hamiltonian acts as a logical operator, potentially offering new schemes for constructing states that do not lose any sensitivity upon application of a noisy channel. We discuss applications of the trade-off relation to sensing using a quantum many-body probe subject to erasure or amplitude-damping noise.

metrology↗
DOE OSTI2020

FAST-1.0 Software Release Document

The purpose of this document is to detail the release of FAST-1.0 and the documents that accompany the development and testing of FAST-1.0. Design requirements were for FAST-1.0 developed in the NRC Statement of Work (SOW) for NRC Agreement Number NRC-HQ-25-14-D-001 and are transcribed to this document as part of the Software Quality Assurance Plan for the FAST Code System, PNNL-28767. FAST 1.0 represents the merger of the FRAPCON and FRAPTRAN codes and as such FAST-1.0 performs steady state and transient fuel performance calculations described below. FAST-1.0 calculates the steady state response of light-water reactor fuel rods during long-term burnup. The code calculates temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include: 1) heat conduction through fuel and cladding to the coolant; 2) cladding elastic and plastic deformation; 3) fuel-cladding mechanical interaction; 4) fission gas release from the fuel and rod internal pressure; and 5) cladding oxidation. FAST is used to perform independent calculations for regulatory evaluations of fuel performance under normal operation and anticipated operational occurrences (AOOs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Deuteron-induced reactions on natural Zr from threshold to 50 MeV: production of 86g Y

Two stacks of thin Zr foils were irradiated with 30 and 50 MeV deuterons, respectively, using the Lawrence Berkeley National Laboratory 88-Inch Cyclotron, and 19 excitation functions for nat Zr(d,x) reactions were measured over a beam energy range of 6.3–47.64 MeV, where the independent cross sections for nat Zr(d,x) 88 Nb and nat Zr(d,x) 86m,g Y were measured for the first time. The well-characterized nat Fe(d,x) 56 Co, nat Ni(d,x) 56 Co, nat Ni(d,x) 58 Co, nat Ni(d,x) 61 Cu, nat Ti(d,x) 46 Sc and nat Ti(d,x) 48 V monitor reactions were used to determine the deuteron beam current throughout the stacks. All cross sections were determined using High Purity Germanium (HPGe) detector γ-ray spectroscopy. A variance minimization technique was employed to simultaneously constrain the deuteron beam currents with multiple monitor reactions, thus reducing systematic uncertainties. An additional 16 channels are reported for reactions on the nickel, titanium, and iron monitor foils, leading to a total of 35 excitation functions, with seven reaction channels reported for the first time in this work. The measured excitation functions are compared to calculations provided by the reaction modeling codes TALYS – 2.0, ALICE – 2020, CoH – 3.5.3 and EMPIRE – 3.2.3, as well as the TENDL – 2023 data library. The degree of agreement between theory and experiments is discussed. The possible production of the important PET radionuclide 86g Y via the nat Zr(d,x) route was critically examined. The physical yields for nat Zr(d,x) 86 Y and other yttrium isotopes produced were calculated and compared to other production pathways. Due to high-level of associated radionuclide impurities, this route cannot deliver 86g Y suitable for medical applications.

86gY↗
DOE OSTI2023

Efforts to enhance reproducibility in a human performance research project

Background: Ensuring the validity of results from funded programs is a critical concern for agencies that sponsor biological research. In recent years, the open science movement has sought to promote reproducibility by encouraging sharing not only of finished manuscripts but also of data and code supporting their findings. While these innovations have lent support to third-party efforts to replicate calculations underlying key results in the scientific literature, fields of inquiry where privacy considerations or other sensitivities preclude the broad distribution of raw data or analysis may require a more targeted approach to promote the quality of research output. Methods: We describe efforts oriented toward this goal that were implemented in one human performance research program, Measuring Biological Aptitude, organized by the Defense Advanced Research Project Agency's Biological Technologies Office. Our team implemented a four-pronged independent verification and validation (IV&V) strategy including 1) a centralized data storage and exchange platform, 2) quality assurance and quality control (QA/QC) of data collection, 3) test and evaluation of performer models, and 4) an archival software and data repository. Results: Our IV&V plan was carried out with assistance from both the funding agency and participating teams of researchers. QA/QC of data acquisition aided in process improvement and the flagging of experimental errors. Holdout validation set tests provided an independent gauge of model performance. Conclusions: In circumstances that do not support a fully open approach to scientific criticism, standing up independent teams to cross-check and validate the results generated by primary investigators can be an important tool to promote reproducibility of results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Preliminary design of control rods in the single-fluid double-zone thorium molten salt reactor (SD-TMSR)

Recent studies on Molten Salt Reactors (MSRs) showed that the excess reactivity at the beginning of the operation is large for many fueling strategies and must be compensated by a reactivity control system. The current work introduces a reliable safety system based on control rods in addition to the online feed system reactivity control in the Single-fluid Double-zone Thorium-based Molten Salt Reactor (SD-TMSR). Three different initial fissile loadings are considered: 233 U, reactor-grade Pu, and transuranic (TRU) elements as a startup fuel. We applied six different absorbing materials to investigate the main operational and safety parameters in the SD-TMSR: natural B 4 C, enriched B 4 C with 90% 10 B, HfB 2 , HfH 1.62 , Eu 2 O 3 , and Gd 2 O 3 . The present work focuses on control rod design, integral and differential control rod worth, shutdown margin, and shadowing effects at steady-state. We employed the SERPENT-2 Monte-Carlo code to calculate the reactivity worth and analyze the performance of the reactivity control system. We showed that 233 U and reactor-grade Pu startup cores maintain adequate shutdown margin with all considered absorbers. Finally, this report proposes a design of control rod clusters that compensate the excess reactivity of the SD-TMSR loaded with different initial fissile material.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗