Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fractional operators”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparison of Tempered and Truncated Fractional Models

Tempered fractional operators are able to model effects that classical partial differential equations cannot quite capture. For example, this includes the super- and subdiffusion effects that are present in hydrology and geophysics models. However, tempered fractional operators are computationally intensive. We analyze a truncated variation of the fractional operators, which are less computationally intensive, in an effort to use them in place of the more complex tempered variation. In particular, we train parameters of the truncated operator using neural networks in order to optimize the difference of the actions of the two operators.

58 GEOSCIENCES↗
DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advances in the long-pulse steady-state high beta H-mode scenario with active controls of divertor heat and particle fluxes in EAST

Since the last IAEA-Fusion Energy Conference, the Experimental Advanced Superconducting Tokamak (EAST) research program has been, in support of ITER and CFETR, focused on development in terms of the long-pulse steady-state (fully noninductive) high beta H-mode scenario with active controls of the stationary and transient divertor heat and particle fluxes. The operational domain of the steady-state H-mode plasma scenario has been significantly extended with ITER-like tungsten mono-block divertor, plasma control and heating schemes. EAST has achieved several important milestones in the development of high β p H-mode scenario and its key physics and technologies. A 60 s-scale long-pulse steady-state high β p H-mode discharge with the major normalized plasma parameters similar to the designed performance of the CFETR 1GW fusion power operation scenario has been successfully established and sustained by pure RF heating and current drive. Several feedback control schemes have been developed for a sustained detachment with good core confinement. This includes control of the total radiation power, target electron temperature, and particle flux measured using divertor Langmuir probes or a combination of the control of target electron temperature and AXUV radiation near the X point. The detachment feedback control schemes have been integrated with small-ELM regimes and high β p scenario via neon seeding, enabling a core and edge compatible integrated high-beta scenario applicable to long-pulse operations. ELM suppression has been achieved using various methods, including resonant magnetic perturbations and impurity seeding. Full suppression of ELMs by using n=4 RMPs has been demonstrated for ITER for the first time in low input torque plasmas in EAST. Finally, EAST has been operated with helium to support the ITER research requirements for the first time. For a long-pulse, high bootstrap current fraction operation, a new lower tungsten divertor with active water-cooling has been installed, along with improvements in the heating and current drive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Variable-spectrum mode control of high poloidal beta discharges

DIII-D experiments demonstrate that high pressure, broad current profile equilibria can be accessed in the high poloidal beta regime by optimizing the MHD mode control poloidal spectrum. A novel, variable spectrum (VS) magnetic feedback scheme implemented using the DIII-D internal non-axisymmetric coils (I-coils) facilitated access to reduced internal inductance $l$ i operation above the no-wall beta limit compared with both no feedback and fixed spectrum feedback. In addition, the VS feedback helped avoid beta collapses caused by marginally unstable resistive wall mode activity. The lower and upper I-coil rows were configured in two independent feedback loops, allowing the feedback field's poloidal spectrum to vary and track changes in the plasma mode structure as the edge safety factor q 95 varied from 11 to 6 during the discharges. The q 95 dependence of the measured phase difference between the lower and upper I-coil rows during VS feedback is qualitatively compatible with ideal MHD simulations of the least-stable plasma kink mode and with plasma response simulations that included kinetic modifications to ideal MHD. The VS feedback approach is a straightforward way to improve resilience to variations in mode structure that occur as plasma parameters change. The demonstrated expansion of the operating space to lower $l$ i is expected to improve the coupling of the plasma kink mode to external fields and beneficial wall eddy currents, and is compatible with high bootstrap fraction operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simultaneous enhancement of tritium burn efficiency and fusion power with low-tritium spin-polarized fuel

This study demonstrates that using spin-polarized deuterium-tritium (D-T) fuel with more deuterium than tritium can increase tritium burn efficiency (TBE) by at least an order of magnitude without compromising fusion power output, compared to unpolarized fuel. Although previous studies show that a low tritium fraction can enhance TBE, this strategy resulted in reduced fusion power density. The surprising improvement in TBE at fixed power reported here is due to the TBE increasing nonlinearly with decreasing tritium fraction but the fusion power density increasing roughly linearly with D-T cross section. A study is performed for an ARC-like tokamak producing 481 MW of fusion power with unpolarized 53:47 D-T fuel, finding the minimum startup tritium inventory (I startup,min ) is 0.69 kg. By spin-polarizing half of the fuel and using a 60:40 D-T mix, I startup,min is reduced to 0.08 kg, and fully spin-polarizing the fuel with a 63:37 D-T mix further reduces I startup,min to 0.03 kg. Some ARC-like scenarios are predicted to achieve plasma ignition with relatively modest spin polarization. These findings indicate that, with advancements in helium divertor pumping efficiency, TBE values of approximately 10%–40% could be achieved using low-tritium-fraction and spin-polarized fuel with minimal power loss. This would dramatically lower tritium startup inventory requirements and reduce the amount of on-site tritium. More generally than just for spin-polarized fuels, increased plasma performance can be used to increase TBE. This strongly motivates the development of spin-polarized fuels and low-tritium-fraction operation for burning plasmas.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Full Scale HEPA Filter Encapsulation in Ultra-High-Performance Grout - Proof of Concept

