Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Higher-derivative corrections to entropy and the weak gravity conjecture in Anti-de Sitter space

We compute the four-derivative corrections to the geometry, extremality bound, and thermodynamic quantities of AdS-Reissner-Nordström black holes for general dimensions and horizon geometries. We confirm the universal relationship between the extremality shift at fixed charge and the shift of the microcanonical entropy, and discuss the consequences of this relation for the Weak Gravity Conjecture in AdS. The thermodynamic corrections are calculated using two different methods: first by explicitly solving the higher-derivative equations of motion and second, by evaluating the higher-derivative Euclidean on-shell action on the leading-order solution. In both cases we find agreement, up to the addition of a Casimir energy in odd dimensions. We derive the bounds on the four-derivative Wilson coefficients implied by the conjectured positivity of the leading corrections to the microcanonical entropy of thermodynamically stable black holes. These include the requirement that the coefficient of Riemann-squared is positive, meaning that the positivity of the entropy shift is related to the condition that c – a is positive in the dual CFT. We discuss implications for the deviation of η/s from its universal value and a potential lower bound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Traversable wormhole solution with a background Kalb–Ramond field

Highlights: • Wormhole solutions of gravity in a Lorentz breaking symmetry scenario are found. • Light deflection by the wormhole in a Lorentz breaking symmetry scenario is calculated. • It is shown the influence of the vacuum expectation value of a Kalb–Ramond field in deflection of light. We obtain a static spherically symmetric wormhole solution due to the vacuum expectation value (vev) of a Kalb–Ramond field. The Kalb–Ramond (vev) is a background tensor field which produces a local Lorentz symmetry breaking (LSB) of spacetime. Considering a non-minimal coupling between the Kalb–Ramond (vev) and the Ricci tensor, we found an exact traversable wormhole solution sustained by matter sources with a negative isotropic pressure. The matter source satisfies the energy conditions at the throat for particular choices of the LSB parameter. Moreover, we employ the Gauss–Bonnet method to analyze the deflection angle of light in the weak limit approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Elucidating Phosphate and Cadmium Cosorption Mechanisms on Mineral Surfaces with Direct Spectroscopic and Modeling Evidence

The simultaneous sorption of cations and anions at the mineral–water interface can substantially alter their individual sorption characteristics; however, this phenomenon lacks a mechanistic understanding. Our study provides direct spectroscopic and modeling evidence of the molecular cosorption mechanisms of the cadmium ion (Cd 2+ ) and phosphate (P) on goethite and layered manganese (Mn) oxide of birnessite, through in situ attenuated total reflection Fourier-transform infrared (ATR-FTIR), P K-edge X-ray absorption near-edge structure (XANES) spectroscopy, and surface complexation modeling. Phosphate synergistically cosorbed with Cd on goethite predominantly through P-bridged ternary complexes (≡Fe–P–Cd) and electrostatic interactions at wide pH conditions. Likewise, P and Cd exhibited synergistic cosorption on birnessite by forming P-bridged ternary complexes (≡Mn–P–Cd) and weak competitive sorption at the layer edge sites. As pH and Cd loading increased, the surface P species transitioned from a binary complex to a ternary complex and/or Cd 3 (PO 4 ) 2 precipitate for both goethite and birnessite. Compared to that in solution at pH 8, the formation of Cd 3 (PO 4 ) 2 was inhibited by the presence of goethite and birnessite, ascribed to the specific adsorption of P and Cd, more pronounced in birnessite due to the stronger sorption of Cd at its vacant sites. Finally, the discovered cosorption mechanisms of P and Cd have important implications for understanding and predicting their mobility and availability in Cd-contaminated settings.

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

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
DOE OSTI2021

Ultrafast x-ray pump x-ray probe transient absorption spectroscopy: A computational study and proposed experiment probing core-valence electronic correlations in solvated complexes

Femtosecond X-ray pump X-ray probe experiments are currently possible at free electron lasers such as the linac coherent light source (LCLS), which opens new opportunities for studying solvated transition metal-ligand complexes. In order to make the most effective use of these kinds of experiments, it is necessary to determine which chemical properties an X-ray pump pulse will measure. In this paper, we have combined electron cascade calculations and excited-state time-dependent density functional theory calculations to predict the initial state prepared by an X-ray pump and the subsequent X-ray probe spectra at the Fe K-edge in the solvated model transition metal complex, K 4 Fe II (CN) 6 . We find several key spectral features that report on the ligand-field splitting, solvent-solute interactions, and the 3p and 3d electron interactions. We then show how these features could be measured in an experiment.

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

Digitalization mapping and assessment process supporting ION strategic transformation activities

