SEARCH · Engineering Papers

Results for “conditionally specified models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Fracture Toughness and Fatigue Crack Growth Rate Testing of Baffle-Former Bolts Harvested from a Westinghouse Two-Loop Downflow Type PWR

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Microstructural Characterization of the Second High Fluence Baffle-Former Bolt Retrieved from a Westinghouse Two-loop Downflow Type PWR

As one of the pressurized water reactor (PWR) internal components, baffle-former bolts (BFBs) are subjected to significant mechanical stress and neutron irradiation from the reactor core during the plant operation. Over the long operation period, these conditions lead to potential degradation and reduced load-carrying capacity of the bolts. In support of evaluating long-term operational performance of materials used in core internal components, the Oak Ridge National Laboratory (ORNL), through the Department of Energy (DOE), Light Water Reactor Sustainability (LWRS) Program, Materials Research Pathway (MRP) has harvested two high fluence BFBs from a commercial Westinghouse two-loop downflow type PWR. The two bolts of interest, i.e. bolts # 4412 and 4416, were withdrawn from service in 2011 as part of a preventative replacement plan. No identification of cracking or potential damage was found for these bolts during their removal in 2011. However, the bolts required a lower torque for removal from the baffle structure than the original torque specified during installation. Irradiation displacement damage levels in the bolts range from 15 to 41 displacements per atom. The goal of this project is to perform detailed microstructural and mechanical property characterization of BFBs following in-service exposures. The information from these bolts will be integral to the LWRS program initiatives in evaluating end of life microstructure and properties. Furthermore, valuable data will be obtained that can be incorporated into model predictions of long-term irradiation behavior and compared to results obtained in high flux experimental reactor conditions. In this report, we present our latest study in FY22 on microstructural characterizations of the second high fluence baffle-former bolt, i.e., bolt # 4412. Analytical electron microscopy and atom probe tomography characterization were performed. The radiation-induced defects in the material add to the large wealth of knowledge for neutron-induced defects in 304/316 grades of stainless steels, specifically for radiation-induced precipitation after high fluence commercial PWR irradiation. The main findings are summarized as follows: 1) The cavity size was considerably larger in the bolt thread section than in the bolt head, with the bolt thread section having a bimodal distribution of cavities greater than ~6 nm in diameter and less than ~3 nm in diameter. The bolt head only had the small-sized cavities. In addition, there was a denuded zone of large cavities near grain boundaries in the thread section of the bolt. 2) Radiation-induced precipitation in the BFB #4412 was highly complex, with the volume fraction, size, and number density of Ni/Si and Cu-rich precipitates depending strongly on the radiation temperature/dose. In many cases, co-precipitates of adjoined clusters were found with Ni/Si-rich precipitates sandwiched between Cu-rich clusters and Mo/Cr/P-rich clusters. 3) Solute segregation out of solution was highest for most solutes in the thread section of the bolt #4412 with the exception of Cu, which experienced more separation out of solution into Curich clusters in the bolt head section. This highlights the difference in the mechanisms for precipitation of Ni/Si clusters, which have the Ni 3 Si phase composition, and precipitation of Cu-rich clusters. 4) There appear to be multiple simultaneous influences that affect the microstructural variation along the length of the bolt that overcomes the ~2X difference in irradiation dose between the bolt head and the bolt thread. The irradiation temperature, thermal/fast neutron ratio variation, potential strain gradient, and exposure to PWR coolant water that each section of the bolt sees may have more influence on the microstructural evolution than the total irradiation dose. The bolt thread and shank, with higher temperature, higher relative fast neutron flux, higher strain, and exposure to coolant but lower dose, underwent more enhanced cavity formation, precipitation, and solute segregation than the bolt head section.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

(SPT-)LSM theorems from projective non-invertible symmetries

Projective symmetries are ubiquitous in quantum lattice models and can be leveraged to constrain their phase diagram and entanglement structure. In this paper, we investigate the consequences of projective algebras formed by non-invertible symmetries and lattice translations in a generalized 1+1 1 + 1 D quantum XY model based on group-valued qudits. This model is specified by a finite group G G and enjoys a projective \mathsf{Rep}(G)× Z(G) 𝖱 𝖾 𝗉 ( G ) × Z ( G ) and translation symmetry, where symmetry operators obey a projective algebra in the presence of symmetry defects. For invertible symmetries, such projective algebras imply Lieb-Schultz-Mattis (LSM) anomalies. However, this is not generally true for non-invertible symmetries, and we derive a condition on G G for the existence of an LSM anomaly. When this condition is not met, we prove an SPT-LSM theorem: any unique and gapped ground state is necessarily a non-invertible weak symmetry protected topological (SPT) state with non-trivial entanglement, for which we construct an example fixed-point Hamiltonian. The projectivity also affects the dual symmetries after gauging \mathsf{Rep}(G)× Z(G) 𝖱 𝖾 𝗉 ( G ) × Z ( G ) sub-symmetries, giving rise to non-Abelian and non-invertible dipole symmetries, as well as non-invertible translations. We complement our analysis with the SymTFT, where the projectivity causes it to be a topological order non-trivially enriched by translations. Throughout the paper, we develop techniques for gauging \mathsf{Rep}(G) 𝖱 𝖾 𝗉 ( G ) symmetry and inserting its symmetry defects on the lattice, which are applicable to other non-invertible symmetries.

