Search NASASearch

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 145 records · Page 8

From Points to Planes: A Workflow for Converting Three‐Dimensional Point Cloud Data Into Discrete Fracture Network Flow and Transport Models

We present the Point cLoud Algorithm for NEtwork Extraction of Discrete Fracture Networks (PLANE-DFN), a point cloud–based algorithm for automatic fracture network extraction designed to support discrete fracture network (DFN) modeling workflows. PLANE-DFN segments three-dimensional fracture planes from raw point cloud data using RANdom SAmple Consensus coupled with statistical outlier removal and density-based clustering to isolate individual fracture features. Each candidate plane is constrained against site-specific structural constraints based on strike and dip. After segmentation, each fracture is converted into a 2-D convex polygon suitable for meshing and simulation. The PLANE-DFN algorithm is validated by comparing geometric and flow and transport data against data from dfnWorks simulations with ensembles of plane-fit networks. We find that the flow and transport in plane-fit networks are comparable to dfnWorks-generated networks when realistic network geometry is maintained. The PLANE-DFN algorithm provides an automated and streamlined workflow to transform point clouds of data into DFN network geometry.

54 ENVIRONMENTAL SCIENCES

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)

Coefficient-to-Basis Network: a fine-tunable operator learning framework for inverse problems with adaptive discretizations and theoretical guarantees

We propose a Coefficient-to-Basis Network (C2BNet), a novel framework for solving inverse problems within the operator learning paradigm. C2BNet efficiently adapts to different discretizations through fine-tuning, using a pre-trained model to significantly reduce computational cost while maintaining high accuracy. Unlike traditional approaches that require retraining from scratch for new discretizations, our method enables seamless adaptation without sacrificing predictive performance. Furthermore, we establish theoretical approximation and generalization error bounds for C2BNet by exploiting low-dimensional structures in the underlying datasets. Our analysis demonstrates that C2BNet adapts to low-dimensional structures without relying on explicit encoding mechanisms, highlighting its robustness and efficiency. To validate our theoretical findings, we conducted extensive numerical experiments that showcase the superior performance of C2BNet on several inverse problems. The results confirm that C2BNet effectively balances computational efficiency and accuracy, making it a promising tool to solve inverse problems in scientific computing and engineering applications.

97 MATHEMATICS AND COMPUTING

Block encodings of discrete subgroups on a quantum computer

We introduce a block encoding method for mapping discrete subgroups to qubits on a quantum computer. This method is applicable to general discrete groups, including crystal-like subgroups such as BI of S U ( 2 ) and V of S U ( 3 ) . We detail the construction of primitive gates—the inversion gate, the group multiplication gate, the trace gate, and the group Fourier gate—utilizing this encoding method for BT and for the first time BI group. We also provide resource estimations to extract the gluon viscosity. The inversion gates for BT and BI are benchmarked on the quantum computer with estimated fidelities of 40 − 4 + 5 % and 4 − 3 + 5 % , respectively. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)

Laser wakefield acceleration driven by a discrete flying focus

Laser wakefield acceleration (LWFA) may enable the next generation of TeV-scale lepton colliders. Reaching such energies will likely require multiple LWFA stages to overcome limitations on the energy gain achievable in a single stage. The use of stages, however, introduces challenges such as alignment, adiabatic matching between stages, and a lower average accelerating gradient. Here, we propose a discrete flying focus that can deliver higher energy gain in a single stage, thereby reducing the number of stages required for a target energy. A sequence of laser pulses with staggered focal points and delays drives a plasma wave in which an electron beam experiences a near-constant accelerating gradient over distances beyond those attainable with a conventional pulse. Simulations demonstrate that a discrete flying focus with a total energy of 150 J can transfer 40 GeV per electron to a 50-pC beam in a single 30-cm stage, corresponding to 50 dephasing lengths.

laser wakefield acceleration

Thermal Management of Discretized Heaters Using CuW Microchannel Heat Sinks and FC3283