The Hanford Tank Waste Treatment and Immobilization Plant (WTP) currently being constructed to treat radioactive waste, includes vitrification facilities for both the high-level waste (HLW) and low activity waste (LAW) fractions. Operation of the WTP will produce contaminated high-efficiency particulate air filters (HEPA), as part of the solid secondary waste (SSW) stream. The HEPA filters receive off-gas from the vessel vent header and primary off-gas treatment system in LAW Facility and remove particulate contaminants including 99 Tc and 129 I salts. The current disposal method for HEPA filters is encapsulation in metal containers using cementitious material (CM) and disposal in the Integrated Disposal Facility (IDF). Results from the 2017 IDF Performance Assessment (PA) WRPS (2018) demonstrated that while compliance is maintained for the 1000 year compliance period mandated by DOE O 435.1 and its accompanying manual, release of constituents from the HEPA filters result in exceedance of the performance objective imposed as the groundwater regulatory limit at later times. For example, at about 1500 years post-closure, solid secondary waste (SSW), including HEPA filters is predicted to become a dominant contributor to 99 Tc release and over the 10,000-year sensitivity analysis period, SSW is the dominant contributor of 129 I release to the groundwater. The estimated release could potentially be reduced if the HEPA filters are not compacted and waste containers could be distributed throughout a large space, thereby diluting the contaminant release. Additionally, a better cementitious material could be used to encapsulate the HEPA filters. One alternative method for disposal of contaminated HEPA filters is encapsulation of the filters in ultra high-performance grout (UHPG). UHPG is a variation of ultra-high-performance concrete (UHPC) is commonly used in the prestressed concrete industry for large structural members. Recent studies of UHPG show it has excellent properties for containing radionuclides such as 99 Tc and 129 I Nichols and Kaplan (2021). This report presents the results of the first attempt to encapsulate a clean, full-size HEPA filter in UHPG and evaluate the effectiveness of the immobilization process and final waste form. A full-scale proof-of-concept simulated waste form was prepared by encapsulating a HEPA filter in a 110-gallon stainless steel (SS) drum using UHPG. A change from Type I/II PC to Type 1L PLC was made after American Rock Products informed the team that they would no longer be using Type I/II by the end of 2024 and the northwest was phasing out Type I/II PC overall. Type I/II PC used in previous studies of UHPG for encapsulation (Nichols and Kaplan 2021). After the UHPG was cured both the scaled mockup and the full-scale simulated waste forms were sectioned for visual examination. UHPG completely encapsulated the filters and bonded to the external surfaces of materials comprising the filters. No cracks were observed in the sectioned waste forms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On Road vs. Off Road Low Load Cycle Comparison

Reducing criteria pollutants while reducing greenhouse gases is an active area of research for commercial on-road vehicles as well as for off-road machines. The heavy duty on-road sector has moved to reducing NOx by 82.5% compared to 2010 regulations while increasing the engine useful life from 435,000 to 650,000 miles by 2027 in the United States (US). An additional certification cycle, the Low Load Cycle (LLC), has been added focusing on part load operation having tight NOx emissions levels. In addition to NOx, the total CO2 emissions from the vehicle will also be reduced for various model years. The off-road market is following with a 90% NOx reduction target compared to Tier 4 Final for 130-560 kW engines along with greenhouse gas targets that are still being established. The off-road market will also need to certify with a Low Load Application Cycle (LLAC), a version of which was proposed for evaluation in 2021. Since the LLAC has not been finalized, this study is being conducted to compare and contrast the LLC for on-road with the LLAC for off-road as there might be some shared learnings. A US off-road production 2023 Fiat Powertrain 13L engine and aftertreatment system was chosen for this work. This engine is used in production for both off-road and on-road products, so it is a good choice for this study. The associated off-road aftertreatment system was aged for more relevant comparisons. The engine calibration was not altered for either of the low load cycles. This study shows that the cycles are quite different in nature as the market needs are different. The LLC includes a large fraction of operation at idle and lower speeds, representing products that use the engine primarily for motive power, where lower vehicle speed means a lower engine speed and load. The LLAC has more time and load spent at high speeds and slightly higher loads. The off-road products represented by this cycle often use the engine to drive auxiliary equipment which means higher parasitic loads and hand/fixed throttle. The comparison will include the use profiles, tailpipe NOx and greenhouse gas emissions (CO2, N2O).

McCarthy, James↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

