Search NASA⌕ Search

SEARCH · Search NASA

Results for “interior point”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the Convergence of Inexact Predictor-Corrector Methods for Linear Programming

Interior point methods (IPMs) are a common approach for solving linear programs (LPs) with strong theoretical guarantees and solid empirical performance. The time complexity of these methods is dominated by the cost of solving a linear system of equations at each iteration. In common applications of linear programming, particularly in machine learning and scientific computing, the size of this linear system can become prohibitively large, requiring the use of iterative solvers, which provide an approximate solution to the linear system. However, approximately solving the linear system at each iteration of an IPM invalidates the theoretical guarantees of common IPM analyses. To remedy this, we theoretically and empirically analyze (slightly modified) predictor-corrector IPMs when using approximate linear solvers: our approach guarantees that, when certain conditions are satisfied, the number of IPM iterations does not increase and that the final solution remains feasible. We also provide practical instantiations of approximate linear solvers that satisfy these conditions for special classes of constraint matrices using randomized linear algebra.

Dexter, Gregory↗
DOE OSTI2024

Riemannian Optimization Applied to AC Optimal Power Flow

The nonlinear, nonconvex AC optimal power flow problem is of growing importance as the nature of the power grid evolves. This problem can be difficult to solve for interior point methods. However, the advent of optimization algorithms over smooth Riemannian manifolds presents an alternative approach. The nonlinear, nonconvex constraints in the AC power flow problem form an embedded submanifold of Euclidean space. In this paper, the authors explore the performance of Riemannian optimization algorithms for the ACOPF problem where the optimization is performed directly on the AC power flow manifold. This is done by using the Julia programming language and the Julia packages PowerModels.jl and Manopt.jl.

AC optimal power flow↗
DOE OSTI2020

Melting Curve and Phase Relations of Fe-Ni Alloys: Implications for the Earth's Core Composition

We discuss that nickel is the second most abundant element in the Earth's core. However, the properties of Fe-Ni alloys are still poorly constrained under planetary cores conditions, in particular concerning the effect of Ni on the melting curve of Fe. Here we show that Ni alloying up to 36 wt% does not affect the melting curve of Fe up to 100 GPa. However, Ni strongly modifies the hexagonal-closed-packed/face-centered-cubic (hcp/fcc) phase boundary, pushing the hcp/fcc/liquid triple point of Fe-20wt%Ni to higher pressures and temperatures. Our results allow constraining the triple point for Fe-10wt%Ni, a composition relevant for the Earth interior, and point out a decrease of the melting temperature at core-mantle boundary by 400 K with respect to pure Fe. A lower amount of light elements than previously predicted is thus required to reduce the crystallization temperature of core materials below that of a peridotitic lower mantle, in better agreement with geochemical observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gradient-Informed Design Optimization of Select Nuclear Systems

In this work, we present a gradient-informed design optimization of nuclear reactor core components based on neutronics objectives with both continuous and discrete materials. The main argument in favor of using gradient-informed design optimization is that it scales well with increasing dimensionality of the design space. First, a challenge problem with 121 free parameters is solved with a gradient-informed method and then with a genetic algorithm. Then, a challenge problem to optimize the flux profile of a simplified assembly with eight axial zones is solved. Both challenge problems are solved using directly calculated derivatives from Tools for Sensitivity and Uncertainty Analysis Methodology Implementation (TSUNAMI) in the SCALE package. Furthermore, we demonstrate how a discrete optimization problem—selection of materials for 121 voxels—can be lifted into a continuous problem with mixed materials. In the continuous space, adjoint-based gradients are well-defined, and gradient descent is applicable. Then, a forcing function is introduced that with the selection of an appropriately sized hyperparameter can be used to guide the optimized continuous solution back into a discrete solution. This paper presents an account of the challenges that were faced when applying a gradient-informed optimization algorithm using a Monte Carlo calculation to estimate the gradient information and compares a gradient descent optimization method to a genetic algorithm optimization of the same geometry. Overall, this work demonstrates the potential use of adjoint-based gradient calculations in design optimization of nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Six-point functions and collisions in the black hole interior

In the eternal AdS black hole geometry, we consider two signals sent from the boundaries into the black hole interior shared between the two asymptotic regions. We compute three different out-of-time-order six-point functions to quantify various properties of the collision of these signals behind the horizons: (i) We diagnose the strength of the collision by probing the two-signal state on a late time slice with boundary operators. (ii) We quantify two-sided operator growth, which provides a dual description of the signals meeting in the black hole interior, in terms of the quantum butterfly effect and quantum circuits. (iii) We consider an explicit coupling between the left and right CFTs to make the wormhole traversable and extract information about the collision product from behind the horizon. At a technical level, our results rely on the method of eikonal resummation to obtain the relevant gravitational contributions to Lorentzian six-point functions at all orders in the GN-expansion. We observe that such correlation functions display an intriguing factorization property. We corroborate these results with geodesic computations of six-point functions in two- and three-dimensional gravity.