Pace, Salvatore D. (ORCID:0000000306093335)↗

DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Cavity evolution and void swelling in dual ion irradiated tempered martensitic steels

Microstructural evolutions in two variants of the F82H 8 wt%Cr tempered martensitic steel under dual (He + plus Fe 3+ ) ion irradiations at 500 °C have been characterized by transmission electron microscopy. Emphasis is on understanding the combined effects of displacements per atom (dpa) and the helium/dpa ratio (He/dpa) on cavity evolution and void swelling (f v ). The dual ion database includes 231 alloy-dpa-He/dpa conditions for irradiations up to ~ 82 dpa and 3700 appm helium. Here, we hypothesize that the large “scatter” observed in the f v data is primarily due to local microstructural variations, especially those that shield regions at shallower depths from the effects of injected interstitial atoms, which would otherwise reduce f v . More generally, our hypothesis is that the highest f v data, at a specified dpa and He/dpa, provides the most appropriate estimate of swelling for fusion neutron irradiation conditions. The dose dependence of f v can be described by an incubation dpa i , marking the onset of significant void growth, and a post-incubation swelling rate, f v ’ = Δf v /Δdpa. Both dpa i and f v ’ decrease with increasing He/dpa. The f v ’ data are quantitatively consistent with a simple model for defect partitioning between evolving bubble, void, and dislocation sinks strengths and biases. Most notably, the model correctly predicts both the incubation dpa i , and the subsequent f v ’, in spite of the fact that they reflect somewhat different physics. Our analysis also suggests that defect recombination is modest at 500 °C. Along with other considerations, this supports the hypothesis that, even at high dpa rates, dual ion irradiations can be used to reasonably emulate void nucleation and growth for a range of fusion-relevant neutron irradiation service conditions. A companion paper incorporates the results of this study in a reduced-order, data-driven void swelling model.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Machine learning assisted phase and size-controlled synthesis of iron oxide particles

Synthesis of iron oxides with specific phases and particle sizes is a crucial challenge in various fields, including materials science, energy storage, biomedical applications, environmental science, and earth science. However, despite significant advances in this area, much of the current palette of particle outcomes has been based on time-consuming trial-and-error exploration of synthesis conditions. The present study was designed to explore a very different approach to 1) predict the outcome of synthesis from specified reaction parameters based on using machine learning (ML) techniques, and 2) correlate sets of parameters to obtain products with desired outcomes by a newly designed recommendation algorithm. To achieve this, four ML algorithms were tested, namely random forest, logistic regression, support vector machine, and k-nearest neighbor. Among the models, random forest outperformed the others, attaining 96% and 81% accuracy when predicting the phase and size of iron oxide particles in the test dataset. Surprisingly, the permutation feature importance analysis revealed that volume, which may strongly relate to pressure, was one of the important features, along with precursor concentration, pH, temperature, and time, influencing the phase and size of iron oxide particles during synthesis. To verify the robustness of the random forest models, prediction and experimental results were compared based on 24 randomly generated methods in additive and non-additive systems not included in the datasets. The predictions of product phase and particle size from the models agreed well with the experimental results. Furthermore, a searching and ranking algorithm was developed to recommend potential synthesis parameters for obtaining iron oxide products with the desired phase and particle size from previous studies in the dataset. Furthermore, this study lays the foundation for a closed-loop approach in materials synthesis and preparation, beginning with suggesting potential reaction parameters from the dataset and predicting potential outcomes, followed by conducting experiments and analyses, and ultimately enriching the dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Differentially Private Generation of Social Networks via Exponential Random Graph Models

Many social networks contain sensitive relational information. One approach to protect the sensitive relational information while offering flexibility for social network research and analysis is to release synthetic social networks at a pre-specified privacy risk level, given the original observed network. In this work, we propose the DP-ERGM procedure that synthesizes networks that satisfy the differential privacy (DP) via the exponential random graph model (EGRM). We apply DP-ERGM to a college student friendship network and compare its original network information preservation in the generated private networks with two other approaches: differentially private DyadWise Randomized Response (DWRR) and Sanitization of the Conditional probability of Edge given Attribute classes (SCEA). The results suggest that DP-EGRM preserves the original information significantly better than DWRR and SCEA in both network statistics and inferences from ERGMs and latent space models. In addition, DP-ERGM satisfies the node DP, a stronger notion of privacy than the edge DP that DWRR and SCEA satisfy.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗

DOE OSTI2023

Agent-based simulation and child protection systems: Rationale, implementation, and verification

