Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time discretization”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Compensatory reserve detects subclinical shock with more expeditious prediction for need of life‐saving interventions compared to systolic blood pressure and blood lactate

Abstract Introduction We conducted a prospective observational study on 205 trauma patients at a level I trauma facility to test the hypothesis that a compensatory reserve measurement (CRM) would identify higher risk for progression to shock and/or need a life‐saving interventions (LSIs) earlier than systolic blood pressure (SBP) and blood lactate (LAC). Methods A composite outcome metric included blood transfusion, procedural LSI, and mortality. Discrete measures assessed as abnormal (ab) were SBP <90 mmHg, CRM <60%, and LAC >2.0. A graded categorization of shock was defined as: no shock (normal [n] SBP [n‐SBP], n‐CRM, n‐LAC); sub‐clinical shock (ab‐CRM, n‐SBP, n‐LAC); occult shock (n‐SBP, ab‐CRM, ab‐LAC); or overt shock (ab‐SBP, ab‐CRM, ab‐LAC). Results Three patients displayed overt shock, 53 displayed sub‐clinical shock, and 149 displayed no shock. After incorporating lactate into the analysis, 86 patients demonstrated no shock, 25 were classified as sub‐clinical shock, 91 were classified as occult shock, and 3 were characterized as overt shock. Each shock subcategory revealed a graded increase requiring LSI and transfusion. Initial CRM was associated with progression to shock (odds ratio = 0.97; p < .001) at an earlier time than SBP or LAC. Conclusions Initial CRM uncovers a clinically relevant subset of patients who are not detected by SBP and LAC. Our results suggest CRM could be used to more expeditiously identify injured patients likely to deteriorate to shock, with requirements for blood transfusion or procedural LSI.

Convertino, Victor A.↗
DOE OSTI2024

Molecular catalyst coordinatively bonded to organic semiconductors for selective light-driven CO 2 reduction in water

The selective photoreduction of CO 2 in aqueous media based on earth-abundant elements only, is today a challenging topic. Here we present the anchoring of discrete molecular catalysts on organic polymeric semiconductors via covalent bonding, generating molecular hybrid materials with well-defined active sites for CO 2 photoreduction, exclusively to CO in purely aqueous media. The molecular catalysts are based on aryl substituted Co phthalocyanines that can be coordinated by dangling pyridyl attached to a polymeric covalent triazine framework that acts as a light absorber. This generates a molecular hybrid material that efficiently and selectively achieves the photoreduction of CO 2 to CO in KHCO 3 aqueous buffer, giving high yields in the range of 22 mmol g -1 (458 μmol g -1 h -1 ) and turnover numbers above 550 in 48 h, with no deactivation and no detectable H 2 . The electron transfer mechanism for the activation of the catalyst is proposed based on the combined results from time-resolved fluorescence spectroscopy, in situ spectroscopies and quantum chemical calculations.

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

Perturbation theory under the truncated Wigner approximation: How system-environment entanglement formation drives quantum decoherence

In discrete quantum systems, quantum decoherence is the disappearance of simple phase relations as a result of interactions with an environment. For many applications, the question is not necessarily how to avoid (inevitable) system-environment interactions, but rather how to design environments that optimally preserve a system's phase relations in spite of such interactions. The formation of system-environment entanglement is a major driving mechanism for decoherence, and a detailed understanding of this process could inform strategies for conserving coherence optimally. This requires scalable, flexible, and systematically improvable quantum dynamical methods that retain detailed information about the entanglement properties of the environment, yet very few current methods offer this combination of features. Here, in this work, we address this need by introducing a theoretical framework wherein we combine the truncated Wigner approximation with standard time-dependent perturbation theory, allowing for computing expectation values of operators in the combined system-environment Hilbert space. We demonstrate the utility of this framework by applying it to the spin-boson model, representative of qubits and simple donor-acceptor systems. For this model, our framework provides an analytical description of perturbative contributions to expectation values. We monitor how quantum decoherence at zero temperature is accompanied by entanglement formation with individual environmental degrees of freedom. Based on this entanglement behavior, we find that the selective suppression of low-frequency environmental modes is particularly effective for mitigating quantum decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Development of an MC&A toolbox for liquid-fueled molten salt reactors with online reprocessing (Final Report)

