Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite precision”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Temporal entanglement transitions in the periodically driven Ising chain

Periodically driven quantum systems can host non-equilibrium phenomena without static analogs, including in their entanglement dynamics. Here, we discover (TET) in a Floquet spin chain, which correspond to a quantum phase transition in the spectrum of the entanglement Hamiltonian and are signaled by dynamical spontaneous symmetry breaking. We identify the symmetry principles underlying these transitions: they appear when the driven Hamiltonian preserves global symmetry (here, ℤ2), the initial state respects this symmetry, and the reduced density matrix carries weight in both subsystem-parity sectors, with TET occurring precisely when the sector weights become equal (given the previous two conditions are also satisfied). Intriguingly, we find these transitions across a broad range of driving frequencies (from adiabatic to high-frequency regime) and independently of drive details, where they manifest as periodic, sharp entanglement spectrum reorganizations marked by the Schmidt-gap closure, a vanishing entanglement echo, and symmetry-quantum-number flips, while remaining invisible to conventional local observables. At high frequencies, the entanglement Hamiltonian acquires an intrinsic timescale decoupled from the drive period, rendering the transitions genuine steady-state features. Finite-size scaling reveals universal critical behavior with correlation-length exponent 𝜈 = 1, matching equilibrium Ising universality despite its emergence from purely dynamical mechanisms decoupled from static criticality. Furthermore, our work establishes TET as novel features in Floquet quantum matter.

36 MATERIALS SCIENCE↗
NASA NTRS2020

Broadband Characterization of Anti-reflection Coated Starshade Optical Edges for Solar Glint Control

A starshade is a large flower-shaped screen designed to enable the direct imaging of exoplanets with a space telescope. The starshade perimeter is composed of sharp, precisely shaped edges to minimize the glint of sunlight into the telescope. Past work has focused on bare edges to minimize the terminal radius. This paper describes the broadband, wide-angle performance of edges coated with a thin multi-layer anti-reflection coating. This coating uses a combination of interference and absorption to reduce the surface reflectivity and to avoid the negative effects associated with a large cross-sectional area. A custom scattered light testbed has been developed to quantify the amount of light scattered from sample edges and to validate Finite-Difference Time-Domain (FDTD) models of the optical scatter. We show that optical edge samples with this coating significantly reduce the solar glint pattern compared to similar uncoated optical edges.

Sheikh, David↗
DOE OSTI2026

Concatenated dual displacement code for continuous-variable quantum error correction

The continuous-variable (CV) Gaussian no-go theorem fundamentally limits the suppression of Gaussian displacement errors using only Gaussian gates and states. Prior studies have employed Gottesman-Kitaev-Preskill (GKP) states as ancillary qumodes to suppress small Gaussian displacement errors. However, when the displacement magnitude becomes large, inevitable lattice-crossing errors arise beyond the correctable range of the GKP state. To address this issue, we concatenate the Gaussian-noise-suppression circuit with an outer analog Steane code that corrects such occasional lattice-crossing events as well as other abrupt displacement errors. Contrary to conventional concatenation, which primarily aims to reduce logical error rates, the Steane-GKP duality in encoding provides complementary protection against displacement errors at different scales: The inner GKP layer employs non-Gaussian resources to suppress continuous Gaussian noise and reduce residual variance, while the outer analog Steane code corrects discrete lattice-crossing events that exceed the GKP correctable range. It is precisely this separation of error-mitigation roles that enables CV error correction. In contrast to prior work on concatenating GKP and repetition codes to establish error correction for discrete qubit/qudit encoding, we provide correction in the continuous encoding space. Analytical studies show that, under infinite squeezing, the concatenated code suppresses the variance of Gaussian displacement errors acting on all qumodes by up to 50%, while enabling unbiased correction of lattice-crossing errors with a success probability determined by the ratio between the residual Gaussian error standard deviation and the lattice-crossing magnitude. Even with finite squeezing, the proposed architecture still provides Gaussian-error suppression and lattice-crossing correction. Moreover, the presence of the outer analog Steane code relaxes the squeezing requirement of the inner GKP states, indicating near-term experimental feasibility. This work establishes a viable route toward fault-tolerant continuous-variable quantum computation and provides insight into the design of concatenated CV error-correcting architectures.

quantum error correction↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Control/structure interaction conceptual design tool

The JPL Control/Structure Interaction Program is developing new analytical methods for designing micro-precision spacecraft with controlled structures. One of these, the Conceptual Design Tool, will illustrate innovative new approaches to the integration of multi-disciplinary analysis and design methods. The tool will be used to demonstrate homogeneity of presentation, uniform data representation across analytical methods, and integrated systems modeling. The tool differs from current 'integrated systems' that support design teams most notably in its support for the new CSI multi-disciplinary engineer. The design tool will utilize a three dimensional solid model of the spacecraft under design as the central data organization metaphor. Various analytical methods, such as finite element structural analysis, control system analysis, and mechanical configuration layout, will store and retrieve data from a hierarchical, object oriented data structure that supports assemblies of components with associated data and algorithms. In addition to managing numerical model data, the tool will assist the designer in organizing, stating, and tracking system requirements.

