Search NASA⌕ Search

SEARCH · Search NASA

Results for “CALCULUS”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗
DOE OSTI2024

Matching Curved Lattices to Anisotropic Tangent Planes

Radial quantization would be the ideal formalism for studying strongly-coupled near-conformal quantum field theories but it requires the ability to perform lattice calculations on static, curved manifolds, specifically a very long cylinder whose cross section is a sphere. Smoothly discretizing the surface of a sphere requires a graph with unequal edge lengths. The geometry of such graphs is well understood since 1961 using Regge Calculus. But, lattice quantum field theories are defined in terms of couplings which appear in the action rather than edge lengths and so the relationship between couplings and lengths must be determined dynamically. A simple example is computing the ratio of spatial to temporal lattice spacings in anisotropic lattice QCD. I will discuss our conjecture that computing anisotropic lattice spacing ratios on affine transformations of regular flat lattices is sufficient to determine coupling assignments on curved lattices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Lectures on Lie Group Analysis: Solving Differential Equations Using Symmetries

These notes are meant to be a supplemental reference for the beginner Lie Group Analyst. It is assumed that the reader has a basic concept of the fundamentals of Lie Group Theory (LGT), e.g. has seen the derivation of the infinitesimal generator and understands the mathematical meaning behind invariance. An excellent reference is Albright et al., “Symmetry Analysis of Differential Equations: A Primer,”. The reader is urged to read at least the first three chapters of that document to be able to follow the outset of Chapter 2 of this document. The reader should also have a general understanding of calculus, ordinary differential equations, and partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Clean Energy Deployment Baseline for the Energy Community and Low-Income Tax Credit Bonuses [Slides]

The Inflation Reduction Act of 2022 introduced, for the first time, place-based federal tax incentives for projects sited in “Energy Communities,” potentially changing the economic calculus of where projects are best sited. Storage projects can qualify for a 10-percentage-point bonus to the Investment Tax Credit (e.g., from 30% to 40%), while wind and solar projects may qualify for either the ITC bonus or a 10% bonus to the Production Tax Credit (e.g., from $\$27.5$ to $\$30.25$/MWh). Energy Communities are areas with historical ties to fossil fuel industries and above average unemployment levels (FFEU), with closed coal mines or power plants, or contaminated properties. They seek to identify locations across the US that could especially benefit from economic revitalization. This report explores how the new federal tax credit incentives are impacting clean energy deployment patterns and establishes historical baselines against which future changes can be compared. We include a few case studies of clean energy projects going specifically to areas that were recently impacted by coal power plant closures to provide concrete examples of investments in Energy Communities. However, this publication does not assess how much of the incentive benefits pass from clean energy developers to hosting communities, nor does it offer a comprehensive view of the economic effects of clean energy deployment on Energy Communities. Key highlights include: - As clean energy projects take multiple years to conceptualize and develop, it is likely too early to see shifts towards Energy Community locations either among newly built projects or those that entered interconnection queues in 2023. - Approximately 35% of onshore wind, 50% of solar, and 60% of storage capacity built in 2023 and the first half of 2024 are located in Energy Communities, making them likely eligible for bonus incentives. While these bonus incentives were not available to projects coming online before 2023, we used 2023 Energy Community definitions to classify whether past projects were built in what is now considered an Energy Community. The deployment levels for 2023-2024 are similar to recent years (2020-2022) for solar and storage but slightly lower for wind. - Clean energy capacity has surged in the interconnection queues over the last few years, with about 45-50% of both recently proposed and total queued capacity being located in Energy Communities. While the amount of capacity in Energy Communities has also grown, its relative share is either stable (solar and storage) or slightly lower (wind) among projects that entered the queue in 2023. - Clean energy projects can be built at lower costs in Energy Communities. The levelized cost of energy after incentives was on average $\$9$/MWh (24%) lower for solar projects and $\$2$/MWh (6%) lower for wind projects built in 2023, relative to projects not located in Energy Communities. Wholesale electricity values at Energy Community locations relative to the rest of the market vary by region. The average value was often higher for wind projects (-$\$3$ to $\$11$/MWh) but lower for solar projects (-$\$6$ to 0/MWh). - Distributed solar that is owned by commercial entities is eligible for the Energy Community bonus and also, potentially, a Low-Income Community bonus. Residential solar installations in qualifying Energy Communities that are third-party owned represent about 10% of the total residential market. Larger commercial and industrial solar installations in Energy Communities make up 17% of the total market in 2023. Nearly 2 GW of distributed solar was built in areas qualifying as Low-Income Communities in 2023, exceeding the available annual program cap of 700 MW. Continued tracking of these trends will be important for system planners, investors, and local communities.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Topological symmetry in quantum field theory