2D Gravity↗
DOE OSTI2021

Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.

42 ENGINEERING↗
DOE OSTI2024

Nb 3 Sn Superconducting Cavities by Bronze Routes for Accelerator Stewardship

This project explored the challenges of making Nb 3 Sn superconducting radio-frequency (SRF) cavities by bronze routes. The SRF research community has prioritized Nb 3 Sn for many reasons, including the possibility to operate linacs for stewardship applications in continuous-wave mode with compact cryogenic systems. Existing Nb 3 Sn approaches use a cavity body made of niobium and apply tin to the polished cavity interior, a starting point that contains the main costs of present niobium cavity technology. Further complexity is then added because a high vacuum and temperature >1100 °C must be used to prevent formation of undesirable Nb-Sn phases, de-gas interstitial contaminants, and maintain the high purity of the Nb body. The project sought a more cost-effective solution for the same performance opportunities by exploiting the Cu-Sn-Nb diffusion reactions used to make Nb 3 Sn superconducting wires for magnets, which can be carried out at ~700 °C and are compatible with much cheaper copper cavity bodies. The project also sought understanding of the materials science and factors that affect superconducting properties. The project has impacted SRF cavity technology by demonstrating feasibility of alternative routes that enable compact, high-power electron accelerators for energy, security, environment, medical, and commercial applications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Reply to “Comment on ‘Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter' ”

In their comment to our paper [1], Bertsch and Stroberg [2] provide three criticisms. Two of these concern the interpretation of our dispersive optical model (DOM) results and their relation to the liquid drop model (LDM) parameters. The third criticism focuses on the potential systematic uncertainties on our DOM results associated with missing three-body contributions. Before addressing these critiques, we want to state that the key message of our paper remains whether or not the DOM results agree with the LDM predictions in the nuclear interior. Here, the key point is that the standard determination of the saturation energy from the LDM is not ideal since the total binding energy has a minimal contribution from the core of the nucleus as pointed out in Figs. 1–3 of our paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

High‐order symmetric cubature rules for tetrahedra and pyramids

Abstract In this article, we present an algorithm to construct high‐order fully symmetric cubature rules for tetrahedral and pyramidal elements, with positive weights and integration points that are in the interior of the domain. Cubature rules are fully symmetric if they are invariant to affine transformations of the domain. We divide the integration points into symmetry orbits where each orbit contains all the points generated by the permutation stars. These relations are represented by equality constraints. The construction of symmetric cubature rules require the solution of nonlinear polynomial equations with both inequality and equality constraints. For higher orders, we use an algorithm that consists of five sequential phases to produce the cubature rules. In the literature, symmetric numerical integration rules are available for the tetrahedron for orders p = 1 – 10, 14 , and for the pyramid up to p = 10 . We have obtained fully symmetric cubature rules for both of these elements up to order p = 20 . Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, weakly singular, and trigonometric test functions over both elements with flat and curved faces. With increase in p , improvements in accuracy is realized, though nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2025

Modeling Co2 Flow Through Faulted/Fractured Reservoirs Using Tedfm in Corner-Point Grids

Interest in underground CO2 storage has increased significantly over the last decade, driven by growing concern about global warming and rising levels of greenhouse gases in the atmosphere. Given that CO2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is considered one of the most direct approaches to achieving the net-zero carbon target. Although CO2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. This work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model CO2 flow in faulted reservoirs, we present the first transient embedded discrete-fracture model for corner-point grids (tEDFM-CPG). By comparing the tEDFM-CPG results with high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO2 flow in faulted/fractured reservoirs. In conclusion, this work presents the use of mixed reality (MR) to efficiently observe CO2 gas migration in the interior of these corner-point grid systems.

02 PETROLEUM↗
DOE OSTI2026

Modeling CO 2 flow through faulted/fractured reservoirs using tEDFM in corner-point grids

The interest in underground CO 2 storage has increased significantly over the last decade because of the rising concern about global warming due to the growing levels of greenhouse gases in the atmosphere. Considering that CO 2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is regarded as one of the most direct approaches to achieving the net zero carbon target. Although CO 2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. Here, this work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model the flow of CO 2 in faulted reservoirs, we present the first transient embedded discrete fracture model for corner-point grids (tEDFM-CPG). By comparing the results of the tEDFM-CPG to high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO 2 flow in faulted/fractured reservoirs. Finally, this work presents the use of mixed reality (MR) to efficiently observe CO 2 gas migration in the interior of these corner-point grid systems.

25 ENERGY STORAGE↗
DOE OSTI2020

High‐order cubature rules for tetrahedra

Summary In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2024