The existing fleet of commercial nuclear power plants (NPPs) are an important asset in the nation’s portfolio of electrical generating resources. Their continued safe and reliable operation are critical to providing a large source of carbon-free electricity to power the nation’s economy. The United States Department of Energy’s (DOE) Light Water Reactor Sustainability (LWRS) Program develops the scientific bases, methods, and tools, for the continued safe and economical operation of the nation's commercial NPPs. The Plant Modernization Pathway within LWRS Program focuses on providing guidance to industry on the full-scale implementation of modernization solutions for NPPs that significantly reduce the technical and financial risks associated with modernization. This research is focused on helping the nuclear industry understand how to digitize and digitalize their NPPs so that they can design their modernization solutions to be scalable, sustainable, and integrated both laterally and horizontally within their organization. That is, this research creates a digital transformation in NPPs by reshaping cultural mindsets and by identifying business efficiencies. In partnership with industry, and using four previously established guiding principles for digitalization, this research supported NPP modernization through assessing readiness for digitalization as a means to achieve integrated operations for nuclear. Specifically, this research created an assessment to review an entire organization’s work processes to gather information about the digitalization health of the plant. The assessment tools were administered to plant employees, and the results were used to develop a digitalization plan. The survey assessment identified the optimal candidate processes that would most benefit from a digitalization initiative which were revealed through analytical frameworks. One analysis calculated mean digitalization health indicator scores for all endorsed activities which allowed the researchers to rank and color code the results for easy identification. Individual health indicator scores are also provided, should our industry partner wish to understand these findings according to their own organizational priorities, business considerations and desired end-state. The results were also analyzed from the perspective that organizations are comprised of different types of innovators (e.g., generators, optimizers, conceptualizers, and implementers), which differentially affects the organization’s ability to comprehend and adapt to change (i.e., opportunities to innovate). Understanding the relative composition of innovator types at an NPP allows them to gather insights into the strengths and weaknesses they have in innovating how work is performed. For the utility that partnered with this research team, the results showed that implementers make up the largest portion of respondents and conceptualizers the smallest portion. Knowing the proportion of innovator types gave this organization insights on how they can effectively implement their innovation solutions. Additionally, the results were analyzed from a technical, economic, and risk perspective to identify and quantify work reduction opportunities (WROs). Recognizing that not all cost-saving opportunities are the same, a Technical, Economic and Risk Assessment (TERA) was performed to evaluate WROs to identify areas of greatest potential and lowest risk. The key results from TERA included a digitalization opportunity score for each activity, and a calculation of potential cost savings. These two outputs formed the bases for calculating a priority index and rank for the activities/processes assessed. From the prioritization calculations, TERA can then help the utility 1) decide what digitalization priorities to invest money in implementing and then 2) calculates how much should be invested in the digitalization initiatives selected to achieve cost savings and/or an acceptable return on investment. Last, onsite interviews revealed several inefficiencies in the standard work processes that occur cross-departmentally that are due to the absence of digitized and digitalized processes. Examples of these include time spent scanning paper documents and then uploading the documents electronically, obtaining signatures, and searching for desired information. This represents a digital but not digitalized process. Over 15 opportunities to improve work processes were identified through this multi-method digitalization assessment. The various analytical assessments used (e.g., TERA, digitalization health indicator scores), as well as discussions with the utility partner, corroborated that all the opportunities identified had a strong potential to make work processes more efficient and to improve overall performance of the NPP.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Comparing the Effects of Side Chain Dipole–Dipole Interactions and Hydrogen Bonding on the Mechanical and Electrical Properties of Poly(3-hexylthiophene)

Poly(3-alkylthiophenes) are simple conjugated polymers with good electrical properties, but tend to be brittle, limiting their application. Here, we investigate how side chain modifications can improve their toughness without compromising their charge transport. Specifically, we compare a weak dipole–dipole interaction of ester groups with a stronger hydrogen-bonding interaction of the hydroxyl groups. Two copolymers with targeted 5 and 10 mol % ester-functionalized side chains were synthesized and hydrolyzed in the solid state to introduce hydroxyl groups. Both hydrolyzed polymers became insoluble in most organic solvents, providing a path to layer-by-layer solution processing. All polymers’ optical, thermal, structural, electrical, and mechanical properties were investigated. Optical spectroscopy revealed that functionalization at these low levels did not significantly affect aggregation or π–π stacking. Differential scanning calorimetry indicated reduced crystallinity in the functionalized copolymers compared to that of P3HT. Grazing-incidence wide-angle X-ray diffraction data showed that functionalization led to an increase in lamellar spacing without altering molecular orientation or π–π stacking. Mechanical testing highlighted that at the lower functionalization level, both ester and hydroxyl groups similarly enhanced toughness, suggesting that disorder, rather than the specific nature of the functional group, was responsible for these improvements. At the higher functionalization level, hydroxyl groups significantly increased the elastic modulus while also negatively impacting charge carrier mobility. This study suggests that with stronger hydrogen-bonding groups, only very small amounts can improve P3HT’s toughness. For a wider tunable range, it may be a better strategy to use weak interactions such as dipole–dipole interactions of ester groups.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Isotropic-to-Nematic Transition in Salt-Free Polyelectrolyte Coacervates from Coarse-Grained Simulations

