Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discretization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Insights into the soft brittle-to-ductile transition from discrete dislocation dynamics

The Brittle-to-ductile transition (BDT) in body centered cubic metals exhibits a soft transition wherein the fracture toughness gradually rises to before the onset of ductility. The resultant brittle-to-ductile transition temperature can be described with an Arrhenius relationship whose activation energy is related to plasticity in the material. To provide further insight into the nature of the BDTT, in this work we utilized a discrete dislo- cation dynamics model with a crack to simulate the BDT and how it depends on the thermally activated nature of plasticity. The interrelationship between the BDT activation energy and the dislocation mobility parameters were determined via the calculation of first order sensitivity coefficients. This analysis allows us to demonstrate that the activation energy for the BDT is directly related to the activation energy for plasticity through an effective stress that defines this relationship. This effective stress physically is the average stress on the dislocations that move out of the crack. Lastly, we are able to show that this effective stress is dictated by the low temperature fracture toughness or cleave energy of the material and the source position, the latter of which can be affected by processing. Collectively, these results provide new insight into what controls the thermal activation of the BDT and what are the important parameters to control it.

36 MATERIALS SCIENCE↗
DOE OSTI2024

All-fiber, discrete arbitrary picket-pulse shaping

Here, a novel all-fiber, hybrid pulse-shaping system based on the temporal stacking of discrete optical picket pulses combined with continuously shaped optical pulses provides 100-ps resolution with the high dynamic range and flexibility required for laser fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Discrete element modeling of irregular-shaped soft pine particle flow in an FT4 powder rheometer

Pine residues are a commonly used biomass feedstock that consists of different anatomical fractions, each with distinct particle characteristics. Here, in this work, an experiment-informed discrete element model is developed to investigate the flowability of pine residues in an FT4 rheometer. Multi-sphere particles with distinct particle attributes are created to model each anatomical fraction type. A systematic analysis of specimens with varying particle characteristics (e.g., anatomical fraction type, particle shape, and size) is conducted to elucidate the relationship between particle attributes and flowability. The results show that stems recorded the highest axial force and torque and, correspondingly, the highest flow energy, which is attributed to their high stiffness and interlocking effect. Increasing their percentage in the mixture increases the flow energy while increasing needles tends to decrease flow energy. Knowledge gained in this study on the flow of anatomical fractions is important for the efficient and robust processing of pine residues.

09 BIOMASS FUELS↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2026

A Gauss-Radau-Laguerre Discrete Variable Representation for Use in Continuum Electron Dynamics

In this work, we detail an implementation, suitable for calculations on highly correlated ionizing systems, of a modified finite element discrete variable representation (FE-DVR) appended with a Gauss-Radau-Laguerre element. The appended element includes exterior complex scaling (ECS) to impose outgoing wave boundary conditions on treatments of processes involving continuum electrons. In this “infinite range” ECS (irECS), the complications that introduce reflections from the end of the grid when the last ECS finite element has finite range are avoided by the use of the Laguerre-weighted exponentially decaying tails, while outgoing wave boundary conditions are still imposed via the ECS transformation. For highly correlated systems in the absence of strong external fields we find that accurate two-electron integrals are essential in this modified FE-DVR. To accurately compute the two-electron integrals over the entire ECS contour, we present a detailed examination of the implications from the boundary terms that arise in a solution of Poisson’s equation with the Radau-Laguerre basis. A boundary term correction is necessary, and when included, the Radau-Laguerre DVR can accurately describe highly correlated states such as the doubly excited states of helium over the entire ECS contour.

elements↗
DOE OSTI2026

Alkali-Metal Interlocking of 2D V 4 O 10 Sheets Defines Discretized Interlayer Shear Relationships

Low-dimensional materials manifest structural anisotropy, quantum confinement, and tightly bound excitonic states, which make them attractive building blocks that can be assembled within three-dimensional laterally stitched heterostructures, stacked van der Waals solids, and complex moiré superlattices. Ion intercalation in the galleries between layered materials provides a means of modifying interlayer separation and coupling, but it is also known to drive the shearing of the layers. In this article, we explore the distinct ligand coordination environments afforded by vanadyl oxygens of singular [V 4 O 10 ] sheets and examine how the size, polarizability, and stoichiometry of Group I cations sandwiched between such layers determine the interlocking of the sheets in stacked structures. Based on the topochemical insertion of alkali-metal ions into the layered λ-V 2 O 5 , we identify seven types of guest ion coordination sites discretized into four distinct regimes of interlayer shear in units of half octahedral widths. The coordination preferences of intercalated cations govern how they interlock 2D [V 4 O 10 ] sheets and engender specific shear conformations. We present evidence that static and dynamic disorder in guest ion arrangement modulate the magnetic structure of the intercalated compounds based on electrostatic polarization, localization of charge and spin density, and lattice distortion. The results illustrate the use of topochemical ion insertion to modulate stacking relationships and magnetic transition characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Energy Transport and Conversion Above a Bright Discrete Auroral Arc

