Search NASA⌕ Search

SEARCH · Search NASA

Results for “local solution bounds”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimization-based, property-preserving finite element methods for scalar advection equations and their connection to Algebraic Flux Correction

In this paper, we continue our efforts to exploit optimization and control ideas as a common foundation for the development of property-preserving numerical methods. Here we focus on a class of scalar advection equations whose solutions have fixed mass in a given Eulerian region and constant bounds in any Lagrangian volume. Our approach separates discretization of the equations from the preservation of their solution properties by treating the latter as optimization constraints. This relieves the discretization process from having to comply with additional restrictions and makes stability and accuracy the sole considerations in its design. A property-preserving solution is then sought as a state that minimizes the distance to an optimally accurate but not property-preserving target solution computed by the scheme, subject to constraints enforcing discrete proxies of the desired properties. Furthermore, we consider two such formulations in which the optimization variables are given by the nodal solution values and suitably defined nodal fluxes, respectively. A key result of the paper reveals that a standard Algebraic Flux Correction (AFC) scheme is a modified version of the second formulation obtained by shrinking its feasible set to a hypercube. In conclusion, we present numerical studies illustrating the optimization-based formulations and comparing them with AFC

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On the continuity of solutions to anisotropic elliptic operators in the limiting case

Abstract We show that local weak solutions to anisotropic elliptic equations with bounded and measurable coefficients, whose prototype is are continuous. Our analysis exploits a particular intrinsic geometry that, within the limiting case on the exponents, allows for a reduction of the oscillation of the solution.