Single-phase cooling using microchannel heat sinks (MCHS) has become a popular approach for overcoming the thermal challenges associated with high-powered microelectronic devices. Thermal management is one of the largest barriers to higher power densities in electronics and frequently limits overall device performance. The implementation of forced convective cooling via single-phase liquid cooling in MCHS reduces the thermal resistance resulting in lower device temperatures at high-power conditions, which can decrease the package size and extend the lifespan of devices. The goal of the work described in this article was to investigate practical cooling solutions for laser diode bars. This study examined the effectiveness of a copper tungsten (CuW) microchannel heat sink paired with a dielectric coolant (FC3283) for dissipating both discrete and uniform heat fluxes up to 600 W/cm2 across a 0.25-cm2 surface area through a numerical and experimental study. CuW was chosen as the MCHS material because it is thermal expansion matched to GaAs, which is a common laser diode substrate. FC3283 serves as a dielectric coolant that is compatible with power electronics cooling. The microchannels utilized in this work were approximately $365$ $μ$ m deep and $100$ $μ$ m wide resulting in a hydraulic diameter of $160$ $μ$ m. The study investigated the cooling performance across a variety of applied power loads and flow rates to determine optimal effectiveness. The resulting thermal resistance ranged from 0.15 cm2 K/W at the highest flow rate to 0.26 cm2 K/W at the lowest flow rate. Further, the resulting thermal performance from this study demonstrates the importance of considering discrete heat sources separately from uniform heat sources and proved that the unique combination of CuW microchannels with FC3283 can be a promising cooling option toward future advancements of laser diode bars and other high-power microelectronics.

42 ENGINEERING

Generic Discretization Library

The GenDiL library is a collection of C++ software abstractions designed to discretize and solve partial differential equations (PDEs) for high-performance computing (HPC) applications. Its primary focus is on modern C++ generic programming, which helps ensure portability across various hardware architectures. The central idea behind the library is to provide building blocks for numerical algorithms-such as discretization methods and iteration patterns-so that domain experts can focus on the math, rather than the low-level details of hardware or implementation. By defining abstractions for data types, iteration over computational grids, and scheduling of operations, the library isolates the high-level PDE algorithms from the platform-specific optimizations needed to achieve efficient performance.