Simulation models are an important tool used in health care and other disciplines to support operational research and decision-making. In the child protection literature, simulation models are an under-utilized source of research evidence. Here, in this paper, we describe the rationale for and the development of an agent-based simulation of a child protection system in the US. Using the investigation, prevention service, and placement histories of 600,000 children served in an urban child welfare system, we walk the reader through the development of a prototype known as OSPEDALE. The governing equations built into OSPEDALE probabilistically simulate the onset of investigations. Then, drawing from empirical survival distributions, the governing equations trace the probability of subsequent interactions with the system (recurrence of maltreatment, service referrals, and placement) conditional on the characteristics of children, their assessed risk level, and prior child protection system involvement. As an initial test of OSPEDALE's utility, we compare empirical admission counts with counts generated from OSPEDALE. Though the verification step is admittedly simple, the comparison shows that OSPEDALE replicates the empirical count of new admissions closely enough to justify further investment in OSPEDALE. Management of public child protection systems is increasingly research evidence-dependent. The emphasis on research evidence as a decision-support tool has elevated evidence acquired through randomized clinical trials. Though important, the evidence from clinical trials represents only one type of research evidence. Properly specified, simulation models are another source of evidence with real-world relevance.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Modeling DOE Standard Canister Configurations with Updated Surface Chemistry

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗

DOE OSTI2025

Mineral Scale Formation during Crossflow Reverse Osmosis at Constant Flux and Constant Transmembrane Pressure Conditions

Mineral scale formation on membrane surfaces is a significant challenge in reverse osmosis water purification. Laboratory fouling experiments are typically run such that the transmembrane pressure (TMP) is fixed, and the permeate flux decreases over time as scales accumulate on the membrane surface. However, this change in flux means that the hydrodynamic conditions at the membrane surface are continuously changing, which could affect crystallization and foulant deposition processes. Operating under constant permeate flux conditions, in contrast, is advantageous because it keeps the hydrodynamic conditions relatively consistent, making it possible to compare how membrane properties (e.g., surface chemistry) affect fouling propensity. Industrial reverse osmosis operations are not run strictly in either constant TMP or constant flux mode; while they may start at a constant TMP, feed pressure may be periodically adjusted to maintain permeate water production within a specified range. Here, the scarcity of constant permeate flux reverse osmosis scaling experiments reported in the literature frustrates efforts to compare membrane fouling processes under constant TMP and constant flux conditions. For the first time, the evolution of the fouling layer resistance was compared as a function of cumulative permeate volume per membrane area during constant TMP and constant flux reverse osmosis filtrations. Scaling experiments were conducted by challenging commercial reverse osmosis membranes with a model feed solution nearly saturated with calcium sulfate dihydrate (gypsum). At low fluxes, the increase in fouling layer resistance was quantitatively similar for the two operational modes. In contrast, at high fluxes, the fouling layer resistance increased more rapidly in constant flux filtration than in constant TMP filtration. The mechanism of scale formation in constant TMP and constant flux operation was self-limiting and self-reinforcing, respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Estimating internal moisture generation rates in recently constructed, occupied homes in the southeastern United States

Internal moisture generation (IMG), or moisture generated by building occupants via activities such as respiration, cooking, bathing, and cleaning, is a critical input required for design, analysis, and simulation of building enclosure and heating, ventilation, and air conditioning (HVAC) systems. Based on previously published values, ASHRAE Standard 160-2021 provides guidance for estimating IMG rates for moisture control design analysis, which is based on occupied home datasets collected in the 1980s and 1990s. Residential energy use simulation software also utilizes estimates for IMG for energy rating, energy analysis, and code compliance calculations, as specified in ANSI/RESNET/ICC Standard 301. Data quantifying IMG rates in newer homes is useful in determining the continued relevance of current design and simulation guidance, and whether the guidance represents conditions found in new housing stock. ASHRAE Research Project 1844, conducted by the Florida Solar Energy Center (FSEC), a research institute of the University of Central Florida, estimated IMG rates in newer occupied homes built since 2013 in the southeastern US using a moisture balance model approach. Furthermore, the project obtained occupied home data in cooperation with a US Department of Energy Building America research project to characterize indoor air quality (IAQ) in newer US homes, along with presence, functionality, and occupant use of control measures. Full-scale laboratory homes operating with known IMG rates were used to validate a moisture balance model and quantify the accuracy of estimates obtained when the model is applied to data from occupied homes with unknown IMG rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Building a testable shear viscosity across the QCD phase diagram

Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio ηT/w of the Quark Gluon Plasma remains unknown. Here we use the Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate ηT/w at low temperatures and at finite baryon densities ρ B . We then match ηT/w to a QCD-based shear viscosity calculation within the deconfined phase to create a table across {T,μ B } for different cross-over and critical point scenarios at a specified location. We find that these new ηT/w(T,μ B ) values would require initial conditions at significantly larger ρB, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