Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW 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.

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

At least 541 records · Page 30

Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1 Chg.1, Radioactive Waste Management, a Composite Analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This document details the application of the P2R Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the CA. The simulated flow field will support the simulation of fate and transport of contaminants for use in predicting dose as part of the CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predictive Flow Simulation with the P2R Model for the 200-IA-1 Preliminary Remediation Goal Saturated Zone Abstraction

In order to develop and evaluate alternatives for remediation of waste sites constituting the 200-IA-1 Operable Unit (OU), a focused feasibility study (FFS) is currently underway. The 200-IA-1 OU FFS requires preliminary remediation goals (PRGs) protective of groundwater be evaluated which includes estimates of fate and transport of contaminants in the groundwater from multiple sources within the Central Plateau Inner Area. This document details the application of the Plateau-to-River (P2R) Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the 200-IA-1 OU PRG calculation. The simulated flow field will support the simulation of fate and transport of contaminants for use in estimating peak unit-length concentration as part of the 200-IA-1 OU PRG calculation presented in a separate environmental calculation file (ECF) (ECF-200IA1-20-0108, Determination of Preliminary Remediation Goals for 200-IA-1 Source Operable Unit).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-temperature seals for supercritical carbon-dioxide (sCO 2 ) turbines (Final Report)

This is the final report for project DE-FE0031924 titled “High-temperature seals for supercritical carbon-dioxide (sCO 2 ) turbines.” The report provides a summary of the entire project efforts from October 2020 through December 2024 including the high-temperature commercial dry gas seal (DGS) tests and thermal modeling of Task 2, as well as the high-temperature, large-diameter seal design and high-temperature tests of large-diameter seals in Task 3. A key outcome of Task 2 was the testing completion of specially instrumented commercial DGS in the GE-SwRI Apollo sCO 2 compressor (27,000 rpm). Test data from the DGS showed elevated temperatures upwards of 350 o F, which are close to the higher operating temperature limit of the DGS. The temperature measurements provide insight into the expected thermal loads on DGS operating in high-speed sCO 2 compressor and provided test data for validation of an in-house thermal model of the compressor/seal. Under Task 2.0, this report also presents the development of a steady-state conjugate heat-transfer model of the DGS operating in the sCO 2 compressor – a first of its kind model for modeling heat transfer of sCO 2 in an actual operating compressor. The findings of the thermal model show a reasonable match between temperature predictions of the model and the measured temperature data, also pointing out the validity of the approach and assumptions made in modeling the flows, heat transfer coefficients and windage modeling in the rig. Under Task 3.0, this report presents the preliminary design of a large-diameter hybrid face seal (14 inch and 26-inch diameter) for field testing in a land-based GE turbine. The preliminary seal design effort presented in this report under project DE-FE0031924 builds on the development and successful laboratory testing for such large diameter hybrid face seal under the prior DE-FE0024007 project. Key aspects of seal fluid analyses with CFD, mechanical design considerations and assembly considerations in a land-based turbine are presented. Finally, under Task 3.0, this report also presents the continued high-temperature testing of the 14-inch diameter hybrid face seal developed previously under the DE-FE0024007 program. Specifically, test data demonstrating successful non-contact seal operation and seal effective leakage of 0.001-inch with seal inlet temperatures above 700 o F are presented in this report. Successful hybrid seal operation in a laboratory environment for a large diameter (14-inch) seal at temperatures above 700 o F is a major technological milestone for this technology.

01 COAL, LIGNITE, AND PEAT↗

Validation of RANS Models for Predicting Bypass Flow in Prismatic VHTRs

In prismatic Very-High Temperature Reactors (VHTRs), gaps occur between neighboring fuel assemblies referred to as bypass gaps. This work aims to measure the flow distribution between the bypass gaps between neighboring assemblies and the coolant channels within the assemblies. Computational fluid dynamics (CFD) was done using Reynolds-Averaged Navier Stokes (RANS) modeling in STAR-CCM+ to assess the ability of RANS models to accurately capture this phenomena. The RANS results were compared to experimental worked performed by Seoul National University (SNU) and Korea Atomic Energy Research Institute (KAERI) for validation. The CFD results showed good comparison with the experimental data available from SNU and flow split was characterized for multiple configurations, inlet mass flow rates, and bypass gap sizes.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Monthly Nitrate Concentration in a Karst Spring with and without Discrete Conduit Flow