Dudouit, Yohann [Lawrence Livermore National Labor

Hybrid Quantum Networks with Discrete Polarizations and Continuous Quadrature Variables

Final Scientific/Technical Report for DOE Project DE-SC0022069. Quantum communication and computation systems have evolved along two largely independent paradigms: discrete-variable (DV) systems that encode information in qubits such as photon polarizations or photon-number states, and continuous-variable (CV) systems that encode information in optical field quadratures. Each approach offers distinct advantages—DVs provide low error rates and compatibility with single-photon platforms, while CVs support deterministic operations and efficient quantum state manipulation. A fundamental challenge for building a scalable Quantum Internet lies in interfacing these two regimes into a unified hybrid architecture that can coherently distribute and process quantum information across heterogeneous quantum nodes. This project aims to develop and demonstrate a hybrid optical quantum network that seamlessly integrates DV and CV systems. Specifically, we investigate a new class of hybrid entanglement between the discrete polarizations of single photons and the continuous quadrature variables of optical cat states, overcoming incompatibilities in existing DV–CV demonstrations. Using this new entanglement resource, the team investigates a multi-node hybrid quantum local area network (Q-LAN) testbed capable of DV–CV entanglement generation, swapping, and distribution across fiber links.

74 ATOMIC AND MOLECULAR PHYSICS

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD

Emergent discrete space-time crystal of Majorana-like quasiparticles in chiral liquid crystals

Time crystals spontaneously break the time translation symmetry, as recently has been frequently reported in quantum systems. Here we describe the observation of classical analogues of both 1+1-dimensional and 2+1-dimensional discrete space-time crystals in a liquid crystal system driven by a Floquet electrical signal. These classical time crystals comprise particle-like structural features and exists over a wide range of temperatures and electrical driving conditions. The phenomenon-enabling period-doubling effect comes from their topological Majorana-like quasiparticle features, where periodic inter-transformations of co-existing topological solitons and disclinations emerge in response to external stimuli and play pivotal roles. Our discrete space-time crystals exhibit robustness against temporal perturbations and spatial defects, behaving like a time-crystalline analogues of a smectic phase. Our findings show that the simultaneous symmetry breaking in time and space can be a widespread occurrence in numerous open systems, not only in quantum but also in a classical soft matter context.

FOS: Physical sciences

Discrete control of linear distributed systems with application to the deformable primary mirror of a large orbiting telescope

One of the more significant technological problems associated with the orbital operation of large astronomical telescope's is the fabrication and maintenance of the primary mirror surface to the tolerance required for diffraction-limited performance. An interesting approach to the solution of this problem involves continuously measuring and automatically correcting the optical surface of a thin deformable mirror by means of discrete actuators located on its rear surface: The realization of diffraction-limited performance from a telescope in space by this method rests on the ability of the designer to achieve extremely accurate control of a highly complex, interacting, multivariable system. This paper presents the results of a detailed study of the discrete control of linear distributed systems with specific application to the design of a practical controller for a plant representative of a telescope primary mirror for an orbiting astronomical observatory. The problem of controlling the distributed plant is treated by employing modal techniques to represent variations in the optical figure. Distortion of the mirror surface, which arises primarily from thermal gradients, is countered by actuators working against a backing structure to apply a corrective force distribution to the controlled surface. Each displacement actuator is in series with a spring attached to the mirror by means of a pad intentionally introduced to restrict the excitation of high-order modes. Control is then exerted over a finite number (equal.to the number of actuators) of the most significant modes. Through the application of the modal expansion technique the mirror equation of motion is transformed tb a set of uncoupled, linear, time-invariant, ordinary differential equations. The desired dynamic response and static accuracy may then be achieved by the application of classical single-variable design techniques. The formulation of a quadratic performance index which incorporates a measure of image quality permits determination of the trade-off between the-number of actuators and optical purity. A criterion for defining actuator placement and pad size is presented which minimizes the tendency of the controller to excite the unmonitored modes.

Jeremiah F Creedon

An Evaluation of Current NASTRAN Discrete Element Models for Monocoque and Semimonocoque Structures

The current discrete element analysis models in NASTRAN, Levels 8. 1.0 and 11. 1. 0, are evaluated for semi-rnonocoque and monocoque structures. Comparative results are presented with the FORMAT-I11 general-purpose code for a swept wing statics analysis and a Skylab substructure vibration analysis. Good agreement between Level i I. 1.0 and FORMAT-111 is demonstrated; however, a major error in Level 8. 1.0 was discovered on the Skylab vibration problem. The computational performance of NASTRAN was superior to FORMAT, but the vibration rigid format was found to have several inefficiencies. Comparative results are also presented between Level 11. 1. 0 and STACUSS, SLADE, BOND, UNIVALVE, and another special-purpose code. Results for two of the structures analyzed, a cylindrical shell and an elliptic plate, indicate curved discrete elements are significantly more accurate than the plate elements in NASTRAN for the same number of degrees of freedom.

Edward L. Stanton

Discrete ordinates-Monte Carlo coupling: A comparison of techniques in NERVA radiation analysis

In the radiation analysis of the NERVA nuclear rocket system, two-dimensional discrete ordinates calculations are sufficient to provide detail in the pressure vessel and reactor assembly. Other parts of the system, however, require three-dimensional Monte Carlo analyses. To use these two methods in a single analysis, a means of coupling was developed whereby the results of a discrete ordinates calculation can be used to produce source data for a Monte Carlo calculation. Several techniques for producing source detail were investigated. Results of calculations on the NERVA system are compared and limitations and advantages of the coupling techniques discussed.

Lindstrom, D. G.

Discrete sonic jets used as boundary-layer trips at Mach numbers of 6 and 8.5

The effect of discrete three-dimensional sonic jets used to promote transition on a sharp-leading-edge flat plate at Mach numbers of 6 and 8.5 and unit Reynolds numbers as high as 2.5 x 100,000 per cm in the Langley 20-inch hypersonic tunnels is discussed. An examination of the downstream flow-field distortions associated with the discrete jets for the Mach 8.5 flow was also conducted. Jet trips are found to produce lengths of turbulent flow comparable to those obtained for spherical-roughness-element trips while significantly reducing the downstream flow distortions. A Reynolds number based upon secondary jet penetration into a supersonic main flow is used to correlate jet-trip effectiveness just as a Reynolds number based upon roughness height is used to correlate spherical-trip effectiveness. Measured heat-transfer data are in agreement with the predictions.

Stone, D. R.

A design study for the addition of higher order parametric discrete elements to NASTRAN

The addition of discrete elements to NASTRAN poses significant interface problems with the level 15.1 assembly modules and geometry modules. Potential problems in designing new modules for higher-order parametric discrete elements are reviewed in both areas. An assembly procedure is suggested that separates grid point degrees of freedom on the basis of admissibility. New geometric input data are described that facilitate the definition of surfaces in parametric space.

Stanton, E. L.

Development of a container for handling, testing, and storing discrete microelectronic components

A container has been developed for handling, testing, burning-in, and storing discrete microelectronic components without removal from the protective package. The package was designed to accommodate the leadless inverted device and other carrier-mounted active devices and chip-type discrete resistors and capacitors. Before the indicated development, components were handled and tested in various ways, some of which resulted in damage or contamination. The basic design of the container utilizes precision machined printed circuit boards and chemically milled (photoetched) contact springs. Included in this design for protection is an O-ring-sealed cover. Methods of fabrication and limitations of the current hardware are presented. Current applications of and possible extensions to the technology are discussed.

Filip, G. L.