Briggs, Hugh C.↗
DOE OSTI2021

Composites and additive manufacturing for high-field coil supports for stellarators

The utilization of composites layered on an additive manufactured substrate, for the production of coil supports for modular coils in small or middle size experimental stellarators is assessed. The focus of the study is a monolithic coil support comprising the coils of a half-period of a stellarator, somewhat similar to the ones in UST_2 and ARIES-CS stellarators. However, the concept may be applicable to quasi-monolithic coil supports (coil forms of NCSX type) or individual coil casings (W7-X type).Coil supports for stellarators require high precision, stiffness and strength for large contorted parts. Traditionally, monolithic coil supports are produced by steel casting/forging and final machining. This production method and material gives accurate, stiff and strong coil supports, but the method may be expensive for monolithic supports due to the geometrical complexity and required accuracy of the structure. In relation to those matters, this work investigates whether a monolithic coil support comprising an additive manufactured resin substrate, which is externally (outward from the coils) surrounded with a thick layer of fibre-reinforced resin, may achieve enough strength and stiffness under middle/high magnetic fields. Finite element calculations are produced to obtain the direction of the principal stresses and their values in compression and tension at different areas of the monolithic support, which is relevant for anisotropic materials. The feasibility of directional application of (carbon) fibres on the winding surface of the stellarator outward from the coils is experimentally tested on a scaled-down additively manufactured prototype of a monolithic support. The strength and stiffness of the composite structure appears sufficient for common magnetic fields in experimental stellarators, and the 3D-composite design and manufacturing was technically feasible.

Queral, Vicente↗
DOE OSTI2023

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2026

Signs of nonmonotonic finite-volume corrections to 𝑔 𝐴

We study finite-volume (FV) corrections to determinations of 𝑔 𝐴 via lattice quantum chromodynamics (QCD) using analytic results and numerical analysis. We observe that 𝑆⁢𝑈⁡(2) heavy Baryon chiral perturbation theory does not provide an unambiguous prediction for the sign of the FV correction, which is not surprising when one also considers large-𝑁 𝑐 constraints on the axial couplings. We further show that nonmonotonic FV corrections are naturally allowed when one considers either including explicit Δ-resonance degrees of freedom or one works to higher orders in the chiral expansion. We investigate the potential impact of these FV corrections with a precision study of 𝑔 𝐴 using models of FV corrections that are monotonic and nonmonotonic. Using lattice QCD data that is approximately at the 1% level of precision, we do not see significant evidence of nonmonotonic corrections. Looking forward to the next phase of lattice QCD calculations, we estimate that calculations that are between the 0.1% and 1% level of precision may be sensitive to these FV artifacts. Finally, we present an update of the CalLat prediction of 𝑔 𝐴 in the isospin limit with subpercent precision, 𝑔$^{QCD}_{𝐴}$ = 1.2674⁢(96).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Asymptotic domination of cold relativistic MHD winds by kinetic energy flux

We study the conditions which lead to the conversion of most Poynting flux into kinetic energy flux in cold, relativistic hydromagnetic winds. It is shown that plasma acceleration along a precisely radial flow is extremely inefficient due to the near cancellation of the toroidal magnetic pressure and tension forces. However, if the flux tubes in a flow diverge even slightly faster than radially, the fast magnetosonic point moves inward from infinity to a few times the light cylinder radius. Once the flow becomes supermagnetosonic, further divergence of the flux tubes beyond the fast point can accelerate the flow via the 'magnetic nozzle' effect, thereby further converting Poynting flux to kinetic energy flux. We show that the Grad-Shafranov equation admits a generic family of kinetic energy-dominated asymptotic wind solutions with finite total magnetic flux. The Poynting flux in these solutions vanishes logarithmically with distance. The way in which the flux surfaces are nested within the flow depends only on the ratio of angular velocity to poliodal 4-velocity as a function of magnetic flux. Radial variations in flow structure can be expressed in terms of a pressure boundary condition on the outermost flux surface, provided that no external toriodal field surrounds the flow. For a special case, we show explicitly how the flux surfaces merge gradually to their asymptotes. For flows confined by an external medium of pressure decreasing to zero at infinity we show that, depending on how fast the ambient pressure declines, the final flow state could be either a collimated jet or a wind that fills the entire space. We discuss the astrophysical implications of our results for jets from active galactic nuclei and for free pulsar winds such as that believed to power the Crab Nebula.

Begelman, Mitchell C.↗
NASA NTRS1999

Skeletal Structural Consequences of Reduced Gravity Environments

The overall goal of this project is to provide structurally meaningful data on bone loss after exposure to reduced gravity environments so that more precise estimates of fracture risk and the effectiveness of countermeasures in reducing fracture risk can be developed. The project has three major components: (1) measure structural changes in the limb bones of rats subjected to complete and partial nonweightbearing, with and without treatment with ibandronate and periodic full weightbearing; (2) measure structural changes in the limb bones of human bedrest subjects, with and without treatment with alendronate and resistive exercise, and Russian cosmonauts flying on the Mir Space Station; and (3) validate and extend the 2-dimensional structural analyses currently possible in the second project component (bedrest and Mir subjects) using 3-dimensional finite element modeling techniques, and determine actual fracture-producing loads on earth and in space.