Recent interest in complex coacervation between oppositely charged polyelectrolytes (PEs) has been fueled by its relevance to biology in the context of membraneless organelle formation within living cells. For PEs with limited flexibility (such as double-stranded DNA), theoretical treatments and recent experiments have reported the emergence of liquid crystalline order (LCO) within the resulting coacervate phases. In this work, we study the underlying physics of this phenomenon using coarse-grained molecular dynamics simulations of symmetric semiflexible-semiflexible and asymmetric semiflexible-flexible coacervates. By comparing coacervates with the corresponding semidilute solutions of neutral polymers, we demonstrate that the presence of Coulomb interactions in coacervates facilitates orientational ordering, in agreement with theoretical predictions. Quantitative comparisons between our simulations and theory indicate that, for asymmetric nematic coacervates, the strong orientational ordering of stiff polyanions induces a weak ordering of the flexible polycations-an effect that was not anticipated by available theoretical studies. Simulations reveal that, for nematic coacervates, the preferred orientation of the PE chains at the liquid-liquid coacervate-supernatant interface is parallel, and the alignment of semiflexible PEs is homogeneous. The results presented here provide new molecular-level insights into the intracoacervate LCO and will help motivate further experimental and theoretical activities in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation

We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗
DOE OSTI2024

Wormholes without averaging

After averaging over fermion couplings, SYK has a collective field description that sometimes has “wormhole” solutions. We study the fate of these wormholes when the couplings are fixed. Working mainly in a simple model, we find that the wormhole saddles persist, but that new saddles also appear elsewhere in the integration space — “half-wormholes.” The wormhole contributions depend only weakly on the specific choice of couplings, while the half-wormhole contributions are strongly sensitive. The half-wormholes are crucial for factorization of decoupled systems with fixed couplings, but they vanish after averaging, leaving the non-factorizing wormhole behind.

1/N Expansion↗
DOE OSTI2025

Quantifying the Lewis Acidity of Mono-, Di-, and Trivalent Cations in Anhydrous Bis(trifluoromethylsulfonyl)imide Salts

While the bis(trifluoromethylsulfonyl)imide anion (TFSI – ; formula [N(SO 2 CF 3 ) 2 ] – ) has been noted for its practical utility, the use of TFSI – salts as sources of Lewis acidic metal cations for studies of cation-driven tuning effects has not been reported. Here, the effective Lewis acidity of mono-, di-, and trivalent cations (namely, K + , Na + , Li + , Ba 2+ , Ca 2+ , Mg 2+ , Zn 2+ , La 3+ , Y 3+ , Lu 3+ , and Sc 3+ ) in the form of their TFSI – salts is described, along with quantitative comparisons to salts of several other weakly coordinating anions (namely, SO 3 CF 3 – , PF 6 – , and BArF 24 – ). Triphenylphosphine oxide (TPPO) was used as a 31 P NMR probe in titration experiments for quantification of key parameters describing the effective Lewis acidity of the salts in acetonitrile (CH 3 CN) solutions. Notably, the TFSI – salts of di- and trivalent cations were found to display strong binding to TPPO with significant speciation and were found to be quite hygroscopic. Taken together, the measurements demonstrate that TFSI – salts are systematically better/stronger effective Lewis acids than their triflate analogues. And, considering the excellent solubility of TFSI – salts, these materials appear attractive for further use and development in Lewis-acidity-dependent applications, including catalysis and tuning of multimetallic materials.

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

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2024

TCF - Sustainable Well Cement for Geothermal, Thermal Recovery, and Carbon Storage wells