We formulate the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that both the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we show that the adjoint wave propagator preserves the dispersion and compensates the amplitude, while the time-reversed adjoint wave propagator behaves identically as the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit Q parameterization, which avoids the implicit Q in the conventional viscoacoustic/viscoelastic full waveform inversion (Q-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. Here, we find that the lossless propagation kernel dominates the velocity kernel, while the dissipation kernel dominates the attenuation kernel over the dispersion kernel.

58 GEOSCIENCES↗
DOE OSTI2022

BETO 2021 Peer Review - FCIC Task 5 - Preprocessing

In the study of feedstock variability, this presentation specifically looks at preprocessing processes that will created well-defined and homogenous feedstock from variable biomass resources, in the quest to produce homogeneous intermediates that will be converted into market-ready products. The objective of this stage of the project is to develop science-based design and operation principles informed by TEA/LCA that result in predictable, reliable, and scalable performance of preprocessing unit operations (comminution, fractionation, deconstruction, and real-time imaging).

biomass↗
DOE OSTI2021

Connections between nonlocal operators: From vector calculus identities to a fractional Helmholtz decomposition

Nonlocal vector calculus, which is based on the nonlocal forms of gradient, divergence, and Laplace operators in multiple dimensions, has shown promising applications in fields such as hydrology, mechanics, and image processing. In this work, we study the analytical underpinnings of these operators. We rigorously treat compositions of nonlocal operators, prove nonlocal vector calculus identities, and connect weighted and unweighted variational frameworks. We combine these results to obtain a weighted fractional Helmholtz decomposition which is valid for sufficiently smooth vector fields. Our approach identifies the function spaces in which the stated identities and decompositions hold, providing a rigorous foundation to the nonlocal vector calculus identities that can serve as tools for nonlocal modeling in higher dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Self-consistent investigation of density fueling needs on ITER and CFETR utilizing the new Pellet Ablation Module

Abstract Self-consistent modeling using the stability, transport, equilibrium, and pedestal (STEP) workflow in the OMFIT integrated modeling framework (predicting pedestal with EPED, core profiles with TGYRO, current profile with ONETWO, and EFIT for equilibrium) suggests ITER and future devices such as China Fusion Engineering Test Reactor (CFETR) Zhuang et al (2019 Nucl. Fusion 59 112010) will benefit from high-density operation (Greenwald limit fraction f g w ≈ 0.7−1.3). Regimes with an operational density near the Greenwald limit will likely need peaked density profiles so that the pedestal density remains below the Greenwald limit. Peaked density profiles can be achieved with the help of pellet injection. A flexible Pellet Ablation Module (PAM), which predicts the density source based on a comprehensive analytical pellet ablation model, has been developed for predicting pellet fueling for transport studies, and has been incorporated into the STEP workflow for predictive modeling. This workflow is applied to DIII-D and finds good agreement with experiments. On ITER the effect of pellet fueling is examined in an advanced inductive scenario, where a fusion gain of up to Q = 9 is predicted with strong central pellet fueling. On CFETR, with a mid-radius density source, an average of 1.5 × 10 22 electrons s −1 are required to achieve the density and temperature profiles necessary for the 1000 MW advanced scenario with a tritium burn-up fraction of ∼ 3 % .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗
DOE OSTI2025

Gas-Gap Calorimeter Sizing and Rating Tools for Heat Pipe Experimentation

Effective gas-gap calorimeter sizing and rating tools are required to design calorimeters that enable well-defined operating conditions, high cooling powers, and accurate calorimetry. The present report describes two tools utilizing Microsoft Excel and MathWorks MATLAB that enable the sizing and rating of gas-gap calorimeters with a He-Ar binary gas mixture in the gap. The Excel tools are two easy-to-use spreadsheets for calculating heat pipe operating temperature or He-Ar mole fractions given the required inputs. The MATLAB tool includes similar models with more robust solution algorithms along with sub-options for full-vacuum, static gas-gap, and forced circulation in the gas-gap. The accuracy of the developed tools were demonstrated by comparing with existing work in literature. The MATLAB tool was used for a parametric study to determine the gas-gap thickness, coolant gap thickness, coolant flow rate, and coolant inlet temperatures for the testing of a 3/4 in outer diameter heat pipe up to powers of 10 kW and temperatures of 1,000°C. The effects of heat pipe and calorimeter inner shell surface emissivities were investigated, considering the inability to control or accurately measure surface emissivities in most applications. It was found that controlling the He-Ar mole fractions provides a wide operating range for gas-gap thicknesses around ~ 0.042--0.090 in (~ 1.07--2.29 mm) for a surface emissivity range of 0.4--0.8, since conduction heat transfer is a significant fraction of the heat transfer across the gas-gap. In addition, it was found that turbulent flow in the coolant gap is needed at high input powers to prevent shell temperatures from approaching the boiling temperature of the water coolant. The parametric study resulted in choosing stainless steel tubes with an outer diameter and thickness of 1 x 0.065 in as the inner shell, and a 1-1/2 x 0.156 in as the outer shell of the calorimeter. Overall, this report presents the necessary information for the sizing and rating of gas-gap calorimeters with He-Ar mixtures for the testing of high-temperature heat pipes (~ 500--1,000°C).

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