We introduce a definition and framework for internal topological symmetries in quantum field theory, including “noninvertible symmetries” and “categorical symmetries”. We outline a calculus of topological defects which takes advantage of well-developed theorems and techniques in topological field theory. Our discussion focuses on finite symmetries, and we give indications for a generalization to other symmetries. We treat quotients and quotient defects (often called “gauging” and “condensation defects”), finite electromagnetic duality, and duality defects, among other topics. We include an appendix on finite homotopy theories, which are often used to encode finite symmetries and for which computations can be carried out using methods of algebraic topology. Throughout we emphasize exposition and examples over a detailed technical treatment.

Mathematics↗
DOE OSTI2026

QPatLib v1.0 — Measurement-based quantum simulation Pauli string unitary pattern collections

This Zenodo record accompanies the paper “Scalable Measurement-Based Quantum Simulation Patterns for Benchmarking” arXiv.2605.12502 and provides QPatLib v1.0 measurement-pattern datasets in human-readable JSONL together with a ZIP archive of OpenQASM 3.0 circuits used for validation and reproducibility. The patterns and circuits implement Pauli string unitaries for benchmark cases. Cases include all possible string combinations for less than 6 qubits and strings used in Hamiltonians for certain diatomic molecules for 6 or more qubits. Format: Each pattern_*.jsonl file is containins measurement patterns for all subsets for a given model/instance and subset strategy: it begins with a preamble containing model metadata, subset definitions, provenance, and (when feasible) full-pattern test results, followed by one pattern entry per subset. Each subset entry includes a required pattern_ascii field storing the measurement pattern in the measurement-calculus/Graphix standard with signal shifting, written left-to-right in the canonical order nodes → edges → measurements (with signal dependencies) → byproduct corrections (X/Z). The circuits are included as circuit_files.zip. Patterns in this record were validated against the corresponding circuits and checked for causal flow. Codes for generating these patterns can be found at QPatLib repository on Github

Graphix↗
NASA NTRS1962

Interplanetary Trajectory Optimization with Powerlimited Propulsion Systems

A trajectory-optimization process is described in which the optimum­ thrust equations are derived using the calculus of variations. The mag­nitude of the thrust is constrained within an upper and a lower bound, but the thrust direction is arbitrary. This formulation allows both the constant-thrust program and the variable-thrust program to be con­sidered. For the constant-thrust program, certain propulsion-system parameters are optimized for maximum final vehicle mass. This theory has been used to study interplanetary missions to Venus and Mars using a power-limited propulsion system. Both one-way and round­ trip rendezvous trajectories are considered. The analysis employs a two-body inverse-square force-field model of three dimensions. An iterative routine used to solve the two-point boundary-value problem is described in the Appendix.

TRAJECTORY↗
NASA NTRS1962

Optimum Interplanetary Rendezvous Trajectories With Powerlimited Vehicles

The optimum-thrust equations for both variable and constant thrust are presented. These thrust programs are used to generate rendezvous trajectories from the Earth to Mars for various flight times and launch dates during the years 1968-71. The manner in which the propulsion requirements vary with flight time and launch date are considered, and a comparison of vehicle performance using the variable- and constant-thrust programs is presented. The optimization of the pro- pulsion system parameters is discussed, and the existence of optimum launch dates is interpreted in terms of certain transversality conditions derivable from the calculus of variations. A brief comparison of the advanced propulsion vehicle and the ballistic vehicle propulsion requirements is made for Earth-Mars rendezvous trajectories. An appendix considering the analytical basis for this work is included.

INTERPLANETARY TRAJECTORY↗
NASA NTRS1965

Necessary conditions for a multistage bolza- mayer problem involving control variables and having inequality and finite equation constraints

A multiplier rule and analogues of the Weierstrass and Clebsch conditions are developed for a multistage Bolza-Meyer calculus of variations problems. The number of stages is fixed, but partition points defining state boundaries are variable. Discontinuities are allowed in variables finite equations and inequalities, as well as differential equations, all of which involve control variables. An appendix summarizes some of the results obtained by C. H. Denbow, as modified by R. W. hunt, for a generalized Bolza problem. The appendix is independent of the rest of the paper.

Differential equation↗