Magnetically connected observations of particle distributions and luminosity from the Reimei spacecraft are used to examine energy transport and conversion occurring above a discrete auroral arc. By combining imaging and in situ measurements it is shown how transverse electromagnetic and kinetic energy fluxes measured along the spacecraft trajectory converge across geomagnetic field-lines into the acceleration region. It is shown how cross-field energy transport is facilitated by the formation of vortices along the length of the arc. From an integration over the vertical extent of the acceleration region it is shown that the transverse and field-aligned flow of energy into the arc locally supports the dissipation needed to power the electron acceleration observed. Estimates of gradients in electromagnetic and kinetic energy flows show how field-aligned electron energization in the acceleration region is supported by the divergence of Poynting flux along and across the background magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A discretized representation for Monte Carlo simulation of deformed semiflexible chains

In this study, we present a novel orientation discretization approach based on the rhombic triacontahedron for Monte Carlo simulations of semiflexible polymer chains, aiming at enhancing structural analysis through rheo-small-angle scattering (rheo-SAS). Our approach provides a more accurate representation of the geometric features of semiflexible chains under deformation, surpassing the capabilities of traditional lattice structures. Validation against the Kratky–Porod chain system demonstrated superior consistency, underscoring its potential to significantly improve the precision of uncovering geometric details from rheo-SAS data. Further, this approach opens new avenues for investigating the conformations of semiflexible polymers and mechanically induced phase transitions in more complex polymer structures, offering deeper insights into their behavior under various conditions.

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

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗
DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗
DOE OSTI2025

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anisotropic excitonic magnetism from discrete C 4 symmetry in CeRhIn 5

Anisotropy in strongly correlated materials is a central parameter in determining the electronic ground state and is tuned through the local crystalline electric field. This is notably the case in the CeCo x Rh 1 − x In 5 system where the ground-state wave function can provide the basis for antiferromagnetism and/or unconventional superconductivity. We develop a methodology to understand the local magnetic anisotropy and experimentally investigate with neutron spectroscopy applied to antiferromagnetic ( T N = 3.8 K ) CeRhIn 5 , which is isostructural to d -wave superconducting ( T c = 2.3 K ) CeCoIn 5 . Through diagonalizing the local crystal field Hamiltonian with discrete tetragonal C 4 point group symmetry and coupling these states with the random phase approximation, we find two distinct modes polarized along the crystallographic c and a − b planes, agreeing with experiment. The anisotropy and bandwidth, underlying the energy scale of these modes, are tuneable with a magnetic field which we use experimentally to separate in energy single and multiparticle excitations thereby demonstrating the instability of excitations polarized within the crystallographic a − b plane in CeRhIn 5 . We compare this approach to a S eff = 1 2 parametrizations and argue for the need to extend conventional SU(2) theories of magnetic excitations to utilize the multilevel nature of the underlying crystal-field basis states constrained by the local point-group C 4 symmetry. Published by the American Physical Society 2024

Brener, D. J. (ORCID:0000000333400177)↗
DOE OSTI2025

Hidden conformal symmetry of the discrete series scalars in dS 2

In D dimensional de Sitter space, a scalar field has an infinite tower of special tachyonic mass values at which enhanced shift symmetries appear. After modding out by these shift symmetries, these fields correspond to the unitary irreducible representations of the de Sitter group known as the discrete series. We show that in D = 2 these theories have global conformal symmetry. In all but the massless case, these theories have no stress tensor and the conformal symmetry does not act in the usual way on the scalar field. We find the conformal symmetry by explicitly computing the correlators of the shift invariant local operators and showing that they take conformally invariant forms. We also demonstrate how these fields are self-dual in D = 2 , and dual to the shift invariant massive vector fields, which are therefore also conformally invariant. Published by the American Physical Society 2025

Farnsworth, Kara (ORCID:000000020200078X)↗
DOE OSTI2024

Giant Graviton Expansion from Bubbling Geometry: Discreteness from Quantized Geometry