A critical barrier to the deployment of MSRs is the absence of a well-defined nuclear material control and accounting (MC&A) approach, a vital prerequisite to meet NRC licensing requirements as well as facilitating future international exports. Liquid-fueled MSR variants (especially those incorporating online refueling or reprocessing) present a unique set of challenges to traditional nuclear material control and accountancy. Unlike solid-fueled light-water reactors, traditional item-counting methods cannot be applied as an accountancy strategy. Rather, MC&A approaches to MSR variants (including both uranium and thorium fueled designs) are more analogous to bulk material handling facilities (e.g., enrichment and reprocessing); yet further complicating matters is the fact that fuel medium is also highly radioactive. Meanwhile, the space of MSRs covers a broad range of design parameters, including thermal and fast spectra designs; operation in actinide breeder or burner modes, choice of actinide fuel used (e.g., 235 U, 232 Th / 233 U, denatured 233 U, etc.), pool or loop-type configuration, and even different salt chemistry. Each of these design choices introduce significant challenges to MC&A approaches within MSR facilities. We propose to bridge this gap for liquid-fueled MSRs by developing a modular, component-based test framework for evaluating viable process monitoring and MC&A techniques specifically suited to liquid-fueled MSR system variants employing online refueling or reprocessing. This test platform will consist of a toolbox of independent process modules representing discrete physical units (such as the reactor core, off-gas processing, decay tanks, and actinide separation units), each with its own self-contained physics responsive to the input mass flow, along with appropriate measurement models that can be coupled to key flow points. These dynamic physical signatures thus afford the ability to test the viability and efficacy of potential accountancy techniques under the full range of reactor operating conditions. As process modules are connected via mass flows, the result is a reconfigurable, generic MSR mass flow model capable of serving as an MC&A test platform for a broad spectrum of possible MSR configurations. The resulting MSR MC&A toolbox will enable robust assessment of accountancy strategies for this unique facility type, including analysis of physical feedbacks arising both from depletion of the fuel over time as well as from potential off-normal events, including those introduced by equipment failures (e.g., a pump failure) as well as by malicious action (i.e., attempts to divert material). The proposed toolbox both addresses a critical needs area for the MPACT analysis toolkit while leveraging existing MPACT-sponsored tools, especially with respect to simulation of measurement and accountancy techniques for advanced fuel cycle facilities. Beyond enabling new analysis capabilities for MSR systems, the design of this toolbox will be to such to maximize compatibility with existing MPACT tools, such to enhance existing facility MC&A analysis capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Fast yaw optimization for wind plant wake steering using Boolean yaw angles

Abstract. In wind plants, turbines can be yawed into the wind to steer their wakes away from downstream turbines and achieve an overall increase in plant power. Mathematical optimization is typically used to determine the best yaw angles at which to operate the turbines in a plant. In this paper, we present a new heuristic to rapidly determine the yaw angles in a wind plant. In this method, we define the turbine yaw angles as Boolean – either yawed at a predefined angle or nonyawed – as opposed to the typical methods of defining yaw angles as continuous or with fine discretizations. We then optimize which turbines should be yawed with an algorithm that sweeps through the turbines from the most upstream to the most downstream. We demonstrate that our new Boolean optimization method can find turbine yaw angles that perform well compared to a traditionally used gradient-based optimizer for which the yaw angles are defined as continuous. There is less than 0.6 % difference in the optimized power between the two optimization methods for randomly placed turbine layouts and less than a 0.6 % difference in the optimal annual energy production between the two optimization methods for a real wind farm. Additionally, we show that our new method is much more computationally efficient than the traditional method. For plants with nonzero optimal yaw angles, our new method is generally able to solve for the turbine yaw angles 50–150 times faster, and in some extreme cases up to 500 times faster, than the traditional method.

17 WIND ENERGY↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2024

Chemically Driven Multistep Crystallization in the Synthesis of Sodium Yttrium Fluoride Via a Porous, Electrochemically Active Intermediate