Understanding the coupled continuum pipe-flow framework for modeling contaminant transport in karst systems is critical for protecting water resources therein. This study simulated point and non-point source nitrate concentration in a karst spring and investigated the results generated from the flow and transport model with and without discrete conduit flow. CFPv2 and CMT3D models were integrated to address the changes in nitrate concentration at a monthly scale in a karst spring, and the results were compared with that from an equivalent porous media (EPM) model with high-hydraulic conductivity (K) zones set in the main karstified area to represent conduits. The results show that the CFPv2+CMT3D model is able to describe well the recession of nitrate concentration in spring discharge, and the relatively larger deviation (slower nitrate recession) from the observed trend for the EPM model is probably a result of the limitation of utilizing high hydraulic conductivity cells to represent conduit. Moreover, simulated hydraulic heads in poorly karstified areas from the two models both show slight differences from the observations (the head RMSE values of calibration/validation for CFPv2 and MODFLOW models are 0.16 m/0.25 m and 0.26 m/0.17 m, respectively), indicating the inclusion of conduits may not affect the simulation considerably, and the lower the proportion of karstic area, the slight effects brought from the inclusion of conduits in the model. For highly karstified areas, the CFPv2+CMT3D model could provide more accurate results (head RMSE of calibration/validation for CFPv2 and MODFLOW are 0.22 m/0.06 and 0.52 m/0.47 m, respectively), showing the coupled continuum pipe-flow framework may be more appropriate for applying to highly and maturely karstified areas where the variations in the behavior of flow and contaminant transport are more affected by turbulent flow regime.

54 ENVIRONMENTAL SCIENCES↗

Modeling gas breakthrough and flow phenomena through engineered barrier systems using a discrete fracture approach

Compacted clays are being considered to build engineered barrier systems (EBS) intended for the safe isolation of high-level nuclear waste (HLW) and spent nuclear fuel (SNF). The corrosion of the metallic canister containing the HLW/SNF will lead to the generation and buildup of the gas pressure in the more internal part of the clay buffer. This phenomenon would eventually trigger the formation and propagation of fractures in the clay barrier, jeopardizing its safety functions. Here, in this work, we propose to use the fragmentation technique (MFT) to model evolving fractures in clays triggered by gas pressurization. The MFT has been successfully used to model the formation of fractures in concrete, drying cracks in soil, hydraulic and thermo-fractures in rocks. In this work, we extend the MFT to deal with multiphase fluid flow in deformable porous media, and we upgraded a fully coupled computer finite element code using the extended technique. The proposed approach is first verified against analytical solutions and is then applied to model gas breakthrough experiments in clays. A very satisfactory performance of the method is observed in all the analyses, showing the potential of the MFT to tackle multiphase flow problems in deformable porous media with evolving discontinuities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon flow through energycane agroecosystems established post-intensive agriculture

As part of an integrated energy and climate system, biomass production for bioenergy based on the tropical perennial C4 grass energycane can both offset fossil fuels and store soil carbon (C). We measured energycane yields, root biomass, soil C pools, and soil C stocks in a 4 year field trial and modeled C flow from plants to soils in the surface layer of no-till energycane planted after more than a century of intensive sugarcane agriculture. Aboveground yields ranged from 16.7 to 19.0 Mg C/ha over the 4 year trial. Although total C stocks did not significantly differ in the surface layer (approx. 0–20 cm) during the study, C in free and occluded light fractions decreased, whereas C in the mineral-rich dense fraction increased over 4 years. Belowground system inputs, estimated from measurements and informed by convergence in the final soil fraction model, were set to 2.5 Mg C ha -1 year -1 . With this input value, we estimated that surface soils retained photosynthetically fixed C predominantly within the mineral-associated organic matter pool for a mean and median transit time of 177 and 110 years, respectively. Although we did not model C flow to deep soil layers (approx. 0–100 cm), observed C accumulation (11.4 Mg C ha -1 year -1 ) and root growth down to 120 cm suggest that soil processes and resulting C sequestration at the surface are likely to persist deeper into the soil profile. Energycane, as a strong candidate for climate change mitigation and land degradation remediation, showed high biomass yields and allocation of resources to roots, with sequestered soil C expected to persist for over a century.

09 BIOMASS FUELS↗

Methane Emissions Produced by Pneumatic Devices and Produced Water Tanks on Natural Gas Wellsites

