Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

BETO 2021 Peer Review - CAP Process Research

Combined Algal Processing (CAP) is an algae biorefining approach that employs pretreatment and extraction operations to fractionate algal biomass into an organic liquid phase containing lipids, an aqueous hydrolysate phase containing carbohydrates and protein, and a residual solid phase containing insoluble matter. Each of these fractions is upgraded to fuels and/or co-products, including non-isocyanate polyurethanes from unsaturated lipids, fuels from saturated lipids, fuel precursor carboxylic acids from the hydrolysate, and conductive carbons from the solid residue. Techno-economic analysis of multiple CAP configurations suggests that some configurations have a viable pathway to algal biofuels at $2.50/gallon of gasoline equivalent.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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↗
DOE OSTI2020

DEMONSTRATION OF MULTI-PHYSICS EVALUATION OF FUEL ROD BURST PROBABILITY WITH BURNUP EXTENSION UNDER LBLOCA CONDITIONS

Extending the fuel discharge burnup level, e.g., from the current limit of rod averaged discharge burnup limit of 62 GWD/MT to a proposed new limit of 75 GWD/MT, can provide significant economic benefits to the current fleet of operating light water reactors (LWRs). It allows for longer operating cycles and improved resource utilization. The major economic gain of longer operating cycles is attributable to the increased capacity factor resulting from decreased refueling time as a fraction of total operating time, as well as fewer assemblies to be discharged for a given amount of energy produced. The main licensing challenges for higher burnup fuel are to ensure fuel rod safety under design basis accident conditions, especially under large-break loss-of-coolant accident (LBLOCA) and reactivity insertion accident (RIA). In this work, two-year cycle core design for a typical 4-loop pressurized water reactor (PWR) is performed with enrichment increased up to 6% and burnup extended to 75 GWD/MT. The fuel rod burst potential evaluations under large-break loss-of-coolant accident (LBLOCA) conditions are subsequently performed using the multi-physics best estimate plus uncertainty analysis framework LOTUS (LOCA Toolkit for the U.S. LWRs) and the preliminary results are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Methane foam performance evaluation in water-wet unconsolidated porous media: A systematic experimental investigation at elevated pressure and temperature conditions

In the petroleum industry, enhanced oil recovery (EOR) techniques employ foam extensively to establish conformance control in heterogeneous and fractured reservoirs in order to increase the sweep efficiency. In such applications, foam performance evaluation under complex subsurface conditions is pivotal for the effective and optimized deployment of foam treatment. However, there is a scarcity of hydrocarbon gas foam generation and evaluation studies that examine the relationships between foam performance and the critical foam parameters at high-pressure and high-temperature conditions. Here, this study aims at methodically investigating the effects of several foam parameters on methane foam performance in water-wet proppant packs under harsh operating conditions. This is in relation to the technical needs of hydrocarbon foam injection into hydraulically-induced, propped fractures in unconventional oil reservoirs. To this end, a state-of-the-art experimental foam generation apparatus was designed, fabricated, and commissioned. We performed a large number of foam flow experiments on proppant packs using methane gas and different foaming agents at 3500 psi and 115°C. Anionic and amphoteric surfactants were employed to probe the effect of their ionic nature on foam performance. Foam performance sensitivities to various foam generation parameters and operating conditions, such as surfactant concentration, gas fraction, total injection rate, operating pressure, salinity, and proppant pack length were investigated. To this end, steady-state pressure drops across the proppant packs during foam generation and foam's apparent viscosity were measured to quantify the foam performance of surfactants. The results were then analyzed to determine optimum values of the foam parameters and the interplay between these parameters are discussed here. The systematic results achieved from this work are in agreement with the trends available in the literature and provide new insights into complexities of in situ foam generation in water-wet, unconsolidated porous media at extreme reservoir conditions.

03 NATURAL GAS↗
DOE OSTI2020

One-Step Ahead Prediction of Thermal Mixing Tee Sensors with Long Short Term Memory (LSTM) Neural Networks

High-temperature advanced reactors under development, such as sodium fast reactors (SFR) and molten salt cooled reactors (MSCR), are expected to offer lower levelized cost of energy (LCOE) compared to existing light water reactor (LWR’s). In the existing light water reactors (LWR’s), operation and maintenance (O&M) expenses constitute the largest fraction of the total operating cost. Some of the O&M costs are related maintenance of sensors which can fail due to exposure to harsh environment in a reactor. The O&M costs of Advanced Reactor (AR)’s are expected to constitute a significant fraction of the total cost as well, because of high temperature and radiation level in AR are likely to cause material fatigue and premature failure of sensors and components. The O&M costs in AR’s could be reduced through integration of advanced informatics of performance-related sensors into a digital twin designed for reactor monitoring. For example, machine learning (ML) could be employed for real-time validation and correction of performance-related sensors, and reducing the number of performance-related physical sensor units through virtual sensing. As part of the effort, we investigate real-time validation of thermal hydraulic sensors through one-step ahead forecasting of sensor values using long short-term memory (LSTM) recurrent neural networks (RNN). The sensors are installed in a flow loop containing a thermal mixing tee, which is a common experimental model to study thermal fatigue in a thermal hydraulic loop. In addition, nonlinear transients generated in a thermal mixing tee constitute a good challenge data set for training and validation of ML algorithms. Sensors in this study include thermocouples, flow meters, and optical fibers for distributed temperature sensing. In one experiment, measurement data sets were obtained for a loop was filled with water, and in another experiment, measurements were performed on a loop filled with liquid metal Galinstan. We have also conducted preliminary investigation of one-step ahead prediction of fiber optics-based distributed temperature sensing with LSTM networks. In predicting fiber-based temperature measurements, we treated each gauge pitch of the fiber as an independent sensor. Accuracy of one-step ahead forecasting was estimated by calculating root mean square error (RMSE) for the test segment of time series of each sensor. RMSE’s for temperature sensors in water loop were, for the most part, lower than for the same sensors in Galinstan loop. The RMSE’s for flow meters were similar for both loops. The RMSE’s for distributed temperature measured with the fiber optic sensor were similar to those of the point sensors. Results of this study demonstrated the capability of LSTM one-step ahead forecasting with RMSE comparable to uncertainty in sensor measurements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Learning nonlinear operators via DeepONet based on the universal approximation theorem of operators