Two-step crystallization mechanisms based on liquid–liquid phase separations followed by crystallization are commonly observed both in the laboratory and in nature. While this pathway quite often occurs as a result of a chemical reaction, the subsequent nucleation and growth are often considered as separate, discrete events from the reaction itself. We show this mechanism in the aqueous synthesis sodium yttrium fluoride, but by using a combination of experimental techniques and computational modeling, we show an additional step of solid-state chemical diffusion that is essential to the nucleation mechanism. In this system, we observe at least four distinct steps in the crystallization process, including (1) the segregation of aqueous ions into a dense liquid phase, (2) the formation of a metastable amorphous aggregate, (3) the continuous, gradual solid-state diffusion of sodium and fluoride ions into the amorphous aggregate toward a NaYF4 stoichiometry, and (4) the crystallization of a stable cubic sodium yttrium fluoride phase. Unlike previous descriptions of nucleation and growth, we find that the stoichiometry of the final solid phase evolves throughout the crystallization process rather than being determined at the time of the initial separation from solution. Further, this emphasizes that the chemical reaction cannot be assumed to be a separate event from the phase separation and growth, especially in compounds with variable stoichiometry. We also find that the amorphous aggregate that forms prior to the ion incorporation step adopts a porous, gel-like structure, which we isolated and showed to be electrochemically active, allowing for its potential use as a battery anode in lithium and sodium ion batteries, among other potential applications.

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

Vector boson dark matter from trinification

We show how trinification models based on the gauge group SU(3) C x SU(3) L x SU(3) R realized near the TeV scale can provide naturally a variety of dark matter (DM) candidates. These models contain a discrete T parity which may remain unbroken even after spontaneous symmetry breaking. The lightest T-odd particle, which could be a fermion, a scalar, or a gauge boson, can constitute the dark matter of the universe. This framework naturally admits a doublet-singlet fermionic DM, a singlet scalar DM, or a vector boson DM. Here we develop the vector boson DM scenario wherein the DM couples off-diagonally with the usual fermions and vector-like fermions present in the theory. We show consistency of this framework with dark matter relic abundance and direct detection limits as well as LHC constraints. We derive upper limits of 900 GeV on the vector gauge boson DM mass and 4.5 TeV on the vector-like quark masses. We also show the consistency of spontaneous gauge symmetry breaking down to Standard Model times an extra U(1) while preserving the T-parity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Historical review and proof-of-concept future method demonstration of adaptive mesh refinement in nuclear engineering for increased fidelity and computational efficiency

As the nuclear industry's use of computational tool increases, the need for increased fidelity and computational efficiency is well known. While most approaches to increased fidelity rely on applying a fine mesh over the problem domain, a more efficient method is to apply an adaptive mesh refinement (AMR) algorithm to the mesh definition. In the field of nuclear engineering, AMR has previously been used in conjunction with deterministic methods, including: S{sub N} transport methods, Lattice Boltzmann Methods, and COMSOL. The future of AMR in nuclear engineering is to couple it to a Monte Carlo code with the goal of reducing calculation time. A proof-of-concept example yielded positive results for using the gradient of the flux as a refinement criteria. The refinement criteria was varied from 0.01 to 0.10, which yielded a recommended range of 0.01 to 0.04, and the number of refinement iterations was varied from 0 to 7, with diminishing returns seen after 5 iterations. After the success of the proof-of-concept exercise, work began on creating a full program coupling MCNP6.2 and the AMR algorithm in the deal.II library. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Mixed-Integer Linear Programming Formulation with Embedded Machine Learning Surrogates for the Design of Chemical Process Families

In previous work, we introduced process family design. The main idea is to design a platform of common elements, and, allowing us to capture additional cost savings, simultaneously design a family of processes, and reducing both engineering and deployment timelines. We formulate this as an optimization problem, specifically a nonlinear generalized disjunctive program (GDP). We have proposed two approaches for reformulating and solving this problem: one based on full-discretization of the design space and one that uses Machine Learning (ML) surrogates to replace the nonlinear process models. Using ML surrogates to predict required system costs and performance indicators allows us to reformulate the nonlinearities in the GDP generate an efficient MILP formulation. In this work, we apply the ML surrogate approach to two case studies. One case study involves designing a family of carbon capture systems to cover a set of different flue gas flow rates and inlet CO 2 concentrations, where we consider the absorber and stripper as common unit module types. The second case study focuses on a water-desalination process, where we design a family of these processes for a variety of salt concentrations and flow rates. In both of these case studies, we demonstrate a scalable optimization approach that enables the design of multiple processes simultaneously, reducing the time-to-market and overall costs by maximizing the cost savings due to both economies of scale and economies of numbers.

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

Scalable Nano-Scaffold SOFC Anode Architecture Enabling Direct Hydrocarbon Utilization