Ciani, Simone [Department of Mathematics, Universi↗
DOE OSTI2020

Kink band and shear band localization in anisotropic perfectly plastic solids

Shear-driven strain localization has been observed in a wide variety of materials and may take the form of shear bands or kink bands. Based on observations of kink bands in plastically anisotropic metallic nanolaminates and single crystal metals, here we posit that, for the specific case of isochoric deformation, the kinematics of kink band formation are indistinguishable from those of plane strain shear band formation. The only distinction between shear bands and kink bands in these systems would then be that kink bands ‘lock up’ at a particular value of material rotation while shear bands may progress to arbitrarily high strains. In order to investigate whether strong material anisotropy is sufficient to arrest shear localization at a geometry that matches the classic kink band geometry, we model the development of a band of simple shear within an anisotropic perfectly plastic material. The resulting analytical model provides the stress state needed to maintain the kinematics of simple shear as a function of material anisotropy, deformation band orientation, and shear strain (or equivalently, material rotation). It is found that plastic anisotropy can promote either kink band or shear band formation depending on the loading orientation. When the deviatoric stress is positive parallel to the plane of anisotropy, shear localization may progress without bound and a shear band is produced. When the deviatoric stress is negative parallel to the plane of anisotropy, shear localization is arrested after a certain material rotation, resulting in a kink band. Examination of the requisite applied stress state during kink band formation provides an explanation for the experimentally-observed ‘lock up’ geometry. Solutions for the band boundary inclination angle are obtained and used to provide bounds on permissible band angles for both shear bands and kink bands.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Certifiably Correct Range-Aided SLAM

We present the first algorithm capable of efficiently computing certifiably optimal solutions to range-aided simultaneous localization and mapping (RA-SLAM) problems. Robotic navigation systems are increasingly incorporating point-to-point ranging sensors, leading state estimation which takes the form of RA-SLAM. However, the RA-SLAM problem is more difficult to solve than traditional pose-graph SLAM; ranging sensor models introduce additional non-convexity, unlike pose-pose or pose-landmark measurements, a single range measurement does not uniquely determine the relative transform between the involved sensors, and RA-SLAM inference is highly sensitive to initial estimates. Our approach relaxes the RA-SLAM problem to a semidefinite program (SDP), which we show how to solve efficiently using the Riemannian staircase methodology. The solution of this SDP provides a high-quality initialization for our original RA-SLAM problem, which is subsequently refined via local optimization, as well as a lower-bound on the RA-SLAM problem's optimal value. Our algorithm, named certifiably correct RA-SLAM (CORA), applies to problems comprised of arbitrary pose-pose, pose-landmark, and ranging measurements. Evaluation on simulated and real-world marine examples shows that our algorithm frequently produces certifiably optimal RA-SLAM solutions; moreover, even suboptimal estimates are typically within 1-2\% of the optimal value.

Papalia, Alan↗
DOE OSTI2022

Geometrical picture of the electron–electron correlation at the large- D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron–electron correlations and might give a bound on the correlation energy. Furthermore, the results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron–electron correlations.

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

PERKS: a Locality-Optimized Execution Model for Iterative Memory-bound GPU Applications

Iterative memory-bound solvers commonly occur in HPC codes. Typical GPU implementations have a loop on the host side that invokes the GPU kernel as much as time/algorithm steps there are. The termination of each kernel implicitly acts the barrier required after advancing the solution every time step. We propose an execution model for running memory-bound iterative GPU kernels: PERsistent KernelS (PERKS). In this model, the time loop is moved inside persistent kernel, and device-wide barriers are used for synchronization. We then reduce the traffic to device memory by caching subset of the output in each time step in the unused registers and shared memory. PERKS can be generalized to any iterative solver: they largely independent of the solver's implementation. We explain the design principle of PERKS and demonstrate effectiveness of PERKS for a wide range of iterative 2D/3D stencil benchmarks (geomean speedup of 2.12x for 2D stencils and 1.24x for 3D stencils over state-of-art libraries), and a Krylov subspace conjugate gradient solver (geomean speedup of 4.86x in smaller SpMV datasets from SuiteSparse and 1.43x in larger SpMV datasets over a state-of-art library). All PERKS-based implementations available at: https://github.com/neozhang307/PERKS.

Zhang, Lingqi↗
DOE OSTI2023

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. Furthermore, this is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Preparing Distribution Utilities for Utility-Scale Storage and Electric Vehicles: A Novel Analytical Framework

Emerging distributed energy resources (DERs)—such as solar photovoltaics (PV), battery energy storage systems (BESS), and electric vehicles (EVs)—are expected to increase substantially in India in the coming years following policy-driven targets of the Government of India to modernize its electricity system, reduce greenhouse gas emissions (GHGs), and improve air quality. These emerging technologies can pose challenges to distribution utilities, forcing overhauls in planning and operational practices. They can also create challenges in power system infrastructure planning and cause more frequent system operational violations (e.g., network voltage bounds and loading thresholds) if not properly integrated. The impacts on the localized power distribution grid from these emerging technologies manifest in increased infrastructure investments and erratic shifts in demand patterns. These impacts are not yet well understood, and analytic solutions are not readily available. To address these challenges, the National Renewable Energy Laboratory (NREL), in collaboration with BSES Rajdhani Power Ltd. (BRPL), developed an advanced power distribution system impact analysis framework of BRPL’s distribution system. This framework helps analyze the readiness of the power distribution network to accommodate emerging technologies and the potential opportunities they might introduce. The framework has been predominantly set up to evaluate distributed PV, BESS, and EVs. In this collaboration between NREL and BRPL, we developed and evaluated the framework on two distribution feeders in the BRPL territory for various scenarios of BESS and EVs. BESS are evaluated for their effectiveness on the grid to mitigate present and future feeder overloading scenarios, and they are subsequently analyzed for their costs compared to the costs of traditional upgradation measures. Scenarios include assessing the effects of EV density on grid infrastructure upgrades and interlinking EV management with BESS integration.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Grey Galaxies’ as an endpoint of the Kerr-AdS superradiant instability

Kerr-AdS d+1 black holes for d ≥ 3 suffer from classical superradiant instabilities over a range of masses above extremality. We conjecture that these instabilities settle down into Grey Galaxies (GGs) — a new class of coarse-grained solutions to Einstein’s equations which we construct in d = 3. Grey Galaxies are made up of a black hole with critical angular velocity ω = 1 in the ‘centre’ of AdS, surrounded by a large flat disk of thermal bulk gas that revolves around the centre of AdS at the speed of light. The gas carries a finite fraction of the total energy, as its parametrically low energy density and large radius are inversely related. GGs exist at masses that extend all the way down to the unitarity bound. Their thermodynamics is that of a weakly interacting mix of Kerr-AdS black holes and the bulk gas. Their boundary stress tensor is the sum of a smooth ‘black hole’ contribution and a peaked gas contribution that is delta function localized around the equator of the boundary sphere in the large N limit. We also construct another class of solutions with the same charges; ‘Revolving Black Holes (RBHs)’. RBHs are macroscopically charged SO(d, 2) descendants of AdS-Kerr solutions, and consist of ω = 1 black holes revolving around the centre of AdS at a fixed radial location but in a quantum wave function in the angular directions. RBH solutions are marginally entropically subdominant to GG solutions and do not constitute the endpoint of the superradiant instability. Nonetheless, we argue that supersymmetric versions of these solutions have interesting implications for the spectrum of supersymmetric states in, e.g. N = 4 Yang-Mills theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Robustness of vortex-bound Majorana zero modes against correlated disorder

Here, we investigate the effect of correlated disorder on Majorana zero modes (MZMs) bound to magnetic vortices in two-dimensional topological superconductors. By starting from a lattice model of interacting fermions with a p x ±ip y superconducting ground state in the disorder-free limit, we use perturbation theory to describe the enhancement of the Majorana localization length at weak disorder and a self-consistent numerical solution to understand the breakdown of the MZMs at strong disorder. We find that correlated disorder has a much stronger effect on the MZMs than uncorrelated disorder and that it is most detrimental if the disorder correlation length ℓ is on the same order as the superconducting coherence length ξ. In contrast, MZMs can survive stronger disorder for ℓ $\ll$ξ as random variations cancel each other within the length scale of ξ, while an MZM may survive up to very strong disorder for ℓ$\gg$ξ if it is located in a favorable domain of the given disorder realization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient Intermolecular Energy Exchange and Soft-ionization of Water at Nanoplatelet Interfaces

X-ray, energetic photon and electron irradiation can ionize and electronically excite target atoms and molecules. These excitations undergo complicated relaxation and energy transfer processes that ultimately determine the manifold of system responses to the deposited excess energy. In weakly bound gas- and solution-phase samples, intermolecular Coulomb decay (ICD) and electron-transfer mediated decay (ETMD) can occur with neighboring atoms or molecules leading to efficient transfer of the excess energy to the surroundings. In ionic solids such as metal-oxides, intra- and inter-atomic Auger decay produces localized final states that lead to lattice damage and typically removal of cations from the substrate. The relative importance of Auger stimulated damage (ASD) versus ICD and ETMD in micro-solvated nanoparticle interfaces is not known. Though ASD is generally expected, essentially no lattice damage and very efficient ICD and ETMD resulting from ionization and electronic excitation of micro-solvated boehmite (AlOOH) nanoplatelets has been observed. Efficient energy transfer and soft ionization of interfacial molecules is likely a general phenomenon at gas-oxyhydroxide nanoparticle interfaces where the density of states of the ionized chemisorbed species significantly overlap the core hole states of the solid.

Jones, Brant M.↗
DOE OSTI2023

Modeling Nanoscale Ohmics in Carbon Supports of Fuel Cell Cathodes

Here, reducing platinum (Pt) loading in polymer electrolyte fuel cells (PEFCs) while meeting performace requirements is critical to their widespread deployment. However, significant polarization losses manifest at higher current densities in cathodes with lower Pt content. The morphology of the carbon supports in PEFC cathodes affects the location of Pt deposition into the micropits or onto the surface of the carbon support, translating into different kinetic and transport resistances. In this work, we derive an agglomerate scale model that differentiates the sink terms for Pt on the surface and in the pits of carbon supports. We develop an approach to assess the impact of nanoscale ohmic resistance to Pt in the micropits arising from weakly ionic solution in the carbon support on PEFC performance. Effectiveness factors relating the actual reaction rate to the maximum reaction rate (had all the Pt been exposed) are derived and embedded into a one-dimensional catalyst layer model. Parameters in the catalyst layer model are tuned based on experimental local oxygen transport resistances. Subsequently, we estimate bounds for the micropore resistances based on geometric and physical arguments. Lastly, polarization curves are simulated to assess the effect of the micropore resistance in fully-humidified and oxygen-rich environments.

25 ENERGY STORAGE↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Enhanced biochemical sensing with high- Q transmission resonances in free-standing membrane metasurfaces

Optical metasurfaces provide solutions to label-free biochemical sensing by localizing light resonantly beyond the diffraction limit, thereby selectively enhancing light–matter interactions for improved analytical performance. However, high-Q resonances in metasurfaces are usually achieved in the reflection mode, which impedes metasurface integration into compact imaging systems. Here, we demonstrate a metasurface platform for advanced biochemical sensing based on the physics of the bound states in the continuum (BIC) and electromagnetically induced transparency (EIT) modes, which arise when two interfering resonances from a periodic pattern of tilted elliptic holes overlap both spectrally and spatially, creating a narrow transparency window in the mid-infrared spectrum. We experimentally measure these resonant peaks observed in transmission mode (Q ~ 734 at λ ~ 8.8 µm) in free-standing silicon membranes and confirm their tunability through geometric scaling. We also demonstrate the strong coupling of the BIC-EIT modes with a thinly coated PMMA film on the metasurface, characterized by a large Rabi splitting (32 cm -1 ) and biosensing of protein monolayers in transmission mode. Our new photonic platform can facilitate the integration of metasurface biochemical sensors into compact and monolithic optical systems while being compatible with scalable manufacturing, thereby clearing the way for on-site biochemical sensing in everyday applications.

Rosas, Samir [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2023

Reduced dynamic complexity allows structure elucidation of an excited state of KRAS G13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS G13D , which suggests an exchange between two conformational states in solution. Methyl and 31 P NMR spectra of active KRAS G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major Pγ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS G13D and KRAS G13D -RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS G13D , showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

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

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