Primary cementing is the most important operation performed on a subterranean well. Cement, placed in the annulus between the casing and the formations, serves as a hydraulic seal preventing fluids and gas migrations, protecting steel casing from corrosion, and supporting the well structure. Poor cementing jobs can be directly responsible for wells not reaching their full capacity, casing corrosion, compromised well integrity and, in the worst-case scenarios, well collapse. Geothermal, thermal recovery, and carbon storage wells offer environments that are especially difficult for cements to survive while the required lifespan of these wells can be years. In these wells currently used cements cannot provide durable well integrity, and new solutions are necessary. The necessity of stabilizing the electric grid, increasing its flexibility, and providing energy on demand will further expand the market of durable cement solutions for applications in high-temperature underground reservoir thermal energy storage wells (HT RTES). The design challenges of special cement systems for such wells originate from chemistry limitations of currently used well cements, aggressive environments, temperature and repeated shock conditions associated with them, and very weak formations that are not uncommon in these wells. During the last 70 years the most common systems for thermal-well construction have been Ordinary Portland Cement (OPC) (absolute majority of the wells), silica-lime system, and high aluminum cement-containing formulations. The major issues of calcium-silicate hydrates, that form during the hydration of OPC resulting in hardened material, are their poor chemical resistance to acids due to the calcium interactions with acid anions followed by its eventual dissolution, even in mild acids, such as carbonic acid from CO 2 dissolution, and inadequate bonding to the steel causing serious casing corrosion problems (Sugama & Pyatina, 2019). Performance of OPC-based cement may be significantly improved with organic additives; however, those are limited by their temperature stability. An alternative to calcium-silicate cement was developed by BNL and commercialized by Halliburton as ThermaLock TM cement. Calcium-aluminate-phosphate-based chemistry of the material allowed overcoming acid-resistance problems of OPC, especially in CO 2 -rich geothermal environments. In 2012-2015, building upon this experience, BNL developed Thermal Shock Resistant Cement (TSRC) with the support of Geothermal Technology Office (GTO) of DOE. TSRC has high-temperature stable chemistry, superior properties of cement-metal casing bond and corrosion protection (4-times better than currently used OPC-based formulation), significantly outperforms common well-cements in thermal-shock tests, has self-healing properties, and like ThermaLock TM is CO 2 resistant incorporating carbonate ions into stable hydrates (Gill et al., 2012).

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Ion-Specific Effects on PuO 2 Nanoparticle Aggregation and Dissolution in Concentrated Electrolytes

Hydrolytic PuO 2 nanoparticles (NPs) are a dominant aqueous Pu-bearing phase in high ionic strength nuclear wastes, yet their reactivity in nonideal brines remains poorly constrained. We quantify how electrolyte identity and concentration control PuO 2 NP aggregation and ligand-assisted dissolution in acidic, high salinity solutions (NaCl, NaNO 3 , NaClO 4 , Na 2 SO 4 , Na 2 C 2 O 4 up to 5 M). A multitechnique workflow combining liquid scintillation counting (operationally defined aqueous [Pu]), scattering/electrokinetic measurements (aggregate size and zeta potential), and spectroscopy (UV–vis, XPS, Raman) resolves electrolyte-dependent partitioning between colloidal and molecular Pu species. Weakly coordinating anions (ClO 4 – , Cl – , NO 3 – ) largely preserve the (aggregated) nanoparticulate fraction but generate distinct dissolved species at high concentration, i.e., Pu(IV)–nitrato complexes in NaNO 3 and Pu(VI)–chloro complexes in NaCl. In contrast, stronger ligands substantially perturb PuO 2 NP stability: sulfate promotes partial dissolution to Pu(IV)–sulfate complexes at low concentration but reduces aqueous [Pu] at higher sulfate levels via secondary Pu(IV) sulfate formation, whereas oxalate drives strong dissolution to aqueous Pu–oxalate complexes. Aged NPs show similar trends with reduced aqueous fractions and more dominant aggregation mechanisms. In conclusion, these spectroscopically constrained speciation data provide a foundation for incorporating PuO 2 NP reactivity into thermodynamic and reactive transport models for high salinity waste and brine environments.

Aggregation↗
DOE OSTI2022

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pricing and Energy Trading in Peer-to-Peer Zero Marginal-Cost Microgrids

Efforts to utilize 100% renewable energy in community microgrids require new approaches to energy markets and transactions to efficiently address periods of scarce energy supply. In this paper we contribute to the promising approach of peer-to-peer (P2P) energy trading in two main ways: analysis of a centralized, welfare-maximizing economic dispatch that characterizes optimal price and allocations, and a novel P2P system for negotiating energy trades that yields physically feasible and at least weakly Pareto-optimal outcomes. Our main results are 1) that optimal pricing is insufficient to induce agents with batteries to take optimal actions, 2) a novel P2P algorithm to addresses this while keeping private information, 3) a formal proof that this algorithm converges to the centralized solution in the case of two agents negotiating for a single period, and 4) numerical simulations of the P2P algorithm performance with up to 10 agents and 24 periods that show it converges on average to total welfare within 0.1% of the social optimum in on the order of 10s to 100s of iterations, increasing with the number of agents, time periods, and total storage capacity.

batteries↗