This project is based on WVU’s pending patents, technology and aims to design and modify the internal surfaces of the Ni/YSZ anode from currently commercially viable Solid Oxide Fuel Cells (SOFCs) using the additive manufacturing process of Atomic Layer Deposition (ALD). The surface architecture/scaffold added onto the internal surface of the anode possesses an engineered nanostructure but it features only commonly-used oxide conductors and electro-catalyst materials. The surface layer possesses a minimum thickness of ~2-40 nm and is solely designed to control the surface reforming reactions and to increase catalytic activity. Three-dimensional (3D) nano scaffold architectures with the noble metal nano-catalyst, low-cost bimetallic catalytic alloys, and nano-scale ionic conducting oxide fully compatible with the state-of-the-art Ni/YSZ anode, were applied to the internal surface of the entire porous SOFC anode using ALD. In the present work, the surface scaffold architecture is essentially multi-functional at the nano-scale, facilitated by the multiple heterostructured interfaces. It will significantly enhance the power density and cell durability for direct hydrocarbon utilization by (1) increasing the number of electrochemical reaction sites to enhance the hydrogen/hydrocarbon oxidation reactions; (2) reducing carbon formation; (3) mitigating the coarsening of backbone Ni phase and the oxidation attack of Ni from oxidants (e.g., H2O, CO2); and (4) promoting the internal reforming capabilities, especially for natural gas applications. ALD is employed to generate stable anode surface architectures that are uniform, precisely controllable at the atomic scale, and accurately repeatable for processing. The engineered anode surface nano-scaffold architecture was cataloged and analyzed using High-resolution Transmission Electron Microscopy (TEM), and cell power/durability performance assessed via comprehensive electrochemical performance testing with commercial specimens and relevant environments using hydrocarbon fuels. To the best of our knowledge, this project is the First Report on ALD of Ni/YSZ. The actual achievement of this Project includes (1). Successful demonstration of 7 types of ALD layers on Ni/YSZ anode, including Co, Ni, Mn, Pt, Ru, ZrOx and multi-functional nano-composite. (1). Conformal coating and subsequently spontaneously pinning the discrete nano-catalyst, including the precious metal nano-catalyst and the Ni and Co catalysts, on the YSZ surface upon the electrochemical operation in the reduced atmosphere. Those nano-catalysts on the ionic-conducting YSZ provided excellent sites for promoting internal reforming; (2). Demonstrated ALD coating increased both catalytic activity and conductivity of Ni/YSZ. Conformal coating provided dopants and introduced additional electrical conducting pathways on the YSZ ionic conductor. The doped surface layer of YSZ with mixed conductivity thus further introduces the active triple phase boundaries adjacent to the ALD-coated nano-catalysts such as Pt, Co, and Ni that are pinned on the YSZ surface. The nano-composite ALD coating on Ni/YSZ anode has significantly increased cell durability; and (3). ALD coating of Ni/YSZ anode increased the power density of the entire cell by 300%. For a long time, the SOFC performance, such as the power density, was deemed hindered by the cathode. The sluggish oxygen reduction reaction (ORR) in the cathode was deemed as hindering the power density of the SOFCs. For the anode-supported commercial SOFCs, the cell performance is considered to be limited by the cathode's performance. For the first time in the field of SOFC, this project has demonstrated that (1). the performance of commercial SOFCs can be further increased by the ALD coating on Ni/YSZ anode backbone. (2). ALD coating on Ni/YSZ fuel electrodes results in the enhancement of power density, and increased reliability, robustness, and endurance of SOFCs, for their application using both hydrogen and hydrocarbon fuels over the entire operating temperature range of 650-800ºC for the inherently functional commercial cells. (3). ALD coating provides alternative approaches of exsolutions for introducing the stable catalyst onto the internal surface of the Ni/YSZ electrode. ALD coating could be much more versatile than exsolution in employing the catalysts with various chemistries onto the various backbones. (4). Due to the negligible amount of ALD materials coated onto the internal surface of the porous cathode of the as-fabricated cells, a peak power density increase up to 300 % induced by ALD coating was simultaneously achieved in terms of both power density and specific power. (5). The ALD coating developed through this project was applied to both the SOFC and Solid Oxide Electrolysis Cells (SOEC). SOEC’s face a similar but more demanding need to improve the fuel electrode's performance. It opens further research directions for electrocatalytic surface nanoionics with a wide range of chemistry. It will revolutionize our ability to render the formation of a nanostructured electrode that has been constantly pursued yet barely achieved for practical SOFC/SOEC applications. The research is also immediately transformative since both the preliminary data and the proposed work are on the direct implantation of nanoionics into the state-of-the-art inherently functional SOCs. It represents an immediate impact on the commercial sectors in SOC technology since the applied ALD processing is computer-controlled ALD coating using the commercial ALD systems, and it is scalable to both the single cells and SOC stacks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diversion Path Analysis: A Proposed Methodology to Develop an MC&A Approach for Liquid-Fueled Molten Salt Reactors