Shifts in belowground processes along a temperate forest edge

Abstract Context Forests are increasingly fragmented, and as a result most forests in the United States are within one km of an edge. Edges change environmental conditions of the forest—especially radiation, roughness, temperature, and moisture—that can have consequences for plant productivity and ecosystem functions. However, edge effects on aboveground characteristics of plants and the environment are better understood relative to plant roots and soil in the belowground environment. Objectives Our main objectives were to determine if soil C pools and fluxes are higher at the edge relative to other landscape positions, and to understand how specific belowground processes contribute to bulk differences in pools and fluxes. Methods We measured environmental conditions, live and dead fine root traits, soil chemistry, and soil respiration along a 75 m transect from interior forest to meadow in Gaithersburg, MD. Results We observed differences in the soil chemical, biological and hydrological environment between the forest interior, edge and adjacent meadow. In some cases, the forest edge represented a mid-point in environmental or belowground characteristics between the forest interior and meadow ( e.g. , pH, C-to-N ratio [C:N], live fine root biomass, heterotrophic respiration), likely reflecting the change in litter type and quality associated with the transition from grass to woody species. In other cases, neighboring landscape positions were different from the forest edge, which was drier and had higher dead fine root biomass. Although soil C contents were not significantly different across landscape positions, there was a tendency towards higher average soil C content at the edge relative to other landscape positions, suggesting that increased C loss related to root decay and greater soil respiration at the edge relative to the forest interior may have been offset by increased C gain from high plant productivity and subsequent inputs to soil. Conclusions This research provides insight into how forest edge environments may differ from the interior and how concurrent processes above- and belowground may contribute to those differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A “Floatilla” of Airborne Seismometers for Venus

Abstract Barometers floating on high‐altitude balloons in the relatively clement cloud layer on Venus could detect and characterize acoustic waves generated by seismic activity, avoiding the need for high‐temperature electronics required for surface seismology. Garcia et al. (2022, https://doi.org/10.1029/2022GL098844 ) recently demonstrated the detection of low‐frequency sound (infrasound) caused by earthquakes of magnitudes 7.3 and 7.5 from stratospheric balloons nearly 3,000 km away from the epicenter. They provided a preliminary demonstration of earthquake magnitude and location inversion, and the determination of S‐ and Rayleigh wave velocities using only their acoustic signature. Large earthquakes produce low‐frequency seismic waves that penetrate the interiors of planets; their detection at continental‐scale distances from a high‐vantage point demonstrates the feasibility of balloon‐based investigations of Venus' interior. We contextualize these results within the effort to perform seismology on Venus from balloons, discuss its limitations, and share perspectives on open research questions in this area.

Krishnamoorthy, Siddharth↗
DOE OSTI2020

H-Canyon exhaust tunnel LiDAR data request for information review and release

This request is to release LiDAR data files acquired by SRNL R&D Engineering and H-Area Operations and Engineering while performing scans of the interior of the H-Canyon Exhaust tunnel. LiDAR data includes 3-dimensional (3D) point cloud data sets and panoramic digital images of the interior tunnel wall surfaces. Presently two deployments have been completed, the first in November 2019, Fig. 1, and the second in June of 2020. This request is for the release of the data collected during those two deployments. Detailed information on the deployment and data collected can be found in SRS document C-ESR-H-00072, “November 2019 Initial Deployment of LiDAR”. It is planned to perform ongoing scans at approximately 6-month intervals, the purpose of the deployments is to evaluate the usefulness of the data collected to enable quantitative measurements such as tunnel dimensions and rate of surface erosion and as a precursor to a potential deployment of a LiDAR system on the tunnel inspection crawler.

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

Density-related properties from self-interaction corrected density functional theory calculations

The Perdew–Zunger self-interaction correction (PZ-SIC) removes unphysical electron self-interaction from calculations employing standard density functional approximations. Doing so improves many computed properties, bringing them into better agreement with experimental observations or with results from high-level quantum chemistry calculations. However, while PZ-SIC generally corrects in the right direction relative to the corresponding reference values, in many cases, it over-corrects. For this reason, scaled-down versions of PZ-SIC have been proposed and investigated. These approaches have mostly employed exterior scaling in which SIC correction terms are scaled in the same way at every point in space. Recently, a new local, or interior, scaling SIC method was proposed on non-empirical grounds to restore a property of the exact, but unknown, density functional that is broken in PZ-SIC. Here, the scaling at each point depends on the character of the charge density at that point. However, the local scaling can be done in various ways while still restoring the behavior of the exact functional. In this work, we compare and contrast the performance of various interior scaling approaches for addressing over-corrections of calculated molecular dipole moments and atomic polarizabilities and properties that reflect the nature of the electronic charge density.

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

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. Here, we used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

58 GEOSCIENCES↗