It is widely known that neural networks (NNs) are universal approximators of continuous functions. However, a less known but powerful result is that a NN with a single hidden layer can accurately approximate any nonlinear continuous operator. This universal approximation theorem of operators is suggestive of the structure and potential of deep neural networks (DNNs) in learning continuous operators or complex systems from streams of scattered data. Here, in this work, we thus extend this theorem to DNNs. We design a new network with small generalization error, the deep operator network (DeepONet), which consists of a DNN for encoding the discrete input function space (branch net) and another DNN for encoding the domain of the output functions (trunk net). We demonstrate that DeepONet can learn various explicit operators, such as integrals and fractional Laplacians, as well as implicit operators that represent deterministic and stochastic differential equations. We study different formulations of the input function space and its effect on the generalization error for 16 different diverse applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Real-time monitoring of the orthohydrogen fraction in a liquid hydrogen moderator

Here, we present the first direct measurements of the orthohydrogen fraction in an operating liquid hydrogen moderator loop at a high-powered neutron source. We have developed a Raman immersion probe capable of characterizing the orthohydrogen and parahydrogen makeup of a liquid hydrogen loop approximately every ten seconds, which can see orthohydrogen fractions lower than 0.20 without difficulty. At the Spallation Neutron Source, where the time-averaged power deposition in the studied moderator system is around 0.02 W/cm 3 , the orthohydrogen fraction increases with beam power and equilibrates at levels of 0.40 and responds to beam trips in a dynamic fashion. Additionally, we find that without deliberate catalyst systems in place, the relaxation due to the piping in the system has insignificant effect in catalyzing the relaxation of freshly condensed hydrogen toward thermodynamic equilibrium, and that the radiation levels associated with activation in a non-operating spallation neutron source target-moderator-reflector assembly likewise has insignificant effect on the liquid hydrogen kinetics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Overview of recent results from the ST40 compact high-field spherical tokamak

ST40 is a compact, high-field ($B_{\mathrm{T0}}\unicode{x2A7D} 2.1\,\,\,\textrm{T}$) spherical tokamak (ST) with a mission to expand the physics and technology basis for the ST route to commercial fusion. The ST40 research programme covers confinement and stability; solenoid-free start-up; high-performance operating scenarios; and plasma exhaust. In 2022, ST40 obtained central deuterium ion temperatures of $9.6 \pm 0.4\ \textrm{keV}$, demonstrating for the first time that pilot plant relevant ion temperatures can be reached in a compact, high-field ST. Analysis of these high-ion temperature plasmas is presented, including a summary of confinement, transport and microstability characteristics, and energetic particle instabilities. Recent scenario development activities have focused on establishing diverted H-mode plasmas across a range of toroidal fields and plasma currents, along with scenarios with high non-inductive current fractions. In future operations, beginning in 2025, a 1 MW dual frequency (104/137 GHz) electron cyclotron (EC) system will be installed to enable the study of EC and electron Bernstein wave plasma start-up and current drive. Predictive modelling of the potential performance of these systems is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Towards a Unified theory of Fractional and Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green’s identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green’s identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. Here, the proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Unified Theory of Fractional Nonlocal and Weighted Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green's identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green's identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. The proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Autocalibration of the E3SM Version 2 Atmosphere Model Using a PCA-Based Surrogate for Spatial Fields

Global Climate Model tuning (calibration) is a tedious and time-consuming process, with high-dimensional input and output fields. Experts typically tune by iteratively running climate simulations with hand-picked values of tuning parameters. Many, in both the statistical and climate literature, have proposed alternative calibration methods, but most are impractical or difficult to implement. We present a practical, robust, and rigorous calibration approach on the atmosphere-only model of the Department of Energy's Energy Exascale Earth System Model (E3SM) version 2. Our approach can be summarized into two main parts: (a) the training of a surrogate that predicts E3SM output in a fraction of the time compared to running E3SM, and (b) gradient-based parameter optimization. To train the surrogate, we generate a set of designed ensemble runs that span our input parameter space and use polynomial chaos expansions on a reduced output space to fit the E3SM output. We use this surrogate in an optimization scheme to identify values of the input parameters for which our model best matches gridded spatial fields of climate observations. To validate our choice of parameters, we run E3SMv2 with the optimal parameter values and compare prediction results to expertly-tuned simulations across 45 different output fields. This flexible, robust, and automated approach is straightforward to implement, and we demonstrate that the resulting model output matches present day climate observations as well or better than the corresponding output from expert tuned parameter values, while considering high-dimensional output and operating in a fraction of the time.

54 ENVIRONMENTAL SCIENCES↗