Nuclear material control and accounting (MC&A) is a critical element of both the US Nuclear Regulatory Commission (NRC) and US Department of Energy (DOE)’s domestic safeguards and security requirements. NRC licensees are required, under Title 10 of the Code of Federal Regulations (10 CFR) Part 74 to establish and maintain an MC&A program that captures and records the quantities and locations of special nuclear material (SNM) at the facility. Along with physical protection, MC&A is a key element of domestic nuclear material safeguards that enables the NRC to ensure that SNM is controlled and accounted for. SNM, per 10 CFR Part 74, refers to plutonium, 233 U, and uranium enriched in the isotope 233 U or 235 U, but does not include source material. Periodic physical inventories, coupled with material balance evaluations, are effective and demonstrated tools to account for and detect theft or diversion of SNM in facilities containing SNM in bulk material form (i.e., not in discrete, countable items). Historically in the United States, these types of facilities have included fuel fabrication, conversion, and enrichment facilities. In comparison, reactors have relied on item counting of assemblies and control of SNM while in containment (e.g., a sealed reactor pressure vessel) because, to date, reactor fuel has been in item form. In liquid-fueled molten salt reactors (MSRs), unlike traditional light water reactors (LWRs) or bulk facilities, bulk SNM quantities can change significantly during operation as a result of depletion and transmutation. This introduces challenges to the use of traditional periodic physical inventories and material balance evaluations to detect theft or diversion of SNM in reactors that use SNM in bulk material form. Liquid-fueled (i.e., salt-fueled) MSR facilities are MSRs that use SNM within a salt eutectic as the fuel. The SNM is in a bulk material form any time it is outside of fresh or spent fuel storage containers. Some examples of when SNM will be in bulk form in the facility are during addition of fuel to the reactor system, while fuel is circulating in operation, and while fuel is in a drain tank. Periodic physical inventories and material balance evaluations can likely be effectively applied to many portions of an MSR facility, including all areas where depletion and transmutation are not significantly changing the quantities of SNM within the control area. Within an MSR facility, this would include fresh fuel receipt and loading, waste streams that may contain SNM, irradiated fuel storage outside of the reactor core, and any irradiated fuel processing that may happen after SNM has been removed from the reactor. All of these process steps could rely on measurements of SNM quantities compared with documented inventories. Any discrepancies from predicted (i.e., book) inventories and measured inventories could be quantified as inventory differences, consistent with traditional MC&A guidance from the NRC (e.g., in NUREG-1065 Revision 2, NUREG-2159 Revision 1, and RG 5.29 Revision 2). Within the reactor system, additions and removals to the book inventory include depletion of the SNM (e.g., fission of 235 U), which complicates the use of physical inventories. SNM control, however, can also likely be effectively applied to detect theft of SNM throughout a liquid-fueled MSR facility. To complement these approaches, prior technical reports have identified that a diversion path analysis may be a useful, risk-informed, and performance-based tool to determine suitable elements of an MC&A approach for the reactor system within a liquid-fueled MSR facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Redefining Resource Adequacy for Modern Power Systems: A Report of the Redefining Resource Adequacy Task Force