The superconformal index of half-BPS states in N = 4 supersymmetric Yang-Mills with gauge group U⁡(N) admits an expansion in terms of giant gravitons, J N (q) = J ∞ ⁡(q)⁢Σ$^{∞}_{m=0}$ q m⁢N ⁢ J^ m ⁡(q), where m is the number of giant gravitons and J ∞ ⁡(q) is the graviton index. The expansion can be viewed as the implementation of trace relations for finite N. We derive this expansion directly in supergravity from the class of half-BPS solutions due to Lin, Lunin, and Maldacena in type IIB supergravity. The moduli space of these configurations can be quantized using covariant quantization methods. We show how this quantization leads to the precise expression for the expansion in terms of giant gravitons. Our proposal provides a derivation of the giant graviton expansion directly in terms of quantized supergravity degrees of freedom, and it recovers discrete data via quantum geometries that are classically nonsmooth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Fault-Tolerant Operation of Bosonic Qubits with Discrete-Variable Ancillae

Fault-tolerant quantum computation with bosonic qubits often necessitates the use of noisy discrete-variable ancillae. In this work, we establish a comprehensive and practical fault-tolerance framework for such a hybrid system and synthesize it with fault-tolerant protocols by combining bosonic quantum error correction (QEC) and advanced quantum control techniques. We introduce essential building blocks of error-corrected gadgets by leveraging ancilla-assisted bosonic operations using a generalized variant of path-independent quantum control. Using these building blocks, we construct a universal set of error-corrected gadgets that tolerate a single-photon loss and an arbitrary ancilla fault for four-legged cat qubits. Notably, our construction requires only dispersive coupling between bosonic modes and ancillae, as well as beam-splitter coupling between bosonic modes, both of which have been experimentally demonstrated with strong strengths and high accuracy. Moreover, each error-corrected bosonic qubit is comprised of only a single bosonic mode and a three-level ancilla, featuring the hardware efficiency of bosonic QEC in the full fault-tolerant setting. We numerically demonstrate the feasibility of our schemes using current experimental parameters in the circuit-QED platform. Finally, we present a hardware-efficient architecture for fault-tolerant quantum computing by concatenating the four-legged cat qubits with an outer qubit code utilizing only beam-splitter couplings. Our estimates suggest that the overall noise threshold can be reached using existing hardware. These developed fault-tolerant schemes extend beyond their applicability to four-legged cat qubits and can be adapted for other rotation-symmetrical codes, offering a promising avenue toward scalable and robust quantum computation with bosonic qubits. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Endogenous Interface Pricing for Consistent Transmission–Distribution Co-Optimization With Discrete Distribution Controls

This paper proposes an endogenous interface pricing model for day-ahead transmission–distribution co-optimization that co-determines the interface locational marginal price (LMP) and the transmission–distribution exchange, ensuring price–dispatch consistency while optimally scheduling discrete distribution controls. The formulation couples a DC optimal power flow (OPF) with a branch-flow AC OPF that schedules distributed energy resources (DERs), tap-changer settings, capacitor banks (CBs), and multi-period energy storage systems (ESSs) under feeder voltage and current limits, and is solved as a mixed-integer second-order cone program (MISOCP). In a T14–D33 system, coordinated device scheduling recovers about 90% of the distribution-to-transmission export achievable in a reference case that ignores distribution network (DN) limits, while satisfying a 1.05 p.u. voltage upper bound. In a T39–D34/D37/D123 system, a sequential decoupled benchmark produces interface LMP distortions up to 12.5% and a 7.28% mismatch in net export energy, whereas the proposed model removes these distortions and the associated settlement mismatches. Second-order cone (SOC) relaxation gaps remain below $10^{-3}$ in all cases.

Noh, Seung-Gil↗
DOE OSTI2025

Discrete-Event Model of WIPP Operations

The Waste Isolation Pilot Plant (WIPP) is the critical component of the Department of Energy's (DOE) Transuranic Radioactive Waste (TRU) disposition infrastructure. Quantifying the operations ofWIPP with a discrete-event model demonstrates the capability to assess efficiency, identity bottlenecks, and improve future operations. Such a model has been developed with the simulation software ExtendSim. This report outlines the structure of that model, summarizes the model's successful reproduction of the annual amount of emplaced waste, and demonstrates that WIPP is successfully receiving and emplacing waste at a rate consistent with the rate at which the waste arrives. The model serves as a first step toward future production enhancements at WIPP. Those enhancements will rely on close collaboration with the Carlsbad Field Office (CBFO), accurate interpretation and incorporation of the model results, and effective planning with other DOE Environmental Management (DOE-EM) entities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