In recent years, more scrutiny has been placed on the release of greenhouse gasses (GHGs) by the oil and gas industry in North America as climate concerns increase. In the coming years, stricter regulations proposed by the Environmental Protection Agency (EPA) concerning the release of GHGs will further direct attention to GHG emissions quantification, modeling, and overall reduction. Active natural gas wellsites employ production equipment that routinely emits methane during normal operation. Production equipment often includes, but is not limited to, gas production units (GPUs), pneumatic controllers (PCs) and actuators powered by produced natural gas, and produced water storage tanks. This work details direct methane emission measurements taken at active natural gas wellsites originating from intermittent pneumatic controllers and produced water storage tanks. Efforts to model methane flows originating from these sources are also detailed, and those results are compared to direct measurement. The direct methane emissions measurement campaign collected flow data from two active dry wellsites in the Marcellus shale region. These measurement campaigns captured, in total, several weeks of continuous data from these sites. Data included multiple pneumatic device vent flowrate channels, a single channel monitoring total emissions vented from the produced water storage tanks, and corresponding site weather data. Periodic composition measurements were taken from flows vented from produced water tanks. Computational modeling efforts attempted to mimic a given site’s methane emissions behavior given limited input parameters. Direct measurement data and simulation data yielded, through reduction and analysis, apparent methane emissions factors, statistics on pneumatic device actuation events, important relationships between GPU parameters and pneumatic device actuation events, produced water storage tank vent flows and flow compositions, and parameters most contributing to the magnitude and frequency of these fugitive methane emissions. Measurements, considering data from a 72hr period, produced PC emissions factors from 0.008 to 0.03 SCFH per device with an average device emitting 0.026 SCFH at the first site. The second site produced PC emissions factors spanning from 0.05 to 25.47 SCFH with an average device emitting 1.93 SCFH not including the 25.47 SCFH device, which was assumed to be malfunctioning. Two unique water tank emission measurements were made, 618 g CH4 per hour and 32.6 g CH4 per hour. On a water basis, the water tank emissions factors were 0.66 and 0.042 kg CH4 per bbl of water, respectively.

03 NATURAL GAS↗

A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application

In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.

54 ENVIRONMENTAL SCIENCES↗

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗

Dual Impacts of Space Heating Electrification and Climate Change Increase Uncertainties in Peak Load Behavior and Grid Capacity Requirements in Texas

Around 60% of households in Texas currently rely on electricity for space heating. As decarbonization efforts increase, non‐electrified households could adopt electric heat pumps, significantly increasing peak (highest) electricity demand in winter. Simultaneously, anthropogenic climate change is expected to increase temperatures, the potential for summer heat waves, and associated electricity demand for cooling. Uncertainty regarding the timing and magnitude of these concurrent changes raises questions about how they will jointly affect the seasonality of peak demand, firm capacity requirements, and grid reliability. This study investigates the net effects of residential space heating electrification and climate change on long‐term demand patterns and load shedding potential, using climate change projections, a predictive load model, and a direct current optimal power flow (DCOPF) model of the Texas grid. Results show that full electrification of residential space heating by replacing existing fossil fuel use with higher efficiency heat pumps could significantly improve reliability under hotter futures. Less efficient heat pumps may result in more severe winter peaking events and increased reliability risks. As heating electrification intensifies, system planners will need to balance the potential for greater resource adequacy risk caused by shifts in seasonal peaking behavior alongside the benefits (improved efficiency and reductions in emissions).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Power Flow Spectroscopy Diagnostics & Platform Development at the Z Pulsed Power Facility

Large pulsed power accelerators deliver multi-MJ pulses of electrical energy to a variety of high energy density (HED) physics experiments that support stockpile science programs. Understanding the plasma formation mechanisms and resulting electrical power transport (or "power flow") in the vacuum magnetically insulated transmission lines (MITLs) is an important area of ongoing research, and could provide a means to improve the performance of today's pulsed power accelerators while improving confidence in the design options for next-generation pulsed power concepts. Power flow science has been studied for decades, but these studies have not provided a predictive understanding of plasma formation and expansion in MITL systems. Several recent factors in pulsed power system design have generated a renewed (and urgent) interest in developing validated, multi-physics power flow engineering models with increased scrutiny and understanding. Examples of these factors include (i) the use of high inductance experimental configurations that could increase current "loss", (ii) interest in long-pulse applications that require predictable pulse shapes, and (iii) the desire to develop a deeper understanding of how current loss phenomena scale to larger accelerator configurations. This work is directed to support the validation of multi-physics power flow engineering models required to realize pulsed power systems for the NNSA mission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