Today's rapidly increasing levels of wind, solar, storage, and load flexibility require the industry to rethink reliability planning and resource adequacy methods for modern power systems. Periods with a risk of shortfall often no longer coincide with peak demand - reliability risks are less about peak load and more about the daily setting of the sun, extended cloud cover, wind speeds, cold snaps, and heat waves. In addition, demand is increasingly flexible. Key resources are time-sensitive, as batteries need time to recharge and electricity customers can only be asked to provide demand response for just so long. And reliability failures are often correlated - with one another and with the weather. Two driving factors require the industry to reconsider its analytical approach for resource adequacy: (1) Chronological grid operations: The increasing importance of variable renewable resources (such as wind and solar) and of energy-limited resources (storage and demand response) make it essential to understand the full year of chronological operation of the grid. Specific attention must be paid to hourly, seasonal, and inter-annual resource variability. The sequence of the variability is key, as energy-limited resources such as batteries or demand response require either a preceding period or subsequent period of high production to be useful for grid reliability. (2) Correlated events: Historically, resource adequacy analysis focused on shortfalls caused by random, discrete mechanical failures of large generating units. In contrast, shortfalls today are often caused by multiple, correlated events caused by common weather patterns. Resource adequacy analysis must increasingly shift its focus to these correlated events. The redesign of resource adequacy methods will benefit from a set of guiding principles to better allow for sharing of insights and best practices, interregional resource coordination, and a smoother regulatory process for resource procurement. The objective of this report is to move this redesign forward. It provides an overview of key drivers changing the way resource adequacy needs to be evaluated, identifies shortcomings of conventional approaches, and outlines first principles for practitioners to consider as they adapt their approaches. The central message is: what got us here won't get us there.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Differential geometric approaches to momentum-based formulations for fluids [Slides]

This SAND report documents CIS Late Start LDRD Project 22-0311, "Differential geometric approaches to momentum-based formulations for fluids". The project primarily developed geometric mechanics formulations for momentum-based descriptions of nonrelativistic fluids, utilizing a differential geometry/exterior calculus treatment of momentum and a space+time splitting. Specifically, the full suite of geometric mechanics formulations (variational/Lagrangian, Lie-Poisson Hamiltonian and Curl-Form Hamiltonian) were developed in terms of exterior calculus using vector-bundle valued differential forms. This was done for a fairly general version of semi-direct product theory sufficient to cover a wide range of both neutral and charged fluid models, including compressible Euler, magnetohydrodynamics and Euler-Maxwell. As a secondary goal, this project also explored the connection between geometric mechanics formulations and the more traditional Godunov form (a hyperbolic system of conservation laws). Unfortunately, this stage did not produce anything particularly interesting, due to unforeseen technical difficulties. There are two publications related to this work currently in preparation, and this work will be presented at SIAM CSE 23, at which the PI is organizing a mini-symposium on geometric mechanics formulations and structure-preserving discretizations for fluids. The logical next step is to utilize the exterior calculus based understanding of momentum coupled with geometric mechanics formulations to develop (novel) structure-preserving discretizations of momentum. This is the main subject of a successful FY23 CIS LDRD "Structure-preserving discretizations for momentum-based formulations of fluids".

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗
DOE OSTI2021

A Modified Vegetation Photosynthesis and Respiration Model (VPRM) for the Eastern USA and Canada, Evaluated With Comparison to Atmospheric Observations and Other Biospheric Models

Atmospheric CO 2 measurements from a dense surface network can help to evaluate terrestrial biosphere model (TBM) simulations of Net Ecosystem Exchange (NEE) with two key benefits. First, gridded CO 2 flux estimates can be evaluated over regional scales, not possible using flux tower observations at discrete locations for model evaluation. Second, TBM ability to explain atmospheric CO 2 fluctuations due to the biosphere can be directly tested, an important objective for anthropogenic emissions monitoring using atmospheric observations. Here, we customize the Vegetation Photosynthesis and Respiration Model (VPRM) for an eastern North American domain with strong biological activity upwind of urban areas. Parameters are optimized using flux tower observations from a historical database with sites in (and near) the domain. In addition, the respiration model (originally a linear function of temperature) is modified to account for impacts of changing foliage, non-linear temperature, and water stress. Flux estimates from VPRM, the Carnegie-Ames-Stanford Approach (CASA) model and the Simple Biosphere Model v4 (SiB4), are convolved with footprints from atmospheric transport models for evaluation with CO 2 observations at 21 towers in the domain, with roughly half of the towers used here for the first time. Results show that the new respiration model in VPRM helps to correct a growing season sink bias in the atmosphere associated with underestimated summertime respiration using the original model with annual parameters. The new VPRM also better explains fine-scale atmospheric CO 2 variability compared to other TBMs, due to higher resolution diagnostic phenology, the new respiration model, domain-specific parameters, and high-quality input data sets.

54 ENVIRONMENTAL SCIENCES↗