Ruff, Christropher B.↗
NASA NTRS2005

Physics-based Simulation of Human Posture Using 3D Whole Body Scanning Technology for Astronaut Space Suit Evaluation

Over the past few years high precision three-dimensional (3D) full body laser scanners have been developed to be used as a powerful anthropometry tool for quantification of the morphology of the human body. The full body scanner can quickly extract body characteristics in non-contact fashion. It is required for the Anthropometry and Biomechanics Facility (ABF) to have capabilities for kinematics simulation of a digital human at various postures whereas the laser scanner only allows capturing a single static posture at each time. During this summer fellowship period a theoretical study has been conducted to estimate an arbitrary posture with a series of example postures through finite element (FE) approximation and found that four-point isoparametric FE approximation would result in reasonable maximum position errors less than 5%. Subsequent pilot scan experiments demonstrated that a bead marker with a nominal size of 6 mm could be used as a marker for digitizing 3-D coordinates of anatomical landmarks for further kinematic analysis. Two sessions of human subject testing were conducted for reconstruction of an arbitrary postures from a set of example postures for each joint motion for the forearm/hand complex and the whole upper extremity.

Kim, Kyu-Jung↗
DOE OSTI2025

BPS and near-BPS black holes in AdS 5 and their spectrum in $\mathcal{N}$ = 4 SYM

We study quantum corrections in the gravitational path integral around nearly 1/16-BPS black holes in asymptotically AdS 5 × S 5 space, dual to heavy states in 4D N = 4 super Yang-Mills. The analysis provides a gravitational explanation of why 1/16-BPS black holes exhibit an exact degeneracy at large N and why all such states have the same charges, confirming the belief that the superconformal index precisely counts the entropy of extremal black holes. We show the presence of a gap of order N −2 between the 1/16-BPS black holes and the lightest near-BPS black holes within the same charge sector. This is the first example of such a gap for black hole states within the context of AdS 5 holography. We also derive the spectrum of near-BPS states that lie above this gap. Our computation relies on finding the correct version of the N = 2 super-Schwarzian theory which captures the breaking of the SU(1, 1|1) symmetry when the black hole has finite temperature and non-zero chemical potential. Finally, we comment on possible stringy and non-perturbative corrections that can affect the black hole spectrum.

AdS-CFT Correspondence↗
DOE OSTI2026

Intrinsic nonlocality of spin- and polarization-resolved probabilities in strong-field quantum electrodynamics

Spin and polarization are central to precision tests of fundamental physics and for interpreting radiation from astrophysical sources and ultraintense laser-matter experiments. Here, focusing on the fundamental process of nonlinear Compton scattering, we demonstrate that a key assumption underlying current strong-field quantum electrodynamics models, i.e., that emission can be treated as an instantaneous random event sampled from a local differential rate, is inconsistent once emission angles, electron spin, and/or photon polarization are resolved. Namely, even in strictly constant and uniform fields , the resulting fully differential distribution is sign indefinite, yielding negative inferred probabilities. The physical reason is that the photon emission probability builds up over a finite length of the electron trajectory, the formation region, during which the electron direction changes by roughly the same small angle that defines the radiation cone. Therefore, we put forward a new method where we integrate over this formation region analytically to obtain a physically consistent electron spin and photon polarization model. We show that the implementation of our model is compatible with existing Monte Carlo and particle-in-cell workflows. Simulations of a GeV-class electron-laser collision accessible at current petawatt facilities and of emission in a pulsarlike magnetic field are shown to reveal spin and polarization patterns that differ even qualitatively from state-of-the-art local models. In particular, our new model predicts substantial angle-dependent circular photon polarization where the well-known collinear-emission approach yields none, and a pronounced helicity bias in the recoiling electrons absent from current predictions. These findings have direct implications for upcoming strong-field QED experiments and for interpreting polarized radiation from extreme astrophysical environments.

astrophysical electromagnetic fields↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Rotation and kinematics of the premain-sequence stars in Taurus-Auriga with Ca II emission

Radial velocities and v sin i values for the stars in the Taurus-Auriga region that were found to have strong Ca II H and K emission by Herbig, Vrba, and Rydgren 'HVR', (1986) are reported. Most of the velocities are determined to better than 2 km/s precision. The kinematic properties of the Ca II emission stars with strong Li are found to be indistinguishable from conventional T Tauris in Taurus-Auriga, contrary to HVR. These Li-rich stars also rotate like T Tauris. Most of the stars that lack Li are probable or possible members of the Hyades, in the foreground, and are among the brightest and most active stars in that cluster for their spectral types. It is suggested following Jones and Herbig (1979), that the apparent absence of low-mass stars older than 10 Myr in Taurus-Auriga is real, and is due to the finite lifetime of the cloud.

Hartmann, Lee W.↗
NASA NTRS2000

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗
DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